Re: [GRASS-user] Importing Landsat images with different coordinates
This all looks fine. And you said that the destination tiff was smaller than the original? Can you display the tiff files? Maybe it's some disk problem (in the log you sent last time, it seems there was an error reading the file...)? On 05/25/2010 04:10 PM, Jenny Turner wrote: Sorry Micha My original Landsat image: Driver: GTiff/GeoTIFF Files: L71203033_0332601_B10.TIF L71203033_0332601_B10.aux L71203033_0332601_B10.rrd Size is 8161, 7141 Coordinate System is: PROJCS[WGS 84 / UTM zone 29N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-9], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32629]] Origin = (521999.9941792,4414200.000) Pixel Size = (30.000,-30.000) Metadata: AREA_OR_POINT=Point Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 522000.000, 4414200.000) ( 8d44'33.81W, 39d52'39.24N) Lower Left ( 522000.000, 4199970.000) ( 8d44'58.57W, 37d56'49.39N) Upper Right ( 766830.000, 4414200.000) ( 5d52'54.42W, 39d50'9.31N) Lower Right ( 766830.000, 4199970.000) ( 5d57'53.91W, 37d54'29.39N) Center ( 644415.000, 4307085.000) ( 7d20'4.58W, 38d54'2.83N) Band 1 Block=8161x1 Type=Byte, ColorInterp=Gray Min=0.000 Max=255.000 Minimum=0.000, Maximum=255.000, Mean=60.374, StdDev=43.530 Overviews: 2041x1786, 1021x893, 511x447, 256x224, 128x112, 64x56 Metadata: STATISTICS_MINIMUM=0 STATISTICS_MAXIMUM=255 STATISTICS_MEAN=60.374122874202 STATISTICS_MEDIAN=8.5986839406512e-287 STATISTICS_MODE=1.7801025473311e-306 STATISTICS_STDDEV=43.529818998446 LAYER_TYPE=athematic And the obtained file: Driver: GTiff/GeoTIFF Files: destination.TIF Size is 8229, 7218 Coordinate System is: PROJCS[ETRS89 / Portugal TM06, GEOGCS[ETRS89, DATUM[European_Terrestrial_Reference_System_1989, SPHEROID[GRS 1980,6378137,298.2572221010042, AUTHORITY[EPSG,7019]], AUTHORITY[EPSG,6258]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4258]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,39.668258], PARAMETER[central_meridian,-8.1331084], PARAMETER[scale_factor,1], PARAMETER[false_easting,0], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,3763]] Origin = (-54187.411212536622770,23417.314857349840167) Pixel Size = (30.006036956002887,-30.006036956002887) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( -54187.411, 23417.315) ( 8d45'59.52W, 39d52'38.78N) Lower Left ( -54187.411, -193166.260) ( 8d44'57.96W, 37d55'35.73N) Upper Right ( 192732.267, 23417.315) ( 5d52'51.29W, 39d51'26.27N) Lower Right ( 192732.267, -193166.260) ( 5d56'29.99W, 37d54'28.07N) Center ( 69272.428, -84874.473) ( 7d20'4.40W, 38d54'3.73N) Band 1 Block=8229x1 Type=Byte, ColorInterp=Gray Hope this helps you (and mee :D ) Thank you Jenny 2010/5/25 Micha Silver mi...@arava.co.il mailto:mi...@arava.co.il On 05/25/2010 11:17 AM, Jenny Turner wrote: Yes Sure :) Attached to this email you can find the log I was thinking of just the results of gdalinfo. This way we can see if the tiff is being recognized, and if its CRS is is available in the tiff header. Here's an example: ~gdalinfo geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Driver: GTiff/GeoTIFF Files: geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Size is 8577, 7573 Coordinate System is: PROJCS[WGS 84 / UTM zone 36N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,33], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32636]] Origin = (603487.500,3623062.500) Pixel Size = (28.500,-28.500) Metadata: TIFFTAG_XRESOLUTION=72 TIFFTAG_YRESOLUTION=72 TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Re: [GRASS-user] Importing Landsat images with different coordinates
Micha You are absolutely right. I was running it in a shared folder, and for some reason (not space beucase it's not full), i was getting that error. Now, I tried in a Desktop folder and it worked. Thanks 2010/5/26 Micha Silver mi...@arava.co.il This all looks fine. And you said that the destination tiff was smaller than the original? Can you display the tiff files? Maybe it's some disk problem (in the log you sent last time, it seems there was an error reading the file...)? On 05/25/2010 04:10 PM, Jenny Turner wrote: Sorry Micha My original Landsat image: Driver: GTiff/GeoTIFF Files: L71203033_0332601_B10.TIF L71203033_0332601_B10.aux L71203033_0332601_B10.rrd Size is 8161, 7141 Coordinate System is: PROJCS[WGS 84 / UTM zone 29N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-9], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32629]] Origin = (521999.9941792,4414200.000) Pixel Size = (30.000,-30.000) Metadata: AREA_OR_POINT=Point Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 522000.000, 4414200.000) ( 8d44'33.81W, 39d52'39.24N) Lower Left ( 522000.000, 4199970.000) ( 8d44'58.57W, 37d56'49.39N) Upper Right ( 766830.000, 4414200.000) ( 5d52'54.42W, 39d50'9.31N) Lower Right ( 766830.000, 4199970.000) ( 5d57'53.91W, 37d54'29.39N) Center ( 644415.000, 4307085.000) ( 7d20'4.58W, 38d54'2.83N) Band 1 Block=8161x1 Type=Byte, ColorInterp=Gray Min=0.000 Max=255.000 Minimum=0.000, Maximum=255.000, Mean=60.374, StdDev=43.530 Overviews: 2041x1786, 1021x893, 511x447, 256x224, 128x112, 64x56 Metadata: STATISTICS_MINIMUM=0 STATISTICS_MAXIMUM=255 STATISTICS_MEAN=60.374122874202 STATISTICS_MEDIAN=8.5986839406512e-287 STATISTICS_MODE=1.7801025473311e-306 STATISTICS_STDDEV=43.529818998446 LAYER_TYPE=athematic And the obtained file: Driver: GTiff/GeoTIFF Files: destination.TIF Size is 8229, 7218 Coordinate System is: PROJCS[ETRS89 / Portugal TM06, GEOGCS[ETRS89, DATUM[European_Terrestrial_Reference_System_1989, SPHEROID[GRS 1980,6378137,298.2572221010042, AUTHORITY[EPSG,7019]], AUTHORITY[EPSG,6258]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4258]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,39.668258], PARAMETER[central_meridian,-8.1331084], PARAMETER[scale_factor,1], PARAMETER[false_easting,0], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,3763]] Origin = (-54187.411212536622770,23417.314857349840167) Pixel Size = (30.006036956002887,-30.006036956002887) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( -54187.411, 23417.315) ( 8d45'59.52W, 39d52'38.78N) Lower Left ( -54187.411, -193166.260) ( 8d44'57.96W, 37d55'35.73N) Upper Right ( 192732.267, 23417.315) ( 5d52'51.29W, 39d51'26.27N) Lower Right ( 192732.267, -193166.260) ( 5d56'29.99W, 37d54'28.07N) Center ( 69272.428, -84874.473) ( 7d20'4.40W, 38d54'3.73N) Band 1 Block=8229x1 Type=Byte, ColorInterp=Gray Hope this helps you (and mee :D ) Thank you Jenny 2010/5/25 Micha Silver mi...@arava.co.il On 05/25/2010 11:17 AM, Jenny Turner wrote: Yes Sure :) Attached to this email you can find the log I was thinking of just the results of gdalinfo. This way we can see if the tiff is being recognized, and if its CRS is is available in the tiff header. Here's an example: ~gdalinfo geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Driver: GTiff/GeoTIFF Files: geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Size is 8577, 7573 Coordinate System is: PROJCS[WGS 84 / UTM zone 36N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,33], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32636]] Origin =
Re: [GRASS-user] Importing Landsat images with different coordinates
On 05/25/2010 11:17 AM, Jenny Turner wrote: Yes Sure :) Attached to this email you can find the log I was thinking of just the results of gdalinfo. This way we can see if the tiff is being recognized, and if its CRS is is available in the tiff header. Here's an example: ~gdalinfo geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Driver: GTiff/GeoTIFF Files: geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Size is 8577, 7573 Coordinate System is: PROJCS[WGS 84 / UTM zone 36N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,33], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32636]] Origin = (603487.500,3623062.500) Pixel Size = (28.500,-28.500) Metadata: TIFFTAG_XRESOLUTION=72 TIFFTAG_YRESOLUTION=72 TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) AREA_OR_POINT=Point Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 603487.500, 3623062.500) ( 34d 6'16.43E, 32d44'25.92N) Lower Left ( 603487.500, 3407232.000) ( 34d 4'54.01E, 30d47'37.09N) Upper Right ( 847932.000, 3623062.500) ( 36d42'38.00E, 32d41'25.66N) Lower Right ( 847932.000, 3407232.000) ( 36d38'1.81E, 30d44'49.94N) Center ( 725709.750, 3515147.250) ( 35d22'58.18E, 31d44'58.41N) Band 1 Block=8577x1 Type=Byte, ColorInterp=Gray Thanks for your help On Mon, May 24, 2010 at 8:46 PM, Micha Silver mi...@arava.co.il mailto:mi...@arava.co.il wrote: On 05/24/2010 06:26 PM, Jenny Turner wrote: Hi Following this topic I have tested this methodology by converting from UTM 29N to a specific local coordinate system (EPSG:3763) (available at PROJ4) gdalwarp -t_srs EPSG:3763 band_landsat.TIF destination.tif With that I got a segmentation error, one TIFFReadDirectory a lot of TIFFWriteEncodedStrip and the destination.tif has a different size. Can you send the output of: gdalinfo band_landsat.TIF Any tips of what migh be wrong? Thank you Jenny On Tue, May 18, 2010 at 8:42 PM, Micha Silver mi...@arava.co.il mailto:mi...@arava.co.il wrote: Jenny Turner wrote: Greetings Firman, but imagine that I have rasters in UTM for differente zones (e.g. huge countries like USA, etc). How can I produce an automatic importing tool that uses gdalwarp to project an image to WGS84 and import to GRASS? Using gdalwarp you don't have to specify the source coordinate system, *if* the originals have the correct projection defined in their headers. For Landsat images, this is most likely the case, so you should be able to get away with just setting the -t_srs option of gdalwarp and each tile will be re-projected to the WGS84 geographic CRS. Thanks for your help Best regards, Jenny On Mon, May 17, 2010 at 11:10 PM, Firman Hadi jalmibur...@gmail.com mailto:jalmibur...@gmail.com mailto:jalmibur...@gmail.com mailto:jalmibur...@gmail.com wrote: Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org http://www.sigro.org http://www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org mailto:grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918
Re: [GRASS-user] Importing Landsat images with different coordinates
Sorry Micha My original Landsat image: Driver: GTiff/GeoTIFF Files: L71203033_0332601_B10.TIF L71203033_0332601_B10.aux L71203033_0332601_B10.rrd Size is 8161, 7141 Coordinate System is: PROJCS[WGS 84 / UTM zone 29N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-9], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32629]] Origin = (521999.9941792,4414200.000) Pixel Size = (30.000,-30.000) Metadata: AREA_OR_POINT=Point Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 522000.000, 4414200.000) ( 8d44'33.81W, 39d52'39.24N) Lower Left ( 522000.000, 4199970.000) ( 8d44'58.57W, 37d56'49.39N) Upper Right ( 766830.000, 4414200.000) ( 5d52'54.42W, 39d50'9.31N) Lower Right ( 766830.000, 4199970.000) ( 5d57'53.91W, 37d54'29.39N) Center ( 644415.000, 4307085.000) ( 7d20'4.58W, 38d54'2.83N) Band 1 Block=8161x1 Type=Byte, ColorInterp=Gray Min=0.000 Max=255.000 Minimum=0.000, Maximum=255.000, Mean=60.374, StdDev=43.530 Overviews: 2041x1786, 1021x893, 511x447, 256x224, 128x112, 64x56 Metadata: STATISTICS_MINIMUM=0 STATISTICS_MAXIMUM=255 STATISTICS_MEAN=60.374122874202 STATISTICS_MEDIAN=8.5986839406512e-287 STATISTICS_MODE=1.7801025473311e-306 STATISTICS_STDDEV=43.529818998446 LAYER_TYPE=athematic And the obtained file: Driver: GTiff/GeoTIFF Files: destination.TIF Size is 8229, 7218 Coordinate System is: PROJCS[ETRS89 / Portugal TM06, GEOGCS[ETRS89, DATUM[European_Terrestrial_Reference_System_1989, SPHEROID[GRS 1980,6378137,298.2572221010042, AUTHORITY[EPSG,7019]], AUTHORITY[EPSG,6258]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4258]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,39.668258], PARAMETER[central_meridian,-8.1331084], PARAMETER[scale_factor,1], PARAMETER[false_easting,0], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,3763]] Origin = (-54187.411212536622770,23417.314857349840167) Pixel Size = (30.006036956002887,-30.006036956002887) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( -54187.411, 23417.315) ( 8d45'59.52W, 39d52'38.78N) Lower Left ( -54187.411, -193166.260) ( 8d44'57.96W, 37d55'35.73N) Upper Right ( 192732.267, 23417.315) ( 5d52'51.29W, 39d51'26.27N) Lower Right ( 192732.267, -193166.260) ( 5d56'29.99W, 37d54'28.07N) Center ( 69272.428, -84874.473) ( 7d20'4.40W, 38d54'3.73N) Band 1 Block=8229x1 Type=Byte, ColorInterp=Gray Hope this helps you (and mee :D ) Thank you Jenny 2010/5/25 Micha Silver mi...@arava.co.il On 05/25/2010 11:17 AM, Jenny Turner wrote: Yes Sure :) Attached to this email you can find the log I was thinking of just the results of gdalinfo. This way we can see if the tiff is being recognized, and if its CRS is is available in the tiff header. Here's an example: ~gdalinfo geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Driver: GTiff/GeoTIFF Files: geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Size is 8577, 7573 Coordinate System is: PROJCS[WGS 84 / UTM zone 36N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,33], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32636]] Origin = (603487.500,3623062.500) Pixel Size = (28.500,-28.500) Metadata: TIFFTAG_XRESOLUTION=72 TIFFTAG_YRESOLUTION=72 TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) AREA_OR_POINT=Point Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 603487.500, 3623062.500) ( 34d 6'16.43E, 32d44'25.92N) Lower Left ( 603487.500, 3407232.000) ( 34d 4'54.01E, 30d47'37.09N) Upper Right ( 847932.000, 3623062.500) ( 36d42'38.00E, 32d41'25.66N) Lower Right ( 847932.000, 3407232.000) ( 36d38'1.81E, 30d44'49.94N) Center ( 725709.750, 3515147.250) (
Re: [GRASS-user] Importing Landsat images with different coordinates
Hi Following this topic I have tested this methodology by converting from UTM 29N to a specific local coordinate system (EPSG:3763) (available at PROJ4) gdalwarp -t_srs EPSG:3763 band_landsat.TIF destination.tif With that I got a segmentation error, one TIFFReadDirectory a lot of TIFFWriteEncodedStrip and the destination.tif has a different size. Any tips of what migh be wrong? Thank you Jenny On Tue, May 18, 2010 at 8:42 PM, Micha Silver mi...@arava.co.il wrote: Jenny Turner wrote: Greetings Firman, but imagine that I have rasters in UTM for differente zones (e.g. huge countries like USA, etc). How can I produce an automatic importing tool that uses gdalwarp to project an image to WGS84 and import to GRASS? Using gdalwarp you don't have to specify the source coordinate system, *if* the originals have the correct projection defined in their headers. For Landsat images, this is most likely the case, so you should be able to get away with just setting the -t_srs option of gdalwarp and each tile will be re-projected to the WGS84 geographic CRS. Thanks for your help Best regards, Jenny On Mon, May 17, 2010 at 11:10 PM, Firman Hadi jalmibur...@gmail.commailto: jalmibur...@gmail.com wrote: Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org http://www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918 http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
On 05/24/2010 06:26 PM, Jenny Turner wrote: Hi Following this topic I have tested this methodology by converting from UTM 29N to a specific local coordinate system (EPSG:3763) (available at PROJ4) gdalwarp -t_srs EPSG:3763 band_landsat.TIF destination.tif With that I got a segmentation error, one TIFFReadDirectory a lot of TIFFWriteEncodedStrip and the destination.tif has a different size. Can you send the output of: gdalinfo band_landsat.TIF Any tips of what migh be wrong? Thank you Jenny On Tue, May 18, 2010 at 8:42 PM, Micha Silver mi...@arava.co.il mailto:mi...@arava.co.il wrote: Jenny Turner wrote: Greetings Firman, but imagine that I have rasters in UTM for differente zones (e.g. huge countries like USA, etc). How can I produce an automatic importing tool that uses gdalwarp to project an image to WGS84 and import to GRASS? Using gdalwarp you don't have to specify the source coordinate system, *if* the originals have the correct projection defined in their headers. For Landsat images, this is most likely the case, so you should be able to get away with just setting the -t_srs option of gdalwarp and each tile will be re-projected to the WGS84 geographic CRS. Thanks for your help Best regards, Jenny On Mon, May 17, 2010 at 11:10 PM, Firman Hadi jalmibur...@gmail.com mailto:jalmibur...@gmail.com mailto:jalmibur...@gmail.com mailto:jalmibur...@gmail.com wrote: Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org http://www.sigro.org http://www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org mailto:grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918 http://www.surfaces.co.il This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918 http://surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Greetings Firman, but imagine that I have rasters in UTM for differente zones (e.g. huge countries like USA, etc). How can I produce an automatic importing tool that uses gdalwarp to project an image to WGS84 and import to GRASS? Thanks for your help Best regards, Jenny On Mon, May 17, 2010 at 11:10 PM, Firman Hadi jalmibur...@gmail.com wrote: Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Dear Jenny, Yes, I am agree that will be such a huge work :D. Actually I am not familiar with shell script but maybe I can imagine the workflow and perhaps someone in this list can help. I think the algorithm of such script is: 1. Put all the imageries in one folder 2. Use gdalinfo command to read coordinate system on each file (looping command) 3. Rename each file according to its coordinate system (EPSG code) 4. Use gdalwarp command to convert each file (with specific naming convention) to WGS84 I don't know this will help :D Regards, Firman. On 18 May 2010, at 15:29, Jenny Turner wrote: Greetings Firman, but imagine that I have rasters in UTM for differente zones (e.g. huge countries like USA, etc). How can I produce an automatic importing tool that uses gdalwarp to project an image to WGS84 and import to GRASS? Thanks for your help Best regards, Jenny On Mon, May 17, 2010 at 11:10 PM, Firman Hadi jalmibur...@gmail.com wrote: Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Jenny Turner wrote: Greetings Firman, but imagine that I have rasters in UTM for differente zones (e.g. huge countries like USA, etc). How can I produce an automatic importing tool that uses gdalwarp to project an image to WGS84 and import to GRASS? Using gdalwarp you don't have to specify the source coordinate system, *if* the originals have the correct projection defined in their headers. For Landsat images, this is most likely the case, so you should be able to get away with just setting the -t_srs option of gdalwarp and each tile will be re-projected to the WGS84 geographic CRS. Thanks for your help Best regards, Jenny On Mon, May 17, 2010 at 11:10 PM, Firman Hadi jalmibur...@gmail.com mailto:jalmibur...@gmail.com wrote: Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org http://www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918 http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Greetings Thanks It really worked but now I have one particualr question: I need to have an automatic proceeding but, in this case my input images are not in the same coordinate system (e.g. different UTM Zones). Can I apply this methdology without defining source_coordinate system? Thanks Jenny On Mon, May 10, 2010 at 12:21 PM, Firman Hadi jalmibur...@gmail.com wrote: On 10 May 2010, at 18:00, Jenny Turner wrote: Greetings I need to import Landsat images (r.in.gdal) to my location but, my Landsat images are in a specific Coordinate system (WGS84) that is different from mine. Could anyone explain me (briefly ) the steps to import this images? Thanks Jenny Dear Jenny, If your images are in TIFF, you can convert the coordinate system to the one that you have, with gdalwarp. As an example, if my coordinate system is UTM Zone 48S (32748) and my data is in lat-long (4326), I will use gdalwarp command like this: gdalwarp -t_srs EPSG:32748 -s_srs EPSG:4326 source.tif target.tif Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Jenny Turner wrote: I need to import Landsat images (r.in.gdal) to my location but, my Landsat images are in a specific Coordinate system (WGS84) that is different from mine. Could anyone explain me (briefly ) the steps to import this images? (There is also a within-grass solution) It is possible to import the Landsat images in a location using their coordinate system (it might be useful to use this location later to export material to the same source coordinate system). After importing them, reproject (r.proj) the images to the location of your interest (e.g. in a wgs84-defined location). Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
On 10 May 2010, at 18:00, Jenny Turner wrote: Greetings I need to import Landsat images (r.in.gdal) to my location but, my Landsat images are in a specific Coordinate system (WGS84) that is different from mine. Could anyone explain me (briefly ) the steps to import this images? Thanks Jenny Dear Jenny, If your images are in TIFF, you can convert the coordinate system to the one that you have, with gdalwarp. As an example, if my coordinate system is UTM Zone 48S (32748) and my data is in lat-long (4326), I will use gdalwarp command like this: gdalwarp -t_srs EPSG:32748 -s_srs EPSG:4326 source.tif target.tif Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
