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

Atif_Rajan · June 27, 2005, 2:31pm

Dear Reader,

I am attempting overlay a vector map on a raster GEO TIFF
On setting the region ( g.region rast=2405-3.1) and then displaying the raster
map, the raster comes up fine on the graphics display.
However, attempting to display the vector map (d.vect c2405_3),
the following message appears:
" The bounding box of the map outside current region, nothing displayed"

Any words of wisdon as to what I am doing incorrectly or not doing for that
matter?

Thank you.

Atif Rajan

Maciej_Sieczka3 · June 27, 2005, 7:17pm

Atif

I guess that your vector layer doesn't fit your raster layer. Try "g.region vect=c2405_3 -p" and "g.region rast=2405-3.1 -p" to check whether the region extent in both cases is the same.

Maciek

----- Original Message ----- From: "Atif Rajan" <atif@lynxseismicdata.com>
To: <GRASSLIST@baylor.edu>
Sent: Monday, June 27, 2005 4:31 PM
Subject: [GRASSLIST:7322] The bounding box of the map outside current region

Dear Reader,

I am attempting overlay a vector map on a raster GEO TIFF
On setting the region ( g.region rast=2405-3.1) and then displaying the raster
map, the raster comes up fine on the graphics display.
However, attempting to display the vector map (d.vect c2405_3),
the following message appears:
" The bounding box of the map outside current region, nothing displayed"

Any words of wisdon as to what I am doing incorrectly or not doing for that
matter?

Thank you.

Atif Rajan

Maciej_Sieczka3 · June 27, 2005, 8:39pm

Atif

Plaese 'Reply to All' so that the message goes to list archive and other
folks can also help or make use of our conversation.

Hello Maciek, thank you. You are correc. They do not match.
How do I go about matching the region extent?

Are you sure they cover the same area? Are they in the same projection? How
did you get your layers into Grass? What do the v.info and r.info say?

Maciek Sieczka wrote:

Atif

I guess that your vector layer doesn't fit your raster layer. Try
"g.region vect=c2405_3 -p" and "g.region rast=2405-3.1 -p" to check
whether the region extent in both cases is the same.

Maciek

----- Original Message ----- From: "Atif Rajan"

Dear Reader,

I am attempting overlay a vector map on a raster GEO TIFF
On setting the region ( g.region rast=2405-3.1) and then displaying the
raster
map, the raster comes up fine on the graphics display.
However, attempting to display the vector map (d.vect c2405_3),
the following message appears:
" The bounding box of the map outside current region, nothing displayed"

Any words of wisdon as to what I am doing incorrectly or not doing for
that
matter?

Thank you.

Atif Rajan

Atif_Rajan · June 28, 2005, 3:19pm

Hello Maciek, thank you. You are correct. They do not match.
As a FYI, I imported the GEOTIFF (r.in.gdal) with the -o flag and I imported the shape file for the same map (v.in.ogr) with the -o flag.
I was unable to import either one without the -o flag
The region set in the grass database is far larger than that encompassed by this map.
How do I go about matching the region extent?
Thank you.
Atif

Maciek Sieczka wrote:

Atif

I guess that your vector layer doesn't fit your raster layer. Try "g.region vect=c2405_3 -p" and "g.region rast=2405-3.1 -p" to check whether the region extent in both cases is the same.

Maciek

----- Original Message ----- From: "Atif Rajan" <atif@lynxseismicdata.com>
To: <GRASSLIST@baylor.edu>
Sent: Monday, June 27, 2005 4:31 PM
Subject: [GRASSLIST:7322] The bounding box of the map outside current region

Dear Reader,

I am attempting overlay a vector map on a raster GEO TIFF
On setting the region ( g.region rast=2405-3.1) and then displaying the raster
map, the raster comes up fine on the graphics display.
However, attempting to display the vector map (d.vect c2405_3),
the following message appears:
" The bounding box of the map outside current region, nothing displayed"

Any words of wisdon as to what I am doing incorrectly or not doing for that
matter?

Thank you.

Atif Rajan

Maciej_Sieczka3 · June 28, 2005, 5:41pm

From: "Atif Rajan" <atif@lynxseismicdata.com>

Hello Maciek, thank you. You are correct. They do not match.

What doesn't match? Projections, regions, else?

Again, there is information missing:

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

You have to be specific if you want to be helped.

As a FYI, I imported the GEOTIFF (r.in.gdal) with the -o flag and I
imported the shape file for the same map (v.in.ogr) with the -o flag.

I was unable to import either one without the -o flag

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

The region set in the grass database is far larger than that encompassed
by this map.
How do I go about matching the region extent?
Thank you.
Atif

Maciek Sieczka wrote:

Atif

I guess that your vector layer doesn't fit your raster layer. Try
"g.region vect=c2405_3 -p" and "g.region rast=2405-3.1 -p" to check
whether the region extent in both cases is the same.

Maciek

----- Original Message ----- From: "Atif Rajan"

Dear Reader,

I am attempting overlay a vector map on a raster GEO TIFF
On setting the region ( g.region rast=2405-3.1) and then displaying the
raster
map, the raster comes up fine on the graphics display.
However, attempting to display the vector map (d.vect c2405_3),
the following message appears:
" The bounding box of the map outside current region, nothing displayed"

Any words of wisdon as to what I am doing incorrectly or not doing for
that
matter?

Thank you.

Atif Rajan

Ed_Davison1 · October 3, 2005, 8:03pm

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?

|----------------------------------------------------------------------------|
| 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