[mapserver-users] WMS from netCDF returns white/transparent image
Hello, I am trying to generate WMS services from layers in a NetCDF file. WCS services already work well...While a WMS call does return an image (ex. PNG), it is all white or transparent, whereas a WCS call to the same layer returns an image with a range of pixel values. I have run gdalinfo on the layer to get the data range, and have tried using the data range to set the PROCESSING DIRECTIVE to scale the output (e.g. PROCESSING SCALE=311.76470947266 ,423.52941894531 ). I have also used the SCALE=AUTO directive, but both just return a white image. When I try using WMS on this layer in QGIS, it is correctly projected, but invisible (zoom to extents on the imported WMS layer is correct). Below is the output from gdalinfo to get the original data range, and then the mapfile entry for the layer to see how I have set this up gdalinfo NETCDF:flat_sample.nc:I_lw Band 1 Block=170x1 Type=Float32, ColorInterp=Undefined Min=311.765 Max=423.529 Minimum=311.765, Maximum=423.529, Mean=344.247, StdDev=16.599 NoData Value=9.96920996838686905e+36 Metadata: ipw_desc=incoming thermal (long-wave) radiation NETCDF_DIM_time=14 NETCDF_VARNAME=I_lw standard_name=downwelling_longwave_flux_in_air STATISTICS_MAXIMUM=423.52941894531 STATISTICS_MEAN=344.24709697293 STATISTICS_MINIMUM=311.76470947266 STATISTICS_STDDEV=16.598777681661 units=W m-2 LAYER DATA 'NETCDF:/geodata/watershed-data/8c/8c18e10f-44b0-4462-bf4f-1c27c984aece/flat_sample.nc:I_lw' EXTENT -116.142905556 43.729406 -116.13758 43.732697 METADATA wcs_label downwelling_longwave_flux_in_air layer_title downwelling_longwave_flux_in_air raster_selected yes wms_metadataurl_type FGDC-STD-001-1998 ows_keywordlist wcs_formats GEOTIFF_16 background no wcs_bands_description ISNoBal Timesteps 1-16 ows_abstract incoming thermal (long-wave) radiation wcs_rangeset_name Timesteps wms_metadataurl_href http://129.24.63.85//apps/vwp/datasets/05e73f4f-51d3-417c-9757-e376b5969a9c/metadata/FGDC-STD-001-1998.xml; wcs_bands_values 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 legend_display yes queryable no wcs_rangeset_label netCDF flat_sample.nc time_sensitive no ows_srs epsg:26911 wcs_bands_name Timestep wcs_bands_rangeitem _bands ows_title incoming thermal (long-wave) radiationstatic no wms_metadataurl_format text/xml wcs_bands_label Timestep Number wcs_rangeset_axes bands wms_encoding UTF-8 wcs_enable_request * layer_data_units W m-2 annotation_name downwelling_longwave_flux_in_air: 2015-06-19 20:30:29.202558 END # METADATANAME downwelling_longwave_flux_in_air PROCESSING CLOSE_CONNECTION=DEFER PROCESSING SCALE=311.76470947266 ,423.52941894531 #PROCESSING “SCALE=AUTO” PROJECTION init=epsg:26911 END # PROJECTION STATUS ON TILEITEM location TYPE RASTER UNITS METERS END # LAYER It seems I am missing something herewould appreciate any insight into how to return a valid image in WMS Thanks, Bill ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] WMS services from a NetCDF file layer
Hello, I am using mapserver to provide WCS and WMS services from individual layers (and their constituent bands) contained in a CF-formatted NetCDF file. I have successfully implemented WCS services, but am continuing to have problems with WMS services. While calls to return a single raster grid (derived from a NetCDF layer) as a WMS request returns a correctly projected image, the image is entirely white, and does not seem to correctly reinterpret original data values into a range of pixel color values. I would appreciate any information on what I might be missing. Thanks! -- ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] mapserver and netcdf subdatasets/bands
Hello all, I am trying to provide wcs and wms services from specific bands within a netcdf file. The documentation appears to be sparse online, and was wondering whether anyone can point me to some good information about how to set up a layer in the mapfile that points to a specific band in a netcdf file Thanks much ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] Return pixel values from raster layers served by mapserver
I guess I am a little unclear as to how to proceed. The PG temporal
index already works for WMS calls. I can pass a time stamp in the URL
request, and I get the correct image. Do I still need a template for
returning the x/y query? Or is there another way to return the values
via a script that sequentially runs through a list of times?
On 01/06/2015 09:37 AM, Lime, Steve D (MNIT) wrote:
If you want to use one x/y to query multiple rasters, that should be doable in
straight MapServer although I guess it depends on how things are organized. I
don't know how things work from a query perspective if you use your tile index
as just a spatial index and not a temporal index. I'd get the one layer version
working then expand it - make sure the PG tileindex works as expected with a
query.
-Original Message-
From: William Hudspeth [mailto:bhudsp...@edac.unm.edu]
Sent: Tuesday, January 06, 2015 10:02 AM
To: Lime, Steve D (MNIT); mapserver-users@lists.osgeo.org
Subject: Re: [mapserver-users] Return pixel values from raster layers served by
mapserver
I have not yet tried a TEMPLATEideally, I am looking for a method to script
a series of calls to mapserver to return a list of pixel values for a time
series of rasters for a given geographic coordinate. In this case, imagine a
time series of air pollution rasters where I want to collect a sequence of
values for a given city or location. I want to be able to use these data points
to build a time series graphAny ideas would be appreciated...
Thanks!
On 01/06/2015 08:44 AM, Lime, Steve D (MNIT) wrote:
Ok, so this just a PostGIS-based tiled layer to access GeoTIFF's. You should
just be able to do a normal MapServer query.
e.g.
http://your.server/cgi-bin/mapserv?map=your.mapmapxy=x+yqlayer=pm25_
monthlymode=query
You'd need to set a TEMPLATE for the pm25_monthly layer. Have you tried that?
WMS should work then too but it's simpler to debug using a basic MapServer
query.
Steve
-Original Message-
From: mapserver-users-boun...@lists.osgeo.org
[mailto:mapserver-users-boun...@lists.osgeo.org] On Behalf Of Bill
Hudspeth
Sent: Tuesday, January 06, 2015 8:59 AM
To: mapserver-users@lists.osgeo.org
Subject: [mapserver-users] Return pixel values from raster layers
served by mapserver
Hello,
I previously wrote about this issue, but have not found a solution yet. I am
running Mapserver 6.4.1, and am using Postgres/Postgis to display a time series
of geotiff images. I would like to know the best way for a user to submit a
given x,y coordinate pair and return the pixel value for a given image. Is this
possible using only the mapserver api? or must I absolutely use mapscript?
Currently, I use a PostgreSQL/PostGIS time-enabled table to store the
information about the location of the geotiffs and the timestamp associated
with them:
pm25_monthly=# \d pm25_monthly_global
Table
public.pm25_monthly_global
Column |Type |
Modifiers
--+-+-
--+-+--
--+-+
oid | integer | default nextval('pm25_monthly_seq'::regclass)
rundate | timestamp with time zone |
datetime| timestamp with time zone|
location| text|
timezone| timestamp without time zone |
the_geom | geometry(MultiPolygon,4326) |
Indexes:
pm25_monthly_global_oid_key
UNIQUE CONSTRAINT, btree (oid)
The 'location' field simply gives the full path to the relevant geotiff on the
file system. And, I define the index in the mapfile:
#**LAYER - TILE INDEX FOR PM2.5 RASTERS (Time sensitive) LAYER
STATUS ON
NAME pm25_time_index
TYPE POLYGON
DATA the_geom from pm25_monthly_global using unique oid using SRID=4326
METADATA
ows_title PM25 TIME INDEX
ows_srs EPSG:4326
ows_extent -180 -90 180 90 #[minx][miny][maxx][maxy]
wms_timeformat -MM-DDTHH:MM:SS
wms_timeextent 2011-01-01T00:00:00/2012-12-01T00:00:00
wms_timeitem timezone
wms_timedefault 2011-01-01T00:00:00
wcs_resolution '0:10 0:10'
END
#CONNECTION user=wilbur dbname=pm25_monthly host=localhost
password=ii
#CONNECTIONTYPE postgis
CONNECTIONTYPE POSTGIS
CONNECTION host=127.0.0.1 dbname=pm25_monthly user=wilbur password=EUR2_live
port=5432
DATA the_geom from pm25_monthly_global using unique oid using srid=4326
END
#***LAYER - MONTHLY PM2.5 (Time
Sensitive)*
LAYER
NAME 'pm25_monthly'
TYPE RASTER
STATUS ON
DEBUG OFF
DUMP TRUE
#TRANSPARENCY 50
#FILTER
METADATA
ows_keywordlist 'PM2.5 (ug/m3), Global' #wms_keywordlist
(Keywordlist)
ows_description Global monthly mean PM 2.5 (ug/m3)concentration
ows_name Global monthly mean
Re: [mapserver-users] Return pixel values from raster layers served by mapserver
I have not yet tried a TEMPLATEideally, I am looking for a method to
script a series of calls to mapserver to return a list of pixel values
for a time series of rasters for a given geographic coordinate. In this
case, imagine a time series of air pollution rasters where I want to
collect a sequence of values for a given city or location. I want to be
able to use these data points to build a time series graphAny ideas
would be appreciated...
Thanks!
On 01/06/2015 08:44 AM, Lime, Steve D (MNIT) wrote:
Ok, so this just a PostGIS-based tiled layer to access GeoTIFF's. You should
just be able to do a normal MapServer query.
e.g.
http://your.server/cgi-bin/mapserv?map=your.mapmapxy=x+yqlayer=pm25_monthlymode=query
You'd need to set a TEMPLATE for the pm25_monthly layer. Have you tried that?
WMS should work then too but it's simpler to debug using a basic MapServer
query.
Steve
-Original Message-
From: mapserver-users-boun...@lists.osgeo.org
[mailto:mapserver-users-boun...@lists.osgeo.org] On Behalf Of Bill Hudspeth
Sent: Tuesday, January 06, 2015 8:59 AM
To: mapserver-users@lists.osgeo.org
Subject: [mapserver-users] Return pixel values from raster layers served by
mapserver
Hello,
I previously wrote about this issue, but have not found a solution yet. I am
running Mapserver 6.4.1, and am using Postgres/Postgis to display a time series
of geotiff images. I would like to know the best way for a user to submit a
given x,y coordinate pair and return the pixel value for a given image. Is this
possible using only the mapserver api? or must I absolutely use mapscript?
Currently, I use a PostgreSQL/PostGIS time-enabled table to store the
information about the location of the geotiffs and the timestamp associated
with them:
pm25_monthly=# \d pm25_monthly_global
Table
public.pm25_monthly_global
Column |Type |
Modifiers
--+-+---
--+-+
oid | integer | default nextval('pm25_monthly_seq'::regclass)
rundate | timestamp with time zone |
datetime| timestamp with time zone|
location| text|
timezone| timestamp without time zone |
the_geom | geometry(MultiPolygon,4326) |
Indexes:
pm25_monthly_global_oid_key
UNIQUE CONSTRAINT, btree (oid)
The 'location' field simply gives the full path to the relevant geotiff on the
file system. And, I define the index in the mapfile:
#**LAYER - TILE INDEX FOR PM2.5 RASTERS (Time sensitive) LAYER
STATUS ON
NAME pm25_time_index
TYPE POLYGON
DATA the_geom from pm25_monthly_global using unique oid using SRID=4326
METADATA
ows_title PM25 TIME INDEX
ows_srs EPSG:4326
ows_extent -180 -90 180 90 #[minx][miny][maxx][maxy]
wms_timeformat -MM-DDTHH:MM:SS
wms_timeextent 2011-01-01T00:00:00/2012-12-01T00:00:00
wms_timeitem timezone
wms_timedefault 2011-01-01T00:00:00
wcs_resolution '0:10 0:10'
END
#CONNECTION user=wilbur dbname=pm25_monthly host=localhost
password=ii
#CONNECTIONTYPE postgis
CONNECTIONTYPE POSTGIS
CONNECTION host=127.0.0.1 dbname=pm25_monthly user=wilbur password=EUR2_live
port=5432
DATA the_geom from pm25_monthly_global using unique oid using srid=4326
END
#***LAYER - MONTHLY PM2.5 (Time Sensitive)*
LAYER
NAME 'pm25_monthly'
TYPE RASTER
STATUS ON
DEBUG OFF
DUMP TRUE
#TRANSPARENCY 50
#FILTER
METADATA
ows_keywordlist 'PM2.5 (ug/m3), Global' #wms_keywordlist
(Keywordlist)
ows_description Global monthly mean PM 2.5 (ug/m3)concentration
ows_name Global monthly mean PM2.5 - Concentration Values (ug/m3)
ows_label Global PM 2.5 (ug/m3)
ows_srs EPSG:4326 EPSG:900913
ows_extent -180 -90 180 90 #[minx][miny][maxx][maxy]
ows_formats GEOTIFF
ows_nativeformat 8-bit GeoTIF
ows_timeitem timezone
wms_onlineresource
http://localhost/cgi-bin/mapserv?map=/var/www/pm25_mapserver/public_html/pm25_ogc_services.map;
wms_abstract 'Global monthly mean PM 2.5 (ug/m3)concentration'
#wms_abstract(Abstract)
wms_opaque '0' #wms_opaque(layer)
wms_title Global monthly mean PM 2.5 (ug/m3)concentration
ows_timeformat -MM-DDTHH:MM:SSZ
ows_timeextent 2011-01-01T00:00:00Z/2012-12-01T00:00:00Z
ows_timedefault 2011-01-01T00:00:00Z
wms_enable_request *
wcs_resolution 0:10 0:10
wcs_bandcount 1
wcs_size 3600 1800
layer_name pm25_monthly
END #METADATA
PROJECTION
init=epsg:4326 #WGS84 Decimal Degrees
END #PROJECTION
TILEITEM location
TILEINDEX pm25_time_index
END #LAYER
Thanks
[mapserver-users] Extract raster point values with mapserver
Hello,
This is how Postgresql makes the data available. PostgresQL stores the
information about the location of the geotiffs and the timestamp
associated with them:
pm25_monthly=# \d pm25_monthly_global
Table public.pm25_monthly_global
Column |Type | Modifiers
--+-+---
oid | integer | default
nextval('pm25_monthly_seq'::regclass)
rundate | timestamp with time zone|
datetime| timestamp with time zone|
location | text|
timezone| timestamp without time zone |
the_geom | geometry(MultiPolygon,4326) |
Indexes:
pm25_monthly_global_oid_key UNIQUE CONSTRAINT, btree (oid)
The 'location' field simply gives the full path to the relevant geotiff
on the file system. And, I define the index in the mapfile:
#*LAYER - TILE INDEX FOR
PM2.5 RASTERS (Time sensitive) **
LAYER
STATUS ON
NAME pm25_time_index
TYPE POLYGON
DATA the_geom from pm25_monthly_global using unique oid using SRID=4326
METADATA
ows_title PM25 TIME INDEX
ows_srs EPSG:4326
ows_extent -180 -90 180 90 #[minx][miny][maxx][maxy]
wms_timeformat -MM-DDTHH:MM:SS
wms_timeextent 2011-01-01T00:00:00/2012-12-01T00:00:00
wms_timeitem timezone
wms_timedefault 2011-01-01T00:00:00
wcs_resolution '0:10 0:10'
#wcs_resolution '0.10 0.10'
END
#CONNECTION user=wilbur dbname=pm25_monthly host=localhost
password=ii
#CONNECTIONTYPE postgis
CONNECTIONTYPE POSTGIS
CONNECTION host=127.0.0.1 dbname=pm25_monthly user=wilbur
password=EUR2_live port=5432
DATA the_geom from pm25_monthly_global using unique oid using srid=4326
END
#*LAYER - MONTHLY PM2.5
(Time Sensitive)*
#http://localhost/cgi-bin/mapserv?map=/var/www/pm25_mapserver/public_html/pm25_ogc_services.mapSERVICE=WMSVERSION=1.1.1REQUEST=GetMapLAYERS=pm25_monthlySTYLES=SRS=EPSG:4326BBOX=-180,-90,180,90WIDTH=1600HEIGHT=800FORMAT=image/pngSLD=http://localhost/pm25_mapserver/public_html/pm25_sld.xmlTIME=2011-09-01T00:00:00Z
LAYER
NAME 'pm25_monthly'
TYPE RASTER
STATUS ON
DEBUG OFF
DUMP TRUE
#TRANSPARENCY 50
#FILTER
METADATA
ows_keywordlist 'PM2.5 (ug/m3), Global' #wms_keywordlist
(Keywordlist)
ows_description Global monthly mean PM 2.5 (ug/m3)concentration
ows_name Global monthly mean PM2.5 - Concentration Values
(ug/m3)
ows_label Global PM 2.5 (ug/m3)
ows_srs EPSG:4326 EPSG:900913
ows_extent -180 -90 180 90 #[minx][miny][maxx][maxy]
ows_formats GEOTIFF
ows_nativeformat 8-bit GeoTIF
ows_timeitem timezone
wms_onlineresource
http://localhost/cgi-bin/mapserv?map=/var/www/pm25_mapserver/public_html/pm25_ogc_services.map;
wms_abstract 'Global monthly mean PM 2.5
(ug/m3)concentration' #wms_abstract(Abstract)
wms_opaque '0' #wms_opaque(layer)
wms_title Global monthly mean PM 2.5 (ug/m3)concentration
ows_timeformat -MM-DDTHH:MM:SSZ
ows_timeextent 2011-01-01T00:00:00Z/2012-12-01T00:00:00Z
ows_timedefault 2011-01-01T00:00:00Z
wms_enable_request *
wcs_resolution 0:10 0:10
#wcs_resolution 0.10 0.10
wcs_bandcount 1
wcs_size 3600 1800
layer_name pm25_monthly
END #METADATA
PROJECTION
init=epsg:4326 #WGS84 Decimal Degrees
END #PROJECTION
TILEITEM location
TILEINDEX pm25_time_index
END #LAYER
Thanks
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] Extract raster point values with mapserver
Hello, I am running Mapserver 6.4.1, and am using Postgres/Postgis to display a time series of geotiff images. I would like to know the best way for a user to submit a given x,y coordinate pair and return the pixel value for a given image. Is this possible using only the mapserver api? or must I absolutely use mapscript? Suggestions and examples welcome, Thanks much ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] Mapserv executable and path on Ubuntu linux
Hello, I have written my first mapfile, and I am trying to view a layer. I am using a time-enabled raster layer indexed in Postgis. I am a little confused about building the correct mapserv URL and referencing my mapfile. I am running Ubuntu linux 14.01, and trying to view my map from localhost. I have a path like this where I have the mapfile: /var/www/pm25/public_html The permissions on the mapfile are 755 with the owner myself, and the group is www-data. I am not sure how to construct a URL to point to mapserv and the mapfile to see a map Thanks -- William Hudspeth, PhD Analyst-Programmer III Earth Data Analysis Center MSC01 1110, Bandelier Hall West 212 1 University of New Mexico Albuquerque, NM 87131-0001 Tel: (505) 277-3622, ext. 249 Email: bhudsp...@edac.unm.edu ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] Problem with cascading WMS reprojection
Hello, I am trying to reproject a cascaded NASA MODIS NDVI raster. The original raster is served as epsg:4326 from the NASA servers. It is using version=1.3.0 and CRS-CRS:84. I am cascading it through my own mapserver client/server and trying to reproject as epsg:26913 (utm zone13), where the MAP projection and extents are in 26913. My MAP projection looks like this: MAP . EXTENT 93549 3464749 694099 4065299 UNITS METERS . END #END MAP OBJECT PROJECTION init=epsg:26913 #UTM-Zone 13N END #END MAP PROJECTION OBJECT and, my LAYER definition looks like this: LAYER NAME MOD13A2_E_NDVI TYPE RASTER CONNECTION http://neowms.sci.gsfc.nasa.gov/wms/wms?; CONNECTIONTYPE WMS STATUS ON METADATA wms_nameMOD13A2_E_NDVI wms_server_version 1.1.1 wms_srs EPSG:4326 wms_format image/jpeg wms_latlonboundingbox -180 -90 180 90 wms_transparent FALSE END #METADATA PROJECTION init=epsg:4326 #WGS84 Decimal Degrees END #PROJECTION END #LAYER Interestingly, even though the MAP projection is NOT set to EPSG:4326, I can make a successful request using epsg:4326. In all cases, the LAYER metadata uses EPSG:4326. So, while this request is successful: http://edacwms.unm.edu/cgi-bin/mapfiles/epht_environment/epht_enviro_wms?version=1.1.1SERVICE=WMSREQUEST=GetMapFORMAT=image/jpegSTYLES=SRS=EPSG:4326BBOX=-120.000,26.000,-97.000,44.000Width=600Height=470Layers=MOD13A2_E_NDVI This one is NOT successful: http://edacwms.unm.edu/cgi-bin/mapfiles/epht_environment/epht_enviro_wms?version=1.1.1SERVICE=WMSREQUEST=GetMapFORMAT=image/jpegSTYLES=SRS=EPSG:26913BBOX=93549,3464749,694099,4065299Width=600Height=470Layers=MOD13A2_E_NDVI where I get an error: msWMSLoadGetMapParams(): WMS server error. Invalid SRS given : SRS must be valid for all requested layers. Any insight would be greatly appreciated... Bill ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
VS: [mapserver-users] Problem with cascading WMS reprojection
Thank you Jukka, Your solution worked! ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] Problems with time-enabled WMS (WMS-T) in a cascaded layer
Hello, I am trying to make time-enabled WMS-T requests to a cascaded layer that originates from the NASA MODIS services. I can successfully make calls to directly to the NASA service as below: http://neowms.sci.gsfc.nasa.gov/wms/wms?version=1.1.1service=WMSrequest=GetMapSRS=EPSG:4326FORMAT=image/jpegSTYLES=BBox=-120.000,26.000,-97.000,44.000Width=600Height=470Layers=MOD13A2_E_NDVITIME=2000-02-18 ...other applicable dates include 2000-05-09, 2000-03-22, and 2011-02-02... The NASA getCapabilities URL for this layer is: http://neo.sci.gsfc.nasa.gov/FGDCMetadata?datasetId=MOD13A2_E_NDVI and the NASA getCapabilites URL for the entire WMS service is: http://neowms.sci.gsfc.nasa.gov/wms/wms?version=1.3.0service=WMSrequest=GetCapabilities When I try to pass a TIME argument to the cascaded layer, I only get the image for the default time as specified in the NASA getCapabilities call to the WMS service My cascaded URL looks like this: http://edacwms.unm.edu/cgi-bin/mapfiles/epht_environment/epht_enviro_wms?version=1.1.1SERVICE=WMSREQUEST=GetMapFORMAT=image/jpegSTYLES=SRS=EPSG:26913BBOX=-1011825.41,2963772.82,1141463.33,4903085.27Width=600Height=470Layers=MOD13A2_E_NDVITIME=2000-05-09 and the LAYER definition looks like this: LAYER NAME MOD13A2_E_NDVI TYPE RASTER CONNECTION http://neowms.sci.gsfc.nasa.gov/wms/wms?; CONNECTIONTYPE WMS STATUS ON METADATA wms_name MOD13A2_E_NDVI wms_server_version1.1.1 wms_srs EPSG:4326 wms_abstract Vegetation Index NDVI 16 day - Terra/MODIS wms_title Vegetation Index NDVI 16 day - Terra/MODIS #wms_time 2011-02-02 wms_formatimage/jpeg wms_connectiontimeout 60 wms_latlonboundingbox -180 -90 180 90 wms_transparent FALSE wms_timeextent 2000-02-18,2000-03-06,2000-03-22,2000-04-07,2000-04-23,2000-05-09,2000-05-25,2000-06-10,2000-06-26,2000-07-12,2000-07-28,2000-08-13,2000-08-29,2000-09-14,2000-09-30,2000-10-16,2000-11-01,2000-11-17,2000-12-03,2000-12-19,2001-01-01,2001-01-17,2001-02-02,2001-02-18,2001-03-06,2001-03-22,2001-04-07,2001-04-23,2001-05-09,2001-05-25,2001-06-10,2001-06-26,2001-07-12,2001-07-28,2001-08-13,2001-08-29,2001-09-14,2001-09-30,2001-10-16,2001-11-01,2001-11-17,2001-12-03,2001-12-19,2002-01-01,2002-01-17,2002-02-02,2002-02-18,2002-03-06,2002-03-22,2002-04-07,2002-04-23,2002-05-09,2002-05-25,2002-06-10,2002-06-26,2002-07-12,2002-07-28,2002-08-13,2002-08-29,2002-09-14,2002-09-30,2002-10-16,2002-11-01,2002-11-17,2002-12-03,2002-12-19,2003-01-01,2003-01-17,2003-02-02,2003-02-18,2003-03-06,2003-03-22,2003-04-07,2003-04-23,2003-05-09,2003-05-25,2003-06-10,2003-06-26,2003-07-12,2003-07-28,2003-08-13,2003-08-29,2003-09-14,2003-09-30,2003-10-16,2003-11-01,2003-11-17,2003-12-03,2003-12-19,2004-01-01,2004-01-17,2004-02-02,2004-02-18,2004-03-06,2004-03-22,2004-04-07,2004-04-23,2004-05-09,2004-05-25,2004-06-10,2004-06-26,2004-07-12,2004-07-28,2004-08-13,2004-08-29,2004-09-14,2004-09-30,2004-10-16,2004-11-01,2004-11-17,2004-12-03,2004-12-19,2005-01-01,2005-01-17,2005-02-02,2005-02-18,2005-03-06,2005-03-22,2005-04-07,2005-04-23,2005-05-09,2005-05-25,2005-06-10,2005-06-26,2005-07-12,2005-07-28,2005-08-13,2005-08-29,2005-09-14,2005-09-30,2005-10-16,2005-11-01,2005-11-17,2005-12-03,2005-12-19,2006-01-01,2006-01-17,2006-02-02,2006-02-18,2006-03-06,2006-03-22,2006-04-07,2006-04-23,2006-05-09,2006-05-25,2006-06-10,2006-06-26,2006-07-12,2006-07-28,2006-08-13,2006-08-29,2006-09-14,2006-09-30,2006-10-16,2006-11-01,2006-11-17,2006-12-03,2006-12-19,2007-01-01,2007-01-17,2007-02-02,2007-02-18,2007-03-06,2007-03-22,2007-04-07,2007-04-23,2007-05-09,2007-05-25,2007-06-10,2007-06-26,2007-07-12,2007-07-28,2007-08-13,2007-08-29,2007-09-14,2007-09-30,2007-10-16,2007-11-01,2007-11-17,2007-12-03,2007-12-19,2008-01-01,2008-01-17,2008-02-02,2008-02-18,2008-03-06,2008-03-22,2008-04-07,2008-04-23,2008-05-09,2008-05-25,2008-06-10,2008-06-26,2008-07-12,2008-07-28,2008-08-13,2008-08-29,2008-09-14,2008-09-30,2008-10-16,2008-11-01,2008-11-17,2008-12-03,2008-12-19,2009-01-01,2009-01-17,2009-02-02,2009-02-18,2009-03-06,2009-03-22,2009-04-07,2009-04-23,2009-05-09,2009-05-25,2009-06-10,2009-06-26,2009-07-12,2009-07-28,2009-08-13,2009-08-29,2009-09-14,2009-09-30,2009-10-16,2009-11-01,2009-11-17,2009-12-03,2009-12-19,2010-01-01,2010-01-17,2010-02-02,2010-02-18,2010-03-06,2010-03-22,2010-04-07,2010-04-23,2010-05-09,2010-05-25,2010-06-10,2010-06-26,2010-07-12,2010-07-28,2010-08-13,2010-08-29,2010-09-14,2010-09-30,2010-10-16,2010-11-01,2010-11-17,2010-12-03,2010-12-19,2011-01-01,2011-01-17,2011-02-02 #wms_timeextent 2000-02-18,2000-03-06/2000-12-19/P16D,2001-01-01/2001-12-19/P16D,2002-01-01/2002-12-19/P16D,2003-01-01/2003-12-19/P16D,2004-01-01/2004-02-18/P16D,2004-03-06/2004-12-19/P16D,2005-01-01/2005-12-19/P16D,2006-01-01/2006-12-19/P16D,2007-01-01/2007-12-19/P16D,2008-01-01/2008-02-18/P16D,2008-03-06/2008-12-19/P16D,2009-01-01/2009-12-19/P16D,2010-01-01/2010-12-19/P16D,2011-01-01/2011-02-02/P16D wms_timeitem TIME
[mapserver-users] Projection issues with cascading WMS requests
Hello, I am trying to cascade a NASA dataset into my local WMS server. Requests to the NASA dataset, as in the example below, require an CRS designation of CRS=CRS:84 http://neowms.sci.gsfc.nasa.gov/wms/wms?version=1.3.0service=WMSrequest=GetMapCRS=CRS:84FORMAT=image/jpegSTYLES=BBox=-180.000,-90.000,180.000,90.000Width=600Height=300Layers=MOD13A2_E_NDVITIME=2000-04-07 If I try to access the cascaded layer with the following WMS request: http://edacwms.unm.edu/cgi-bin/mapfiles/epht_environment/epht_enviro_wms?WMTVER=1.1.1SERVICE=WMSREQUEST=GetMapFORMAT=image/jpegSTYLES=SRS=EPSG:4326BBOX=-90.0,-180.0,90,180.0Width=600Height=300Layers=MODIS_16d_ndviTIME=2000-04-07 (also tried changing the order of BBOX params to --- BBOX=-180.0,-90.0,180.0,90.0) I get an error: msWMSLoadGetMapParams(): WMS server error. Invalid SRS given : SRS must be valid for all requested layers. My layer is defined as: LAYER NAME 'MODIS_16d_ndvi' TYPE RASTER STATUS ON DEBUG OFF DUMP TRUE #TRANSPARENCY 50 #FILTER METADATA ows_keywordlist 'Terra/MODIS, Vegetation Index [NDVI],16 day ' #wms_keywordlist (Keywordlist) ows_description Vegetation Index [NDVI] (16 day - Terra/MODIS) ows_name Vegetation Index [NDVI] (16 day - Terra/MODIS) ows_label Terra/MODIS 16 day NDVI ows_srs CRS:84 ows_extent -180 -90 180 90 #[minx][miny][maxx][maxy] ows_formats GEOTIFF ows_nativeformat 8-bit GeoTIF ows_timeitem timezone wms_abstract 'Vegetation Index [NDVI] (16 day - Terra/MODIS)' #wms_abstract(Abstract) wms_opaque '0' #wms_opaque(layer) wms_title Vegetation Index [NDVI] (16 day - Terra/MODIS) ows_timeformat -MM-DD ows_timeextent 2000-02-18/2011-02-02 ows_timeitem timezone ows_timedefault 2011-01-01 wcs_resolution 0:6 0:6 wcs_bandcount 1 #wcs_size 2700 319 layer_name MODIS_16d_ndvi END #METADATA PROJECTION init=epsg:4326 #WGS84 Decimal Degrees END #PROJECTION END #LAYER Any help appreciated on how to get around the projection issues... Thanks, Bill ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
