[GRASS-dev] [GRASS GIS] #352: r.in.gdal region troubles in LL

[GRASS GIS]

#352: r.in.gdal region troubles in LL
-------------------------+--------------------------------------------------
 Reporter:  neteler      |       Owner:  grass-dev at lists.osgeo.org
     Type:  defect       |      Status:  new                      
 Priority:  major        |   Milestone:  6.4.0                    
Component:  default      |     Version:  unspecified              
 Keywords:               |    Platform:  Unspecified              
      Cpu:  Unspecified  |  
-------------------------+--------------------------------------------------
 The GeoTIFF file http://aloboaleu.googlepages.com/npp.tif (3.3MB) is
 causing problems in the region settings when imported with r.in.gdal:

 {{{
 gdalinfo npp.tif
 Driver: GTiff/GeoTIFF
 Files: npp.tif
 Size is 1440, 572
 Coordinate System is:
 GEOGCS["WGS 84",
     DATUM["WGS_1984",
         SPHEROID["WGS 84",6378137,298.2572235630016,
             AUTHORITY["EPSG","7030"]],
         AUTHORITY["EPSG","6326"]],
     PRIMEM["Greenwich",0],
     UNIT["degree",0.0174532925199433],
     AUTHORITY["EPSG","4326"]]
 Origin = (-180.000000000000000,85.000000000114397)
 Pixel Size = (0.250000000000200,-0.250000000000200)
 Metadata:
   AREA_OR_POINT=Area
 Image Structure Metadata:
   INTERLEAVE=BAND
 Corner Coordinates:
 Upper Left  (-180.0000000,  85.0000000) (180d 0'0.00"W, 85d 0'0.00"N)
 Lower Left  (-180.0000000, -58.0000000) (180d 0'0.00"W, 58d 0'0.00"S)
 Upper Right ( 180.0000000,  85.0000000) (180d 0'0.00"E, 85d 0'0.00"N)
 Lower Right ( 180.0000000, -58.0000000) (180d 0'0.00"E, 58d 0'0.00"S)
 Center      (   0.0000000,  13.5000000) (  0d 0'0.00"E, 13d30'0.00"N)
 Band 1 Block=1440x1 Type=Float32, ColorInterp=Gray

 g.region n=90 s=-90 w=-180 e=180 res=0:05 -p
 projection: 3 (Latitude-Longitude)
 zone:       0
 datum:      wgs84
 ellipsoid:  wgs84
 north:      90N
 south:      90S
 west:       180W
 east:       180E
 nsres:      0:05
 ewres:      0:05
 rows:       2160
 cols:       4320
 cells:      9331200

 r.in.gdal npp.tif out=npp
 Projection of input dataset and current location appear to match
  100%
 <npp> created
 r.in.gdal complete.

 g.region rast=npp -p
 projection: 3 (Latitude-Longitude)
 zone:       0
 datum:      wgs84
 ellipsoid:  wgs84
 north:      85N
 south:      58S
 west:       180W
 east:       179:59:59.999999W
 nsres:      0:15
 ewres:      0
 rows:       572
 cols:       1440
 cells:      823680

 g.region -p
 ERROR: region for current mapset is invalid
        line 9: <e-w resol: 0>
        run "g.region"
 }}}

 Apparently the global wrap-around recognition is failing for 180W/E.

 Markus

-- 
Ticket URL: <http://trac.osgeo.org/grass/ticket/352>
GRASS GIS <http://grass.osgeo.org>