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_03320000601_B10.TIF > L71203033_03320000601_B10.aux > L71203033_03320000601_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",500000], > PARAMETER["false_northing",0], > UNIT["metre",1, > AUTHORITY["EPSG","9001"]], > AUTHORITY["EPSG","32629"]] > Origin = (521999.999999999941792,4414200.000000000000000) > Pixel Size = (30.000000000000000,-30.000000000000000) > Metadata: > AREA_OR_POINT=Point > Image Structure Metadata: > INTERLEAVE=BAND > Corner Coordinates: > Upper Left ( 522000.000, 4414200.000) ( 8d44'33.81"W, 39d52'39.24"N) > Lower Left ( 522000.000, 4199970.000) ( 8d44'58.57"W, 37d56'49.39"N) > Upper Right ( 766830.000, 4414200.000) ( 5d52'54.42"W, 39d50'9.31"N) > Lower Right ( 766830.000, 4199970.000) ( 5d57'53.91"W, 37d54'29.39"N) > Center ( 644415.000, 4307085.000) ( 7d20'4.58"W, 38d54'2.83"N) > 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.66825833333333], > PARAMETER["central_meridian",-8.133108333333334], > 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.52"W, 39d52'38.78"N) > Lower Left ( -54187.411, -193166.260) ( 8d44'57.96"W, 37d55'35.73"N) > Upper Right ( 192732.267, 23417.315) ( 5d52'51.29"W, 39d51'26.27"N) > Lower Right ( 192732.267, -193166.260) ( 5d56'29.99"W, 37d54'28.07"N) > Center ( 69272.428, -84874.473) ( 7d20'4.40"W, 38d54'3.73"N) > 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",500000], >> PARAMETER["false_northing",0], >> UNIT["metre",1, >> AUTHORITY["EPSG","9001"]], >> AUTHORITY["EPSG","32636"]] >> Origin = (603487.500000000000000,3623062.500000000000000) >> Pixel Size = (28.500000000000000,-28.500000000000000) >> 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.43"E, 32d44'25.92"N) >> Lower Left ( 603487.500, 3407232.000) ( 34d 4'54.01"E, 30d47'37.09"N) >> Upper Right ( 847932.000, 3623062.500) ( 36d42'38.00"E, 32d41'25.66"N) >> Lower Right ( 847932.000, 3407232.000) ( 36d38'1.81"E, 30d44'49.94"N) >> Center ( 725709.750, 3515147.250) ( 35d22'58.18"E, 31d44'58.41"N) >> 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> 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> 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>> 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 >>>> >>>> >>>> >>> >>> This mail was received via Mail-SeCure System. >>> >>> >>> >>> -- >>> Micha Silver >>> Arava Development Co. +972-52-3665918 >>> http://surfaces.co.il >>> >>> >> >> This mail was received via Mail-SeCure System. >> >> >> _______________________________________________ >> grass-user mailing >> listgrass-u...@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user >> >> This mail was received via Mail-SeCure System. >> >> >> >> >> >> -- >> Micha Silverhttp://www.surfaces.co.il/ >> >> Arava Development Co. +972-52-3665918 >> >> > > This mail was received via Mail-SeCure System. > > > > -- > Micha Silverhttp://www.surfaces.co.il/ > Arava Development Co. +972-52-3665918 > >
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user