[GRASSLIST:8469] Re: The bounding box of the map outside current region

[Ed Davison]

I am having a similar problem.  I used r.in.gdal to import a GeoTIFF and
v.in.ogr to import a shp file.  The error comes when I display the radar
raster and then try to display the gshhs_shp imported vector line file.
The location I am working on is LatLong.

> Are you sure they cover the same area? Are they in the same projection? What 
> do the v.info and r.info say?

GRASS 6.0.0 (global):~/gis/data/gshhs > v.info map=gshhs_shp
+----------------------------------------------------------------------------+
 | Layer:    gshhs_shp                      Organization:
|
 | Mapset:   bfdi533                        Source Date:
|
 | Location: global                         Name of creator:
|
 | Database: /usr/local/storage/gisdata
|
 | Title:
|
 | Map Scale:  1:1
|
 | Map format: native
|

|----------------------------------------------------------------------------|
 |   Type of Map:  Vector (level: 2)
|
 |
|
 |   Number of points:       0               Number of areas:      0
|
 |   Number of lines:        219136          Number of islands:    0
|
 |   Number of boundaries:   0               Number of faces:      0
|
 |   Number of centroids:    0               Number of kernels:    0
|
 |
|
 |   Map is 3D:              0
|
 |   Number of dblinks:      1
|
 |
|
 |   Projection: Latitude-Longitude (zone 0)
|
 |            N: 83.633        S: -78.570
|
 |            E: 360.000        W: 0.000
|
 |            B: 0.000         T: 0.000
|
 |
|
 |   Digitize threshold: 0.00000
|
 |   Comments:
|
 |
|

+----------------------------------------------------------------------------+

GRASS 6.0.0 (global):~/gis/data/gshhs > r.info map=radar2761913
+----------------------------------------------------------------------------+
 | Layer:    radar2761913                   Date: Mon Oct  3 14:30:22
2005    |
 | Mapset:   bfdi533                        Login of Creator: bfdi533
|
 | Location: global
|
 | DataBase: /usr/local/storage/gisdata
|
 | Title:     ( radar2761913 )
|

|----------------------------------------------------------------------------|
 |
|
 |   Type of Map:  raster              Number of Categories: 255
|
 |   Data Type:    CELL
|
 |   Rows:         1600
|
 |   Columns:      3400
|
 |   Total Cells:  5440000
|
 |        Projection: Latitude-Longitude (zone 0)
|
 |            N:  50:22:48N    S:     21:36N   Res: 0:01:04.755
|
 |            E:  66:24:36W    W: 127:34:12W   Res: 0:01:04.757647
|
 |   Range of data:    min =  0 max = 255
|
 |
|
 |   Data Source:
|
 |
|
 |
|
 |
|
 |   Data Description:
|
 |    generated by r.in.gdal
|
 |
|
 |
|

+----------------------------------------------------------------------------+

> That *may* mean they are in different projections. What do "gdalinfo
> your_geotiff" and "ogrinfo -so -al your_shapefile" say?

GRASS 6.0.0 (global):~/gis/data/gshhs >
gdalinfo ../noaa/Radar2761913.tif Driver: GTiff/GeoTIFF
Size is 3400, 1600
Coordinate System is:
GEOGCS["NAD83",
    DATUM["North_American_Datum_1983",
        SPHEROID["GRS 1980",6378137,298.2572221010042,
            AUTHORITY["EPSG","7019"]],
        AUTHORITY["EPSG","6269"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4269"]]
Origin = (-127.570000,50.380000)
Pixel Size = (0.01798824,-0.01798750)
Metadata:
  TIFFTAG_XRESOLUTION=200
  TIFFTAG_YRESOLUTION=200
  TIFFTAG_RESOLUTIONUNIT=1 (unitless)
Corner Coordinates:
Upper Left  (-127.5700000,  50.3800000) (127d34'12.00"W, 50d22'48.00"N)
Lower Left  (-127.5700000,  21.6000000) (127d34'12.00"W, 21d36'0.00"N)
Upper Right ( -66.4100000,  50.3800000) ( 66d24'36.00"W, 50d22'48.00"N)
Lower Right ( -66.4100000,  21.6000000) ( 66d24'36.00"W, 21d36'0.00"N)
Center      ( -96.9900000,  35.9900000) ( 96d59'24.00"W, 35d59'24.00"N)
Band 1 Block=3400x2 Type=Byte, ColorInterp=Gray

GRASS 6.0.0 (global):~/gis/data/gshhs > ogrinfo -so -al gshhs.shp
INFO: Open of `gshhs.shp'
using driver `ESRI Shapefile' successful.

Layer name: gshhs
Geometry: Line String
Feature Count: 219136
Extent: (0.000000, -78.570396) - (360.000000, 83.633400)
Layer SRS WKT:
GEOGCS["GCS_WGS_1984",
    DATUM["WGS_1984",
        SPHEROID["WGS_1984",6378137.0,298.257223563]],
    PRIMEM["Greenwich",0.0],
    UNIT["Degree",0.0174532925199433]]
FNODE_: Integer (11.0)
TNODE_: Integer (11.0)
LPOLY_: Integer (11.0)
RPOLY_: Integer (11.0)
LENGTH: Real (13.3)
GSHHS_: Integer (11.0)
GSHHS_ID: Integer (11.0)

Any help is appreciated.

-- 

--------------------------
Ed Davison