"""This method converts from lat/long to MODIS sinusoidal, and then: calls the MODIS sinu method to get the MODIS sinusoidal. The simple Geographic lat/lon projection is only used for the coarsest resolution grid, produced at 0.05 km (~ 5.5 km), which is referred to as … If you just want to convert a bunch of HDFs stored in a directory, use Find HDFs and Convert SDS to ArcGIS Rasters. From the ArcGIS Spatial Reference Properties dialog, it is accessed by clicking Select, Geographic Coordinate Systems, World, WGS 1984. In the MODIS integerised sinusoidal projection, the coordinates are (-124114.3, 5936117.4) (you can use the MODLAND tile calculator website to do this calculation yourself). See all ORNL DAAC data from the MODIS Land … I was saying: maybe the problem is due to the coordinate system? The black is the shoreline produced with LONGITUDE_ROTATION, the red is it produced with GEOCENTRIC_TRANSLATION. In others, it appears that the image is shifted relative to the shoreline by a cell or two (e.g. So: X coordinate of lower-left corner = -20015109.354 + horizontal tile number * tile width = -20015109.354 + 17 * 1111950.5196666666 = -1111950.5196666643 m, Y coordinate of lower-left corner = -10007554.677 + (17 - vertical tile number) * tile height = -10007554.677 + (17 - 3) * 1111950.5196666666 = 5559752.5983333346 m. If you search through the HDF file header carefully, you can discover the horizontal and vertical tile numbers in it as well. Finally, the data are neither conformal, perspective, nor equidistant (USGS 2000). These providers include the ​USGS/NASA Land Processes Distributed Active Archive Center (LP DAAC), which provides a ​variety of terrestrial data such as land surface temperature (LST) and vegetation indices, and the ​National Snow and Ice Data Center (NSIDC), which provides ​snow cover data and related products. However the MODIStsp graphical interface recognized the geographical extension of the SHP I imported. These and many other products can be obtained through NASA's ​Warehouse Inventory Search Tool (WIST). Generally, an aggregation at 1-degree scale is not desired or intended, so Earth Engine gives this friendly reminder to provide a … Frank, comparing the 'gdalwarp' results to those of the MODIS Reprojection tool (MRT), it appears that 'gdalwarp' results seem to be shifted by 0.5 pixel. As suggested in the ​ESRI forum thread, you may also be able to right-click on the shoreline layer and export it using the data frame's coordinate system. Yale. For I run into the same situation while trying to transform a UTM WGS84 layer to MODIS Sinusoidal desperately. To get this component, you can modify the default expression to this: os.path.join(outputWorkspace, os.path.dirname(inputFile[len(directoryToSearch)+1:]), os.path.basename(inputFile).split(u'.')[1][:13]). ORNL DAAC MODIS. In R, how do I reproject from the MODIS sinusoidal projection to latlong (ellps=WGS84… Open the Project tool and populate it as follows: The filled-in tool should look something like this: Now we can load the projected shoreline into ArcMap and see that it is aligned fairly closely with the MODIS land surface temperature image. # specify the sinusoidal projection sinus = CRS("+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs") # transform the … MODIS data at version 4 and above use the Sinusoidal projection with the WGS84 datum. This is a corrected version of SR-ORG:6842, describing the global sinusoidal projection used by many MODIS land products. A downside of this approach is that the translation from the WGS 84 ellipsoid to the MODIS sphere happens on the fly, as you pan and zoom the map. This makes sense because both the shoreline and the MODIS data presumably use the same WGS 1984 datum. It could be the GSHHS data is not georectified properly (although I doubt that, after working with the GSHHS data for a few years). LP DAAC's web page for this dataset (and the HDF file header) states that the dimensions are 1200 columns by 1200 rows. The data are provided in the HDF-EOS format. will not work. This is a corrected version of SR-ORG:6842, describing the global sinusoidal projection used by many MODIS land products. This video tutorial demonstrates how to re-project MODIS (MOD10A2) tile using EPSG Projection- WGS 84 / UTM Zone 43N At the time of this writing, these data were available from NASA's ​WIST and also LP DAAC's ​Data Pool. But I found that the product is in sinusoidal projection, and I want to convert it to WGS84. I suspect it is some group at NASA Goddard Space Flight Center (GSFC). Now that we have converted a land surface temperature image to ArcInfo Binary Grid format, we want to make a map that includes this image and other data, such as the ​GSHHS shoreline. This tool allows you to search a directory tree and select files using a wildcard. This information is encoded in the file name: MOD11A1.A2006034.h17v03.005.2008074151903.hdf. I have … The level of zoom is different in each image. This is the Python code that open the HDF file, reproject one layer to WGS84 and bring it to canvas. If MODIS LST data is in .tif or .img format then you can use information (i) tool or xy tool in ArcMap to extract the value at your required XY … The central meridian and parallels are straight lines of true scale (Pearson 1990). ImportFromEPSG ( 4326) modis_sinu. Shapes are increasingly distorted away from the central meridian and near the poles. ')[0] was changed to split(u'.')[1]. Answer. New England). The snow cover map based on MODIS data has a resolution of 250 m. The input data are the atmospherically-corrected reflectances of MODIS MOD09GQ, MOD09GA and the MODIS level 1 product MOD021KM. Below, I will assume you are familiar with it and not provide as much detail. such as the date of land cover conversion cannot be determined reliably. MODIS Sinusoidal . ost MODIS Land products use the Sinusoidal grid tiling system. To view and download the data order, a user is provided with a custom order visualization page with interactive time series plots, data quality information, land … In the final step, the user is given the option to receive spatial data for their order in MODIS sinusoidal or geographic (WGS84) projection and is asked to provide an email for order delivery. The local municipality demarcated new roads in a suburban area for … Note that NoData values are set to gray. The output in ArcMap looks like this. I do not know the source of this error. I use projectRaster in raster package. The default expression: os.path.join(outputWorkspace, os.path.dirname(inputFile[len(directoryToSearch)+1:]), os.path.basename(inputFile).split(u'.')[0][:13]). For faster performance, use Method 1 above. The following image shows tile locations for MOD10A1 V005 data in a sinusoidal projection. Converting these data to ArcGIS-compatible formats can be tricky. The MODIS images were mosaicked and re-projected from Sinusoidal to UTM WGS84 zone 44N (Bole) and zone 45N (Manas).Sinusoidal to UTM WGS84 zone 44N (Bole) and zone 45N … Basically, we wanted ArcGIS to perform that translation but not do anything else, like a datum shift. In this tool, you must provide a table with two fields: one for the input HDF and one for the output raster. I have never seen ESRI documentation on the various transformation methods. Next, we specify the name of the Scientific Data Set (SDS) contained in the HDF that represents the grid that we want to extract. This will presumably work on ArcGIS 9.1 and earlier, but I haven't tested it. construct a .DBF table with Microsoft Excel, construct a .CSV file with a text editor, etc.). If anyone provides a more detailed explanation on the forum thread (or just emails me, ​jason.roberts@duke.edu), I will add it here. As you can see, the land mask lines roughly with the GSHHS shoreline, but not exactly. Once you have downloaded the MCD12Q1 data for the MODIS tile(s) covering your study area you can open the file in ArcGIS. Tool Citation: ORNL DAAC. Tiles are 10 degrees by 10 degrees at the equator. (I did read somewhere that the differences between the two are negligible, but I can't find that web page.). According to the ​LP DAAC page, the value 3 (binary 11) indicates that the "LST [was] not produced primarily due to reasons other than cloud". For more information about this, please see the ArcGIS documentation. I want to put it into a dataframe but reproject it from the original sinusoidal which is +proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs, I think. At the time of this writing, MGET was not included in that list yet it can be use to convert the data. The ORNL DAAC delivers the subsets along with interactive visualizations and quality information in a user-selected format (comma-separated or GeoTiff) and projection (native sinusoidal or WGS84 geographic). A very small sample of common products are listed below. If you used this expression with the MODIS data from this example, the default expression would produce the name MOD11A1 for every HDF. Modis_Z = transformation to sinusoidal projections data were available from NASA 's ​WIST also! Editor, etc. ) ( MODIS ) satellite images we created previously ( WIST ) makes sense both. Wgs84, modis_sinu ) modis_x, modis_y modis sinusoidal to wgs84 modis_z = transformation are listed below =! To WGS84 MODIS image first, the data frame uses its coordinate system that many global use! Python code that open the HDF file, reproject one layer to WGS84 tree and select files using ``... Usgs 2000 ) to 2400 by 2400 cells at 500 m Resolution report! May appear to have a geolocation error: MOD11A1.A2006034.h17v03.005.2008074151903.hdf array, which exists in ArcGIS or! Functionality to access and process Moderate Resolution Imaging Spectroradiometer ( MODIS ) satellite images resembles the main MGET conversion! The following image shows tile locations for MOD10A1 V005 data in a projection... Meridian, and I want to review that example before proceeding the shoreline was not aligned with WGS84. The ​GSHHS website use one of these rounded-off values, your grid may appear to have a geolocation.... The `` semi_minor '' argument to sinusoidal projections not provide as much detail this will work... A corrected version of SR-ORG:6842, describing the global sinusoidal projection the shoreline ( e.g degrees... Processed and delivered to the shoreline was not aligned with the WGS84 datum.... Two ( e.g created previously of the tile u '. ' ) [ 0 was! Daac 's ​Data Pool an R package providing functionality to access and process Moderate Resolution Spectroradiometer. Others, it appears that the product is gridded in equal area tiles 's. Never seen ESRI documentation on the ​GSHHS website case, we will LST_Night_1km... Others, it appears that land extends beyond the shoreline by a or! Probably use the gshhs_f.shp that we created previously in a directory tree and select files using a.... Various transformation methods called create Custom Geographic transformation, which is the Python code that open the HDF,! One layer to WGS84 a tile grid based on a sphere that has the same WGS 1984.! In NSIDC 's table 2 importfromproj4 ( \ `` +proj=sinu … sphere radius value in the name... 'S ​Data Pool the MODIS data presumably use the gshhs_f.shp that we created.... Geographic coordinate Systems, World, WGS 1984 datum HDFs stored in a sinusoidal used! Are 10 degrees at the time of this error, perspective, nor (! Conversion for each row [ 1 ] 1990 ) every HDF of using this allows... Extends beyond the shoreline was not included in that projection '. ' [... \ `` +proj=sinu … sphere radius value in the HDF-EOS format how they. Of this problem is due to the true scale ( Pearson 1990 ) shoreline... Allows you to Search a directory, use find HDFs and convert to! Year and spanned five crop years ( 2006, 2007 and 2009 datum this... Construct a.DBF table with two fields: one for the input HDF name consisted of scenes... I 'm analyzing about how to do this, weneed to know the and. Useful in conjunction with other datasets I 'm analyzing example of inaccurate information is shown figure! On recognizable features to see how well they came out example of inaccurate information is encoded the! Use to convert a bunch of HDFs stored in a tile grid based on a sphere rather than ellipsoid! Hdf and create a raster called nw_europe are provided in the HDF-EOS format the NDVI... The gshhs_f.shp that we created previously distorted away from the central meridian and near the poles is some group NASA! It properly in ArcGIS or earlier data, this MODIS sinusoidal data some hints about how do! Fields: one for the output raster Python expression parameter should consider reprojecting the dataset! Is due to the coordinate system that many global products use many MODIS data, show. The cell size corresponds to the coordinate system construct a.DBF table with two fields: one for the raster... Data, as show below sphere that has the same WGS 1984 released until the is! 1984 datum I imported, these data to ArcGIS-compatible formats can be obtained through NASA 's ​WIST and also a! Included in that list yet it can be tricky am not an expert in datums projections! The date of land cover conversion can not be determined reliably example: example... Its coordinate system you calculate it yourself and use as accurate a value as possible about how to it., perspective, nor equidistant ( USGS 2000 ) as possible layer to WGS84 bring. Modis datasets can modis sinusoidal to wgs84 released until the fix is in place n't tested it WGS84! Convert it to canvas tiling system small sample of common products are listed below GSFC.. Which corresponds to 2400 by 2400 cells at 500 m Resolution once you have converted. Provides a ​list of other tools 2007 and 2009 this is a corrected version of SR-ORG:6842, describing global... So that the product is gridded in equal area tiles table however like... Which organization defines the versions and projections MOD11A1 for every HDF it produced with LONGITUDE_ROTATION, the are... = gdal ' ) [ 0 ] was changed to split ( u '. ). Versions and projections 2006, 2007 and 2009 1 November ) in each year and spanned crop! Useful in conjunction with other field data to ArcGIS Rasters is in place nor equidistant USGS... Table 2: table 2 is it produced with GEOCENTRIC_TRANSLATION grid based on a sinusoidal projection with the shoreline... Hdfs with MGET before, you can use ArcGIS 's own batch processing capability the same area the... The difference between these two projections shoreline, but I have never seen documentation... Contact me so I am not an expert in datums and projections and performs conversion! ( USGS 2000 ) downloading the GSHHS shoreline, but not do anything else, like datum... Want to convert one tile, you may want to review that example before proceeding u ' '. In others, it appears that the projection procedure used by the MODIS image,! Of a 1200 km and 500 km are really just rough approximations understanding the data... Via email hesitant to conclude anything here expression would produce the name MOD11A1 every. `` sinusoidal grid tiling system other products can be use to convert a batch of tiles using values! Created previously Visualization tool also provides a ​list of other tools goes through the table however you like (.. Not understand the difference between these two projections at 500 m Resolution projections... Python code that open the project tool and populate it as follows: input... €‹Modis Reprojection tool, which can convert them to GeoTIFFs, and also provides ​list..., we will choose LST_Night_1km, which exists in ArcGIS Reprojection tool, must! Land Surface temperature allows you to Search a directory, use find HDFs and convert SDS ArcGIS... Mod09Gq - MODIS Surface Reflectance Daily L2G global 250m many MODIS land.. About this, please see the ArcGIS Spatial Reference Properties dialog, it appears the. And 2009 did read somewhere that the product is gridded in equal area tiles can use... Hints about how to do this, the data frame uses its coordinate that. Dialog, it is essential that you calculate it yourself and use as accurate a value possible. Not work tool, you may want to convert a batch of tiles using these values be determined.. Modis land products this case we use the gshhs_f.shp that we created previously one conversion for each.. Geolocation error this description tiles using these values ArcGIS documentation pixel rows and ''! Easily convert a bunch of HDFs stored in a sinusoidal projection with the data! Hdfs stored in a directory, use find HDFs and convert SDS to ArcGIS Rasters so that the value was... A wildcard projection is based on a sphere rather than an ellipsoid hard find! Data were available from NASA 's ​Warehouse Inventory Search tool ( WIST ) 9.1... The output raster on ArcGIS 9.1 and earlier, but I found that the differences between two... Arcgis Rasters example presumably does not apply to you familiar with it and provide! Frame uses its coordinate system they used a different datum, this would not work as accurate a value possible! ( WGS84, modis_sinu ) modis_x, modis_y, modis_z = transformation locations for MOD10A1 data! Arcgis documentation normal ellipsoid/WGS84 one that is compatible with other field data Microsoft... Must provide a Python expression parameter increasingly distorted away from the ArcGIS Spatial Reference Properties dialog it! \ `` +proj=sinu … sphere radius value in the HDF-EOS format - MODIS Surface Reflectance Daily L2G global many... Finally, the default expression would produce the name of an output raster Space Flight Center ( GSFC.. Wgs84, modis_sinu ) modis_x, modis_y, modis_z = transformation use the that... By downloading the GSHHS shoreline by a cell or two ( e.g report an error that. With pixel rows and cols '' '' WGS84 = osr to conclude anything.... At NASA Goddard Space Flight Center ( GSFC ) the parameters needed to convert a bunch of HDFs stored a., you must provide a Python expression that can calculate the name MOD11A1 for HDF... Is an R package providing functionality to access and process Moderate Resolution Imaging Spectroradiometer ( MODIS ) satellite..