/* * Author : Gérald FENOY * * Copyright 2010-2011 Fondazione Edmund Mach. All rights reserved. * * Permission is hereby granted, free of charge, to any person obtaining a copy * of this software and associated documentation files (the "Software"), to deal * in the Software without restriction, including without limitation the rights * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell * copies of the Software, and to permit persons to whom the Software is * furnished to do so, subject to the following conditions: * * The above copyright notice and this permission notice shall be included in * all copies or substantial portions of the Software. * * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN * THE SOFTWARE. */ /** * Cross platform definition of layerObj->class */ #ifndef WIN32 #define CLASS class #else #define CLASS _class #endif #include "service_internal_ms.h" #include "server_internal.h" #include "response_print.h" /** * Get a list of configuration keys having a corresponding mandatory ows_*. * Map composed by a main.cfg maps name as key and the corresponding * MapServer Mafile Metadata name to use * see doc from here : * - http://mapserver.org/ogc/wms_server.html * - http://mapserver.org/ogc/wfs_server.html * - http://mapserver.org/ogc/wcs_server.html * * @return a new map containing a table linking a name of a configuration key * to a corresponding mandatory ows_* keyword (ie. "fees" => "ows_fees"). */ map* getCorrespondance(){ map* res=createMap("encoding","ows_encoding"); addToMap(res,"abstract","ows_abstract"); addToMap(res,"title","ows_title"); addToMap(res,"keywords","ows_keywordlist"); addToMap(res,"fees","ows_fees"); addToMap(res,"accessConstraints","ows_accessconstraints"); addToMap(res,"providerName","ows_attribution_title"); addToMap(res,"providerSite","ows_service_onlineresource"); addToMap(res,"individualName","ows_contactperson"); addToMap(res,"positionName","ows_contactposition"); addToMap(res,"providerName","ows_contactorganization"); addToMap(res,"role","ows_role"); addToMap(res,"addressType","ows_addresstype"); addToMap(res,"addressCity","ows_city"); addToMap(res,"addressDeliveryPoint","ows_address"); addToMap(res,"addressPostalCode","ows_postcode"); addToMap(res,"addressAdministrativeArea","ows_stateorprovince"); addToMap(res,"addressCountry","ows_country"); addToMap(res,"phoneVoice","ows_contactvoicetelephone"); addToMap(res,"phoneFacsimile","ows_contactfacsimiletelephone"); addToMap(res,"addressElectronicMailAddress","ows_contactelectronicmailaddress"); // Missing Madatory Informations addToMap(res,"hoursOfService","ows_hoursofservice"); addToMap(res,"contactInstructions","ows_contactinstructions"); return res; } /** * Add width and height keys to an output maps containing the maximum width * and height for displaying the full data extent. * Restriction to an image having a size of 640x480 (width * height) * * @param output * @param minx the lower left x coordinate * @param miny the lower left y coordinate * @param maxx the upper right x coordinate * @param maxy the upper right y coordinate */ void setMapSize(maps* output,double minx,double miny,double maxx,double maxy){ double maxWidth=640; double maxHeight=480; double deltaX=maxx-minx; double deltaY=maxy-miny; double qWidth; qWidth=maxWidth/deltaX; double qHeight; qHeight=maxHeight/deltaY; #ifdef DEBUGMS fprintf(stderr,"deltaX : %.15f \ndeltaY : %.15f\n",deltaX,deltaY); fprintf(stderr,"qWidth : %.15f \nqHeight : %.15f\n",qWidth,qHeight); #endif double width=deltaX*qWidth; double height=height=deltaY*qWidth; if(deltaXcontent,"width",sWidth); addToMap(output->content,"height",sHeight); } } /** * Add a Reference key to an output containing the WMFS/WFS/WCS request for * accessing service result * * @param m the conf maps containing the main.cfg settings * @param tmpI the specific output maps to add the Reference key */ void setReferenceUrl(maps* m,maps* tmpI){ outputMapfile(m,tmpI); map *msUrl=getMapFromMaps(m,"main","mapserverAddress"); if(msUrl==NULL){ errorException (m, _("Unable to find any mapserverAddress defined in the main.cfg file"), "InternalError", NULL); exit(-1); } map *msOgcVersion=getMapFromMaps(m,"main","msOgcVersion"); map *dataPath=getMapFromMaps(m,"main","dataPath"); map *sid=getMapFromMaps(m,"lenv","usid"); map* format=getMap(tmpI->content,"mimeType"); map* rformat=getMap(tmpI->content,"requestedMimeType"); map* width=getMap(tmpI->content,"width"); map* height=getMap(tmpI->content,"height"); map* protoMap=getMap(tmpI->content,"msOgc"); map* versionMap=getMap(tmpI->content,"msOgcVersion"); char options[3][5][25]={ {"WMS","1.3.0","GetMap","layers=%s","wms_extent"}, {"WFS","1.1.0","GetFeature","typename=%s","wcs_extent"}, {"WCS","1.1.0","GetCoverage","coverage=%s","wcs_extent"} }; int proto=0; if(rformat==NULL){ rformat=getMap(tmpI->content,"mimeType"); } if(strncasecmp(rformat->value,"text/xml",8)==0) proto=1; if(strncasecmp(rformat->value,"image/tiff",10)==0) proto=2; if(protoMap!=NULL){ if(strncasecmp(protoMap->value,"WMS",3)==0) proto=0; else{ if(strncasecmp(protoMap->value,"WFS",3)==0) proto=1; else proto=2; } } char *protoVersion=options[proto][1]; if(proto==1){ if(msOgcVersion!=NULL) protoVersion=msOgcVersion->value; if(versionMap!=NULL) protoVersion=versionMap->value; } map* extent=getMap(tmpI->content,options[proto][4]); map* crs=getMap(tmpI->content,"crs"); int hasCRS=1; if(crs==NULL){ crs=getMapFromMaps(m,"main","crs"); if(crs==NULL){ crs=createMap("crs","epsg:4326"); hasCRS=0; } } char layers[128]; sprintf(layers,options[proto][3],tmpI->name); char* webService_url=(char*)malloc((strlen(msUrl->value)+strlen(format->value)+strlen(tmpI->name)+strlen(width->value)+strlen(height->value)+strlen(extent->value)+256)*sizeof(char)); if(proto>0){ sprintf(webService_url, "%s?map=%s/%s_%s.map&request=%s&service=%s&version=%s&%s&format=%s&bbox=%s&crs=%s", msUrl->value, dataPath->value, tmpI->name, sid->value, options[proto][2], options[proto][0], protoVersion, layers, rformat->value, extent->value, crs->value ); } else{ sprintf(webService_url, "%s?map=%s/%s_%s.map&request=%s&service=%s&version=%s&%s&width=%s&height=%s&format=%s&bbox=%s&crs=%s", msUrl->value, dataPath->value, tmpI->name, sid->value, options[proto][2], options[proto][0], protoVersion, layers, width->value, height->value, rformat->value, extent->value, crs->value ); } if(hasCRS==0){ freeMap(&crs); free(crs); } addToMap(tmpI->content,"Reference",webService_url); free(webService_url); } /** * Set projection for a layer in a MAPFILE using Authority Code and Name if * available or fallback to proj4 definition if available or fallback to * default EPSG:4326 * * @param output the output maps * @param m the opened mapObj * @param myLayer the layerObj * @param pszProjection a char* containing the SRS definition in WKT format */ void setSrsInformations(maps* output,mapObj* m,layerObj* myLayer, char* pszProjection){ OGRSpatialReferenceH hSRS; map* msSrs=NULL; hSRS = OSRNewSpatialReference(NULL); if( pszProjection!=NULL && strlen(pszProjection)>1){ if(OSRImportFromWkt( hSRS, &pszProjection ) == CE_None ){ char *proj4Str=NULL; if(OSRGetAuthorityName(hSRS,NULL)!=NULL && OSRGetAuthorityCode(hSRS,NULL)!=NULL){ char tmpSrs[20]; sprintf(tmpSrs,"%s:%s", OSRGetAuthorityName(hSRS,NULL),OSRGetAuthorityCode(hSRS,NULL)); msLoadProjectionStringEPSG(&m->projection,tmpSrs); msLoadProjectionStringEPSG(&myLayer->projection,tmpSrs); char tmpSrss[256]; sprintf(tmpSrss,"EPSG:4326 EPSG:900913 %s",tmpSrs); msInsertHashTable(&(m->web.metadata), "ows_srs", tmpSrss); msInsertHashTable(&(myLayer->metadata), "ows_srs", tmpSrss); #ifdef DEBUGMS fprintf(stderr,"isGeo %b\n\n",OSRIsGeographic(hSRS)==TRUE); #endif if(output!=NULL){ if(OSRIsGeographic(hSRS)==TRUE) addToMap(output->content,"crs_isGeographic","true"); else addToMap(output->content,"crs_isGeographic","false"); addToMap(output->content,"crs",tmpSrs); } } else{ OSRExportToProj4(hSRS,&proj4Str); if(proj4Str!=NULL && strlen(proj4Str)>0){ #ifdef DEBUGMS fprintf(stderr,"PROJ (%s)\n",proj4Str); #endif msLoadProjectionString(&(m->projection),proj4Str); msLoadProjectionString(&(myLayer->projection),proj4Str); if(output!=NULL){ if(OSRIsGeographic(hSRS)==TRUE) addToMap(output->content,"crs_isGeographic","true"); else addToMap(output->content,"crs_isGeographic","false"); } free(proj4Str); } else{ msLoadProjectionStringEPSG(&m->projection,"EPSG:4326"); msLoadProjectionStringEPSG(&myLayer->projection,"EPSG:4326"); if(output!=NULL){ addToMap(output->content,"crs_isGeographic","true"); } } if(output!=NULL){ addToMap(output->content,"crs","EPSG:4326"); addToMap(output->content,"real_extent","true"); } msInsertHashTable(&(m->web.metadata),"ows_srs", "EPSG:4326 EPSG:900913"); msInsertHashTable(&(myLayer->metadata),"ows_srs","EPSG:4326 EPSG:900913"); } } } else{ if(output!=NULL){ msSrs=getMap(output->content,"msSrs"); } if(msSrs!=NULL){ msLoadProjectionStringEPSG(&m->projection,msSrs->value); msLoadProjectionStringEPSG(&myLayer->projection,msSrs->value); char tmpSrs[128]; sprintf(tmpSrs,"%s EPSG:4326 EPSG:900913",msSrs->value); msInsertHashTable(&(m->web.metadata),"ows_srs",tmpSrs); msInsertHashTable(&(myLayer->metadata),"ows_srs",tmpSrs); }else{ msLoadProjectionStringEPSG(&m->projection,"EPSG:4326"); msLoadProjectionStringEPSG(&myLayer->projection,"EPSG:4326"); msInsertHashTable(&(m->web.metadata),"ows_srs","EPSG:4326 EPSG:900913"); msInsertHashTable(&(myLayer->metadata),"ows_srs","EPSG:4326 EPSG:900913"); } if(output!=NULL){ addToMap(output->content,"crs",msSrs->value); addToMap(output->content,"crs_isGeographic","true"); } } OSRDestroySpatialReference( hSRS ); } /** * Set the MAPFILE extent, the the ows_extent for the layer, add wms_extent and * wfs_extent to the output maps and call setMapSize. * * @param output the specific output * @param m the mapObj * @param myLayer the layerObj * @param minX the lower left x coordinate * @param minY the lower left y coordinate * @param maxX the upper right x coordinate * @param maxY the upper right y coordinate * @see setMapSize */ void setMsExtent(maps* output,mapObj* m,layerObj* myLayer, double minX,double minY,double maxX,double maxY){ msMapSetExtent(m,minX,minY,maxX,maxY); #ifdef DEBUGMS fprintf(stderr,"Extent %.15f %.15f %.15f %.15f\n",minX,minY,maxX,maxY); #endif char tmpExtent[1024]; sprintf(tmpExtent,"%.15f %.15f %.15f %.15f",minX,minY,maxX,maxY); #ifdef DEBUGMS fprintf(stderr,"Extent %s\n",tmpExtent); #endif msInsertHashTable(&(myLayer->metadata), "ows_extent", tmpExtent); if(output!=NULL){ map* test=getMap(output->content,"real_extent"); if(test!=NULL){ pointObj min, max; projectionObj tempSrs; min.x = m->extent.minx; min.y = m->extent.miny; max.x = m->extent.maxx; max.y = m->extent.maxy; char tmpSrsStr[1024]; msInitProjection(&tempSrs); msLoadProjectionStringEPSG(&tempSrs,"EPSG:4326"); msProjectPoint(&(m->projection),&tempSrs,&min); msProjectPoint(&m->projection,&tempSrs,&max); sprintf(tmpExtent,"%.3f,%.3f,%.3f,%.3f",min.y,min.x,max.y,max.x); map* isGeo=getMap(output->content,"crs_isGeographic"); #ifdef DEBUGMS fprintf(stderr,"isGeo = %s\n",isGeo->value); #endif if(isGeo!=NULL && strcasecmp("true",isGeo->value)==0) sprintf(tmpExtent,"%f,%f,%f,%f", minY,minX, maxY, maxX); addToMap(output->content,"wms_extent",tmpExtent); sprintf(tmpSrsStr,"%.3f,%.3f,%.3f,%.3f",min.x,min.y,max.x,max.y); addToMap(output->content,"wcs_extent",tmpExtent); }else{ sprintf(tmpExtent,"%f,%f,%f,%f",minX, minY, maxX, maxY); map* isGeo=getMap(output->content,"crs_isGeographic"); if(isGeo!=NULL){ #ifdef DEBUGMS fprintf(stderr,"isGeo = %s\n",isGeo->value); #endif if(isGeo!=NULL && strcasecmp("true",isGeo->value)==0) sprintf(tmpExtent,"%f,%f,%f,%f", minY,minX, maxY, maxX); } addToMap(output->content,"wms_extent",tmpExtent); sprintf(tmpExtent,"%.3f,%.3f,%.3f,%.3f",minX,minY,maxX,maxY); addToMap(output->content,"wcs_extent",tmpExtent); } } setMapSize(output,minX,minY,maxX,maxY); } /** * Try to open a vector output and define the corresponding layer in the MAPFILE * * @param conf the conf maps containing the main.cfg settings * @param output the specific output maps * @param m the mapObj */ int tryOgr(maps* conf,maps* output,mapObj* m){ map* tmpMap=getMap(output->content,"storage"); char *pszDataSource=tmpMap->value; /** * Try to open the DataSource using OGR */ OGRRegisterAll(); /** * Try to load the file as ZIP */ OGRDataSourceH poDS1 = NULL; OGRSFDriverH *poDriver1 = NULL; char *dsName=(char*)malloc((8+strlen(pszDataSource)+1)*sizeof(char)); char *odsName=zStrdup(pszDataSource); char *sdsName=zStrdup(pszDataSource); char *demo=".data"; sdsName[strlen(sdsName)-(strlen(demo)-1)]='d'; sdsName[strlen(sdsName)-(strlen(demo)-2)]='i'; sdsName[strlen(sdsName)-(strlen(demo)-3)]='r'; sdsName[strlen(sdsName)-(strlen(demo)-4)]=0; odsName[strlen(odsName)-(strlen(demo)-1)]='z'; odsName[strlen(odsName)-(strlen(demo)-2)]='i'; odsName[strlen(odsName)-(strlen(demo)-3)]='p'; odsName[strlen(odsName)-(strlen(demo)-4)]=0; sprintf(dsName,"/vsizip/%s",odsName); #ifdef DEBUGMS fprintf(stderr,"Try loading %s, %s, %s\n",dsName,odsName,dsName); #endif FILE* file = fopen(pszDataSource, "rb"); FILE* fileZ = fopen(odsName, "wb"); free(odsName); fseek(file, 0, SEEK_END); unsigned long fileLen=ftell(file); fseek(file, 0, SEEK_SET); char *buffer=(char *)malloc(fileLen+1); fread(buffer, fileLen, 1, file); fwrite(buffer,fileLen, 1, fileZ); fclose(file); fclose(fileZ); free(buffer); #ifdef DEBUGMS fprintf(stderr,"Try loading %s",dsName); #endif poDS1 = OGROpen( dsName, FALSE, poDriver1 ); if( poDS1 == NULL ){ fprintf(stderr,"Unable to access the DataSource as ZIP File\n"); setMapInMaps(conf,"lenv","message","Unable to open datasource in read only mode"); OGR_DS_Destroy(poDS1); }else{ #ifdef DEBUGMS fprintf(stderr,"The DataSource is a ZIP File\n"); #endif char** demo=VSIReadDir(dsName); int i=0; zMkdir(sdsName #ifndef WIN32 ,S_IRWXU | S_IRGRP | S_IXGRP | S_IROTH | S_IXOTH #endif ); while(demo[i]!=NULL){ #ifdef DEBUGMS fprintf(stderr,"ZIP File content : %s\n",demo[i]); #endif char *tmpDs=(char*)malloc((strlen(dsName)+strlen(demo[i])+2)*sizeof(char)); sprintf(tmpDs,"%s/%s",dsName,demo[i]); fprintf(stderr,"read : %s\n",tmpDs); VSILFILE* vsif=VSIFOpenL(tmpDs,"rb"); #ifdef DEBUGMS fprintf(stderr,"open : %s\n",tmpDs); #endif VSIFSeekL(vsif,0,SEEK_END); vsi_l_offset size=VSIFTellL(vsif); #ifdef DEBUGMS fprintf(stderr,"size : %d\n",size); #endif VSIFSeekL(vsif,0,SEEK_SET); char *vsifcontent=(char*) malloc(((int)size+1)*sizeof(char)); VSIFReadL(vsifcontent,1,(size_t)size,vsif); char *fpath=(char*) malloc((strlen(sdsName)+strlen(demo[1])+2)*sizeof(char)); sprintf(fpath,"%s/%s",sdsName,demo[i]); int f=zOpen(fpath,O_WRONLY|O_CREAT,S_IRUSR|S_IWUSR|S_IRGRP|S_IWGRP|S_IROTH|S_IWOTH); zWrite(f,vsifcontent,(int)size); close(f); chmod(fpath,S_IRWXU | S_IRGRP | S_IXGRP | S_IROTH); char* tmpP=strstr(fpath,".shp"); if(tmpP==NULL) tmpP=strstr(fpath,".SHP"); if(tmpP!=NULL){ #ifdef DEBUGMS fprintf(stderr,"*** DEBUG %s\n",strstr(tmpP,".")); #endif if( strcmp(tmpP,".shp")==0 || strcmp(tmpP,".SHP")==0 ){ tmpMap=getMap(output->content,"storage"); free(tmpMap->value); tmpMap->value=(char*) malloc((strlen(fpath)+1)*sizeof(char)); sprintf(tmpMap->value,"%s",fpath); pszDataSource=tmpMap->value; #ifdef DEBUGMS fprintf(stderr,"*** DEBUG %s\n",pszDataSource); #endif } } VSIFCloseL(vsif); i++; } } free(sdsName); free(dsName); OGRDataSourceH poDS = NULL; OGRSFDriverH *poDriver = NULL; poDS = OGROpen( pszDataSource, FALSE, poDriver ); if( poDS == NULL ){ #ifdef DEBUGMS fprintf(stderr,"Unable to access the DataSource %s\n",pszDataSource); #endif setMapInMaps(conf,"lenv","message","Unable to open datasource in read only mode"); OGR_DS_Destroy(poDS); OGRCleanupAll(); #ifdef DEBUGMS fprintf(stderr,"Unable to access the DataSource, exit! \n"); #endif return -1; } int iLayer = 0; for( iLayer=0; iLayer < OGR_DS_GetLayerCount(poDS); iLayer++ ){ OGRLayerH poLayer = OGR_DS_GetLayer(poDS,iLayer); if( poLayer == NULL ){ #ifdef DEBUGMS fprintf(stderr,"Unable to access the DataSource Layer \n"); #endif setMapInMaps(conf,"lenv","message","Unable to open datasource in read only mode"); return -1; } /** * Add a new layer set name, data */ if(msGrowMapLayers(m)==NULL){ return -1; } if(initLayer((m->layers[m->numlayers]), m) == -1){ return -1; } layerObj* myLayer=m->layers[m->numlayers]; #ifdef DEBUGMS dumpMaps(output); #endif myLayer->name = zStrdup(output->name); myLayer->tileitem=NULL; myLayer->data = zStrdup(OGR_L_GetName(poLayer)); myLayer->connection = zStrdup(pszDataSource); myLayer->index = m->numlayers; myLayer->dump = MS_TRUE; myLayer->status = MS_ON; msConnectLayer(myLayer,MS_OGR,pszDataSource); /** * Detect the Geometry Type or use Polygon */ if(OGR_L_GetGeomType(poLayer) != wkbUnknown){ switch(OGR_L_GetGeomType(poLayer)){ case wkbPoint: case wkbMultiPoint: case wkbPoint25D: case wkbMultiPoint25D: #ifdef DEBUGMS fprintf(stderr,"POINT DataSource Layer \n"); #endif myLayer->type = MS_LAYER_POINT; break; case wkbLineString : case wkbMultiLineString : case wkbLineString25D: case wkbMultiLineString25D: #ifdef DEBUGMS fprintf(stderr,"LINE DataSource Layer \n"); #endif myLayer->type = MS_LAYER_LINE; break; case wkbPolygon: case wkbMultiPolygon: case wkbPolygon25D: case wkbMultiPolygon25D: #ifdef DEBUGMS fprintf(stderr,"POLYGON DataSource Layer \n"); #endif myLayer->type = MS_LAYER_POLYGON; break; default: myLayer->type = MS_LAYER_POLYGON; break; } }else myLayer->type = MS_LAYER_POLYGON; /** * Detect spatial reference or use WGS84 */ OGRSpatialReferenceH srs=OGR_L_GetSpatialRef(poLayer); if(srs!=NULL){ char *wkt=NULL; OSRExportToWkt(srs,&wkt); setSrsInformations(output,m,myLayer,wkt); free(wkt); } else{ addToMap(output->content,"crs","EPSG:4326"); addToMap(output->content,"crs_isGeographic","true"); msLoadProjectionStringEPSG(&m->projection,"EPSG:4326"); msInsertHashTable(&(m->web.metadata), "ows_srs", "EPSG:4326 EPSG:900913"); msInsertHashTable(&(myLayer->metadata), "ows_srs", "EPSG:4326 EPSG:900913"); } OGREnvelope oExt; if (OGR_L_GetExtent(poLayer,&oExt, TRUE) == OGRERR_NONE){ setMsExtent(output,m,myLayer,oExt.MinX, oExt.MinY, oExt.MaxX, oExt.MaxY); } /** * Detect the FID column or use the first attribute field as FID */ char *fid=(char*)OGR_L_GetFIDColumn(poLayer); if(strlen(fid)==0){ OGRFeatureDefnH def=OGR_L_GetLayerDefn(poLayer); int fIndex=0; for(fIndex=0;fIndexmetadata), "gml_featureid", fid); msInsertHashTable(&(myLayer->metadata), "gml_include_items", "all"); msInsertHashTable(&(myLayer->metadata), "ows_name", output->name); map* tmpMap=getMap(output->content,"title"); if(tmpMap!=NULL) msInsertHashTable(&(myLayer->metadata), "ows_title", tmpMap->value); else msInsertHashTable(&(myLayer->metadata), "ows_title", "Default Title"); if(msGrowLayerClasses(myLayer) == NULL) return -1; if(initClass((myLayer->CLASS[myLayer->numclasses])) == -1) return -1; myLayer->CLASS[myLayer->numclasses]->type = myLayer->type; if(msGrowClassStyles(myLayer->CLASS[myLayer->numclasses]) == NULL) return -1; if(initStyle(myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]) == -1) return -1; /** * Apply msStyle else fallback to the default style */ tmpMap=getMap(output->content,"msStyle"); if(tmpMap!=NULL) msUpdateStyleFromString(myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles],tmpMap->value,0); else{ /** * Set style */ myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.red=125; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.green=125; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.blue=255; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->outlinecolor.red=80; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->outlinecolor.green=80; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->outlinecolor.blue=80; /** * Set specific style depending on type */ if(myLayer->type == MS_LAYER_POLYGON) myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->width=3; if(myLayer->type == MS_LAYER_LINE){ myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->width=3; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->outlinewidth=1.5; } if(myLayer->type == MS_LAYER_POINT){ myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->symbol=1; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->size=15; } } myLayer->CLASS[myLayer->numclasses]->numstyles++; myLayer->numclasses++; m->layerorder[m->numlayers] = m->numlayers; m->numlayers++; } OGR_DS_Destroy(poDS); OGRCleanupAll(); return 1; } /** * Try to open a raster output and define the corresponding layer in the MAPFILE * * @param conf the conf maps containing the main.cfg settings * @param output the specific output maps * @param m the mapObj */ int tryGdal(maps* conf,maps* output,mapObj* m){ map* tmpMap=getMap(output->content,"storage"); char *pszFilename=tmpMap->value; GDALDatasetH hDataset; GDALRasterBandH hBand; double adfGeoTransform[6]; int i, iBand; /** * Try to open the DataSource using GDAL */ GDALAllRegister(); hDataset = GDALOpen( pszFilename, GA_ReadOnly ); if( hDataset == NULL ){ #ifdef DEBUGMS fprintf(stderr,"Unable to access the DataSource %s \n",pszFilename); #endif setMapInMaps(conf,"lenv","message","gdalinfo failed - unable to open"); GDALDestroyDriverManager(); return -1; } #ifdef DEBUGMS fprintf(stderr,"Accessing the DataSource %s %d\n",pszFilename,__LINE__); #endif /** * Add a new layer set name, data */ if(msGrowMapLayers(m)==NULL){ return -1; } if(initLayer((m->layers[m->numlayers]), m) == -1){ return -1; } m->layers[m->numlayers]->index=m->numlayers; layerObj* myLayer=m->layers[m->numlayers]; myLayer->name = zStrdup(output->name); myLayer->tileitem=NULL; myLayer->data = zStrdup(pszFilename); myLayer->index = m->numlayers; myLayer->dump = MS_TRUE; myLayer->status = MS_ON; myLayer->type = MS_LAYER_RASTER; char *title=output->name; tmpMap=getMap(output->content,"title"); if(tmpMap!=NULL) title=tmpMap->value; char *abstract=output->name; tmpMap=getMap(output->content,"abstract"); if(tmpMap!=NULL) abstract=tmpMap->value; msInsertHashTable(&(myLayer->metadata), "ows_label", title); msInsertHashTable(&(myLayer->metadata), "ows_title", title); msInsertHashTable(&(myLayer->metadata), "ows_abstract", abstract); msInsertHashTable(&(myLayer->metadata), "ows_rangeset_name", output->name); msInsertHashTable(&(myLayer->metadata), "ows_rangeset_label", title); /** * Set Map Size to the raster size */ m->width=GDALGetRasterXSize( hDataset ); m->height=GDALGetRasterYSize( hDataset ); /** * Set projection using Authority Code and Name if available or fallback to * proj4 definition if available or fallback to default EPSG:4326 */ const char *tRef=GDALGetProjectionRef( hDataset ); if( tRef != NULL && strlen(tRef)>0 ){ char *pszProjection; pszProjection = (char *) GDALGetProjectionRef( hDataset ); #ifdef DEBUGMS fprintf(stderr,"Accessing the DataSource %s\n",GDALGetProjectionRef( hDataset )); #endif setSrsInformations(output,m,myLayer,pszProjection); }else{ fprintf(stderr,"NO SRS FOUND ! %s\n",GDALGetProjectionRef( hDataset )); } /** * Set extent */ if( GDALGetGeoTransform( hDataset, adfGeoTransform ) == CE_None ){ if( adfGeoTransform[2] == 0.0 && adfGeoTransform[4] == 0.0 ){ double minX = adfGeoTransform[0] + adfGeoTransform[2] * GDALGetRasterYSize(hDataset); double minY = adfGeoTransform[3] + adfGeoTransform[5] * GDALGetRasterYSize(hDataset); double maxX = adfGeoTransform[0] + adfGeoTransform[1] * GDALGetRasterXSize(hDataset); double maxY = adfGeoTransform[3] + adfGeoTransform[4] * GDALGetRasterXSize(hDataset); setMsExtent(output,m,myLayer,minX,minY,maxX,maxY); } } /** * Extract information about available bands to set the bandcount and the * processing directive */ char nBands[2]; int nBandsI=GDALGetRasterCount( hDataset ); sprintf(nBands,"%d",GDALGetRasterCount( hDataset )); msInsertHashTable(&(myLayer->metadata), "ows_bandcount", nBands); if(nBandsI>=3) msLayerAddProcessing(myLayer,"BANDS=1,2,3"); else if(nBandsI>=2) msLayerAddProcessing(myLayer,"BANDS=1,2"); else msLayerAddProcessing(myLayer,"BANDS=1"); /** * Name available Bands */ char lBands[6]; char *nameBands=NULL; for( iBand = 0; iBand < nBandsI; iBand++ ){ sprintf(lBands,"Band%d",iBand+1); if(nameBands==NULL){ nameBands=(char*)malloc((strlen(lBands)+1)*sizeof(char)); sprintf(nameBands,"%s",lBands); }else{ if(iBand<4){ char *tmpS=zStrdup(nameBands); nameBands=(char*)realloc(nameBands,(strlen(nameBands)+strlen(lBands)+1)*sizeof(char)); sprintf(nameBands,"%s %s",tmpS,lBands); free(tmpS); } } } msInsertHashTable(&(myLayer->metadata), "ows_bandnames", nameBands); /** * Loops over metadata informations to setup specific informations */ for( iBand = 0; iBand < nBandsI; iBand++ ){ //int bGotNodata;//, bSuccess; double adfCMinMax[2]/*, dfNoData*/; //int nBlockXSize, nBlockYSize, nMaskFlags; //double /*dfMean, dfStdDev*/; hBand = GDALGetRasterBand( hDataset, iBand+1 ); CPLErrorReset(); GDALComputeRasterMinMax( hBand, FALSE, adfCMinMax ); char tmpN[21]; sprintf(tmpN,"Band%d",iBand+1); if (CPLGetLastErrorType() == CE_None){ char tmpMm[100]; sprintf(tmpMm,"%.3f %.3f",adfCMinMax[0],adfCMinMax[1]); char tmpI[21]; sprintf(tmpI,"%s_interval",tmpN); msInsertHashTable(&(myLayer->metadata), tmpI, tmpMm); map* test=getMap(output->content,"msClassify"); if(test!=NULL && strncasecmp(test->value,"true",4)==0){ /** * Classify one band raster pixel value using regular interval */ int _tmpColors[10][3]={ {102,153,204}, {51,102,153}, {102,102,204}, {51,204,0}, {153,255,102}, {204,255,102}, {102,204,153}, {255,69,64}, {255,192,115}, {255,201,115} }; if(nBandsI==1){ double delta=adfCMinMax[1]-adfCMinMax[0]; double interval=delta/10; double cstep=adfCMinMax[0]; for(i=0;i<10;i++){ /** * Create a new class */ if(msGrowLayerClasses(myLayer) == NULL) return -1; if(initClass((myLayer->CLASS[myLayer->numclasses])) == -1) return -1; myLayer->CLASS[myLayer->numclasses]->type = myLayer->type; if(msGrowClassStyles(myLayer->CLASS[myLayer->numclasses]) == NULL) return -1; if(initStyle(myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]) == -1) return -1; /** * Set class name */ char className[7]; sprintf(className,"class%d",i); myLayer->CLASS[myLayer->numclasses]->name=zStrdup(className); /** * Set expression */ char expression[1024]; if(i+1<10) sprintf(expression,"([pixel]>=%.3f AND [pixel]<%.3f)",cstep,cstep+interval); else sprintf(expression,"([pixel]>=%.3f AND [pixel]<=%.3f)",cstep,cstep+interval); msLoadExpressionString(&myLayer->CLASS[myLayer->numclasses]->expression,expression); /** * Set color */ myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.red=_tmpColors[i][0]; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.green=_tmpColors[i][1]; myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.blue=_tmpColors[i][2]; cstep+=interval; myLayer->CLASS[myLayer->numclasses]->numstyles++; myLayer->numclasses++; } char tmpMm[100]; sprintf(tmpMm,"%.3f %.3f",adfCMinMax[0],adfCMinMax[1]); } } } if( strlen(GDALGetRasterUnitType(hBand)) > 0 ){ char tmpU[21]; sprintf(tmpU,"%s_band_uom",tmpN); msInsertHashTable(&(myLayer->metadata), tmpU, GDALGetRasterUnitType(hBand)); } } m->layerorder[m->numlayers] = m->numlayers; m->numlayers++; GDALClose( hDataset ); GDALDestroyDriverManager(); CPLCleanupTLS(); return 1; } /** * Create a MapFile for WMS, WFS or WCS Service output * * @param conf the conf maps containing the main.cfg settings * @param outputs a specific output maps */ void outputMapfile(maps* conf,maps* outputs){ /** * First store the value on disk */ map* mime=getMap(outputs->content,"mimeType"); char *ext="data"; if(mime!=NULL) if(strncasecmp(mime->value,"application/json",16)==0) ext="json"; map* tmpMap=getMapFromMaps(conf,"main","dataPath"); map* sidMap=getMapFromMaps(conf,"lenv","usid"); char *pszDataSource=(char*)malloc((strlen(tmpMap->value)+strlen(sidMap->value)+strlen(outputs->name)+17)*sizeof(char)); sprintf(pszDataSource,"%s/ZOO_DATA_%s_%s.%s",tmpMap->value,outputs->name,sidMap->value,ext); int f=zOpen(pszDataSource,O_WRONLY|O_CREAT,S_IRUSR|S_IWUSR|S_IRGRP|S_IWGRP|S_IROTH|S_IWOTH); map *gfile=getMap(outputs->content,"generated_file"); if(gfile!=NULL){ readGeneratedFile(conf,outputs->content,gfile->value); } map* sizeMap=getMap(outputs->content,"size"); map* vData=getMap(outputs->content,"value"); if(sizeMap!=NULL){ zWrite(f,vData->value,atoi(sizeMap->value)*sizeof(char)); } else{ zWrite(f,vData->value,(strlen(vData->value)+1)*sizeof(char)); } close(f); addToMap(outputs->content,"storage",pszDataSource); free(pszDataSource); /* * Create an empty map, set name, default size and extent */ mapObj *myMap=msNewMapObj(); free(myMap->name); myMap->name=zStrdup("ZOO-Project_WXS_Server"); msMapSetSize(myMap,2048,2048); msMapSetExtent(myMap,-1,-1,1,1); /* * Set imagepath and imageurl using tmpPath and tmpUrl from main.cfg */ map *tmp1=getMapFromMaps(conf,"main","tmpPath"); myMap->web.imagepath=zStrdup(tmp1->value); tmp1=getMapFromMaps(conf,"main","tmpUrl"); myMap->web.imageurl=zStrdup(tmp1->value); /* * Define supported output formats */ outputFormatObj *o1=msCreateDefaultOutputFormat(NULL,"AGG/PNG","png"); o1->imagemode=MS_IMAGEMODE_RGBA; o1->transparent=MS_TRUE; o1->inmapfile=MS_TRUE; msAppendOutputFormat(myMap,msCloneOutputFormat(o1)); msFreeOutputFormat(o1); #ifdef USE_KML outputFormatObj *o2=msCreateDefaultOutputFormat(NULL,"KML","kml"); o2->inmapfile=MS_TRUE; msAppendOutputFormat(myMap,msCloneOutputFormat(o2)); msFreeOutputFormat(o2); #endif outputFormatObj *o3=msCreateDefaultOutputFormat(NULL,"GDAL/GTiff","tiff"); if(!o3) fprintf(stderr,"Unable to initialize GDAL driver !\n"); else{ o3->imagemode=MS_IMAGEMODE_BYTE; o3->inmapfile=MS_TRUE; msAppendOutputFormat(myMap,msCloneOutputFormat(o3)); msFreeOutputFormat(o3); } outputFormatObj *o4=msCreateDefaultOutputFormat(NULL,"GDAL/AAIGRID","grd"); if(!o4) fprintf(stderr,"Unable to initialize GDAL driver !\n"); else{ o4->imagemode=MS_IMAGEMODE_INT16; o4->inmapfile=MS_TRUE; msAppendOutputFormat(myMap,msCloneOutputFormat(o4)); msFreeOutputFormat(o4); } #ifdef USE_CAIRO outputFormatObj *o5=msCreateDefaultOutputFormat(NULL,"CAIRO/PNG","cairopng"); if(!o5) fprintf(stderr,"Unable to initialize CAIRO driver !\n"); else{ o5->imagemode=MS_IMAGEMODE_RGBA; o5->transparent=MS_TRUE; o5->inmapfile=MS_TRUE; msAppendOutputFormat(myMap,msCloneOutputFormat(o5)); msFreeOutputFormat(o5); } #endif /* * Set default projection to EPSG:4326 */ msLoadProjectionStringEPSG(&myMap->projection,"EPSG:4326"); myMap->transparent=1; /** * Set metadata extracted from main.cfg file maps */ maps* cursor=conf; map* correspondance=getCorrespondance(); while(cursor!=NULL){ map* _cursor=cursor->content; map* vMap; while(_cursor!=NULL){ if((vMap=getMap(correspondance,_cursor->name))!=NULL){ if (msInsertHashTable(&(myMap->web.metadata), vMap->value, _cursor->value) == NULL){ #ifdef DEBUGMS fprintf(stderr,"Unable to add metadata"); #endif return; } } _cursor=_cursor->next; } cursor=cursor->next; } freeMap(&correspondance); free(correspondance); /* * Set mapserver PROJ_LIB/GDAL_DATA or any other config parameter from * the main.cfg [mapserver] section if any */ maps *tmp3=getMaps(conf,"mapserver"); if(tmp3!=NULL){ map* tmp4=tmp3->content; while(tmp4!=NULL){ msSetConfigOption(myMap,tmp4->name,tmp4->value); tmp4=tmp4->next; } } /** * Set a ows_rootlayer_title, */ if (msInsertHashTable(&(myMap->web.metadata), "ows_rootlayer_name", "ZOO_Project_Layer") == NULL){ #ifdef DEBUGMS fprintf(stderr,"Unable to add metadata"); #endif return; } if (msInsertHashTable(&(myMap->web.metadata), "ows_rootlayer_title", "ZOO_Project_Layer") == NULL){ #ifdef DEBUGMS fprintf(stderr,"Unable to add metadata"); #endif return; } /** * Enable all the WXS requests using ows_enable_request * see http://mapserver.org/trunk/development/rfc/ms-rfc-67.html */ if (msInsertHashTable(&(myMap->web.metadata), "ows_enable_request", "*") == NULL){ #ifdef DEBUGMS fprintf(stderr,"Unable to add metadata"); #endif return; } msInsertHashTable(&(myMap->web.metadata), "ows_srs", "EPSG:4326"); if(tryOgr(conf,outputs,myMap)<0) if(tryGdal(conf,outputs,myMap)<0) return ; tmp1=getMapFromMaps(conf,"main","dataPath"); char *tmpPath=(char*)malloc((13+strlen(tmp1->value))*sizeof(char)); sprintf(tmpPath,"%s/symbols.sym",tmp1->value); msInitSymbolSet(&myMap->symbolset); myMap->symbolset.filename=zStrdup(tmpPath); free(tmpPath); map* sid=getMapFromMaps(conf,"lenv","usid"); char *mapPath= (char*)malloc((7+strlen(sid->value)+strlen(outputs->name)+strlen(tmp1->value))*sizeof(char)); sprintf(mapPath,"%s/%s_%s.map",tmp1->value,outputs->name,sid->value); msSaveMap(myMap,mapPath); free(mapPath); msGDALCleanup(); msFreeMap(myMap); }