[OSRS-PROJ] PROJ4 and ISIN (MODIS) wish

[Markus Neteler]

Hi,

the MODIS data from TERRA mission are delivered in a sort
of unusual projection:
 Integerized Sinusoidal (ISIN) projection

"The ISIN projection is unique to the MODIS land products, analogous to the
Sinusoidal projection except that the ISIN projection is centered about 0°
longitude and special coefficients are used to flatten the ellipsoid."

MODIS Reprojection Tool Distribution
http://edc.usgs.gov/programs/sddm/modisdist/index.shtml

I would like to know if there is a chance to add ISIN to PROJ4.

Background:

 Today I managed to compile a tool to import MODIS data into GRASS. However,
 this is currently done in image coordinates. The ISIN boundary coordinates
 of the MODIS maps are delivered, so it would be great to be able to use
 r.proj of GRASS (which utilizes an older version of PROJ4) to reproject
 the MODIS/ISIN to another projection.

 Above mentioned MODIS Reprojection Tool is generally working and writing
 for example UTM in GeoTIFF, but not for the interesting 250m MODIS data
 products. For these data another tool is available which needs IDL (I
 don't want to buy a license, sorry).

 As workaround I got "select_vNsds" running on my linux box which generates
 binary file from HDF (convenient for r.in.bin). Only remaining problem is
 that I cannot set up an ISIN projection in GRASS due to the PROJ4
 limitations to directly import MODIS into an ISIN GRASS location.

Above mentioned MODIS Reprojection Tool is coming with a modified
GTCP package where ISIN is added:

/******************************************************************************
NAME                           ISINFOR.C
 
PURPOSE:    Integerized Sinusoidal Library Functions - library routines to
            perform mapping to and from the Integerized Sinusoidal.  These
            functions perform the mapping from geographic coordinates
            (longitude/latitude) to the map projection coordinates (x/y).

/******************************************************************************
NAME                           ISININV.C
 
PURPOSE:    Integerized Sinusoidal Library Functions - library routines to
            perform mapping to and from the Integerized Sinusoidal.  These
            functions perform the mapping from projection coordinates (x/y)
            to the geographic coordinates (longitude/latitude).

As the tool is provided by USGS it should be public domain. So the code
adapted might be allowed without reimplementing everything (?).
Then GRASS would be one of the first GIS supporting ISIN :-)
Unfortunately I am not able to "hack" the PROJ4 to support ISIN...

Of course I will publish the GRASS module for MODIS (and also ASTER) import,
when things are working. 

Comments are welcome,

 Markus Neteler

----------------------------------------
PROJ.4 Discussion List
See http://www.remotesensing.org/proj for subscription, unsubscription
and other information.