title: Geoserver - How can I make segmented roads appear as one line?

text: I have shp file with segmented roads. I've edited the SLD in geoserver and the final result looks like this : ![enter image description here](http://i.stack.imgur.com/cUkng.png) How can I make the roads look continuous ? This is the roads SLD I use :                       <?xml version="1.0" encoding="UTF-8"?>     <StyledLayerDescriptor xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" version="1.1.0" xmlns:xlink="http://www.w3.org/1999/xlink" xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.1.0/StyledLayerDescriptor.xsd" xmlns:se="http://www.opengis.net/se">       <NamedLayer>         <se:Name>roads</se:Name>         <UserStyle>           <se:Name>roads</se:Name>           <se:FeatureTypeStyle>             <se:Rule>               <se:Name>first</se:Name>               <se:Description>                 <se:Title>first</se:Title>               </se:Description>               <ogc:Filter>                 <ogc:Or>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>footway</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>footpath</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>steps</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>cycleway</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>pedestrian</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>track</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>residential</ogc:Literal>                   </ogc:PropertyIsEqualTo>                 </ogc:Or>               </ogc:Filter>               <se:MinScaleDenominator>9000</se:MinScaleDenominator>               <se:MaxScaleDenominator>30000</se:MaxScaleDenominator>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#646450</se:SvgParameter>                   <se:SvgParameter name="stroke-width">3</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                   <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">11.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#ffffff</se:SvgParameter>                   <se:SvgParameter name="stroke-width">3</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">11.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>                  </se:Rule>           </se:FeatureTypeStyle>                     <se:FeatureTypeStyle>             <se:Rule>               <se:Name>second</se:Name>               <se:Description>                 <se:Title>second</se:Title>               </se:Description>               <ogc:Filter>                 <ogc:Or>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>footway</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>footpath</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>steps</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>cycleway</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>pedestrian</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>track</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>bridleway</ogc:Literal>                   </ogc:PropertyIsEqualTo>                     <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>residential</ogc:Literal>                   </ogc:PropertyIsEqualTo>                 </ogc:Or>               </ogc:Filter>               <se:MinScaleDenominator>30000</se:MinScaleDenominator>               <se:MaxScaleDenominator>150000</se:MaxScaleDenominator>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#c7c7c7</se:SvgParameter>                   <se:SvgParameter name="stroke-width">4</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>             </se:Rule>           </se:FeatureTypeStyle>              <se:FeatureTypeStyle>             <se:Rule>               <se:Name>third</se:Name>               <se:Description>                 <se:Title>third</se:Title>               </se:Description>               <ogc:Filter>                 <ogc:Or>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>footway</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>footpath</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>steps</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>cycleway</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>pedestrian</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>track</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>bridleway</ogc:Literal>                   </ogc:PropertyIsEqualTo>                     <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>residential</ogc:Literal>                   </ogc:PropertyIsEqualTo>                 </ogc:Or>               </ogc:Filter>               <se:MinScaleDenominator>2</se:MinScaleDenominator>               <se:MaxScaleDenominator>9000</se:MaxScaleDenominator>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#505050</se:SvgParameter>                   <se:SvgParameter name="stroke-width">5</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#f8f8f8</se:SvgParameter>                   <se:SvgParameter name="stroke-width">5</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>             </se:Rule>           </se:FeatureTypeStyle>              <se:FeatureTypeStyle>             <se:Rule>               <se:Name>forth</se:Name>               <se:Description>                 <se:Title>forth</se:Title>               </se:Description>               <ogc:Filter>                 <ogc:Or>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>secondary</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>tertiary</ogc:Literal>                   </ogc:PropertyIsEqualTo>                 </ogc:Or>               </ogc:Filter>               <se:MinScaleDenominator>1</se:MinScaleDenominator>               <se:MaxScaleDenominator>20000</se:MaxScaleDenominator>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#505050</se:SvgParameter>                   <se:SvgParameter name="stroke-width">5</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#f8f8f8</se:SvgParameter>                   <se:SvgParameter name="stroke-width">5</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>             </se:Rule>           </se:FeatureTypeStyle>              <se:FeatureTypeStyle>             <se:Rule>               <se:Name>fifth</se:Name>               <se:Description>                 <se:Title>fifth</se:Title>               </se:Description>               <ogc:Filter>                 <ogc:Or>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>secondary</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>tertiary</ogc:Literal>                   </ogc:PropertyIsEqualTo>                 </ogc:Or>               </ogc:Filter>               <se:MinScaleDenominator>20000</se:MinScaleDenominator>               <se:MaxScaleDenominator>200000</se:MaxScaleDenominator>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#505050</se:SvgParameter>                   <se:SvgParameter name="stroke-width">1</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#ffffff</se:SvgParameter>                   <se:SvgParameter name="stroke-width">0.5</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>             </se:Rule>           </se:FeatureTypeStyle>              <se:FeatureTypeStyle>             <se:Rule>               <se:Name>sixth</se:Name>               <se:Description>                 <se:Title>sixth</se:Title>               </se:Description>               <ogc:Filter>                 <ogc:Or>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>primary</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>primary_link</ogc:Literal>                   </ogc:PropertyIsEqualTo>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>fclass</ogc:PropertyName>                     <ogc:Literal>trunk</ogc:Literal>                   </ogc:PropertyIsEqualTo>                 </ogc:Or>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#505050</se:SvgParameter>                   <se:SvgParameter name="stroke-width">4</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#ffe168</se:SvgParameter>                   <se:SvgParameter name="stroke-width">4</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>             </se:Rule>           </se:FeatureTypeStyle>              <se:FeatureTypeStyle>             <se:Rule>               <se:Name>eighth</se:Name>               <se:Description>                 <se:Title>eighth</se:Title>               </se:Description>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>        <ogc:PropertyName>fclass</ogc:PropertyName>                 <ogc:Literal>motorway</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#4c2600</se:SvgParameter>                   <se:SvgParameter name="stroke-width">1</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#fa9e25</se:SvgParameter>                   <se:SvgParameter name="stroke-width">0.904762</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>               <se:TextSymbolizer>              <se:Label>                <ogc:PropertyName>name</ogc:PropertyName>              </se:Label>                  <se:Font>                                 <se:SvgParameter name="font-family">David</se:SvgParameter>                                 <se:SvgParameter name="font-size">9.0</se:SvgParameter>                                 <se:SvgParameter name="font-style">normal</se:SvgParameter>                                 <se:SvgParameter name="font-weight">normal</se:SvgParameter>                             </se:Font>              <se:LabelPlacement>                <se:LinePlacement />              </se:LabelPlacement>                   <se:VendorOption name="followLine">true</se:VendorOption>            </se:TextSymbolizer>             </se:Rule>           </se:FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>


title: Convert shapefile to raster using arcmap 10

text: Hi I have a question regarding how to convert a shapefile to a raster image using arcmap 10. If anyone could help on this matter I would greatly appreciate it and also the functions within the conversion toolbox.


title: Import csv file and use "toggle editing" and "field calculator" without Saving it as a Shapefile.

text: I am importing a .csv file in QGIS and I have also made the *.csvt file where I have specified my data types. I still cannot use "Toggle editing" and "field calculator".Whether I must save it as a shapefile and import it again? Is there any way to work with the .csv file without saving it as a shapefile and import it again?


title: How to change the color of a selected feature with pyQGIS?

text: is it possible to change the color of a selected feature in qgis 1.8? To highlight a feature I use:               vlayer.setSelectedFeatures(selectedList)     box = vlayer.boundingBoxOfSelected()     qgis.utils.iface.mapCanvas().setExtent(box)     qgis.utils.iface.mapCanvas().refresh()      But the selected featue is colored yellow. Is it possible to define for example a red color for this feature?


title: UTFgrid support in QGIS

text: Starting with GDAL 1.10, the MBTiles driver allows to read rasters in the MBTiles format. The driver can also retrieve pixel attributes encoded according to the UTFGrid specification available in some MBTiles file. Thus, installing QGIS (2.0) with gdal 1.10 support will enable QGIS to open MBTiles. This does work (on Ubuntu 13.04 with gdal 1.10.1 and QGIS 2.0 installed). When using the identifier, one get's the rgba band values. However, is there any way to enable UTFgrid support so the user can get the embedded raster 'attribute' values?


title: Projection system differs

text: There is a problem . I want to define projection for a Shape file. The projection system is: "NAD_1983_StatePlane_Montana_FIPS_2500" which is refered in Qgis 1.9 as EPSG:102700 .The same projection system is not available in Qgis 1.8 .Here it is designated as "NAD 83/Mountana(ft)" and has EPSG:2256 and "NAD83/Mountana" and has EPSG:32100.My problem is if I want to define projection in Qgis 1.8 what projection system should I use? EPSG:2256 or EPSG:32100?


title: How to cluster Raster data

text: I have thousands of Raster data concerning the ice coverage at the pole at different times (days) over several years. Now I would like to perform a cluster analysis. I want to split the data into different periods (e.g. 4 seasons) with respect to the ice cover. Can I perform such a cluster analysis in R? And if, how? (I am a beginner in R, so please keep it as simple as possible. Thanks.) Or is it possible to use Grass GIS or QGIS (considering the large amount of data)?


title: How to extract german restaurants from planet.osm?

text: Want to get the data out of http://wiki.openstreetmap.org/wiki/Planet.osm **UPDATE:** i just come back from the site http://wiki.openstreetmap.org/wiki/Planet.osm which has got some files that belong to the openstreetmap-Project. i downloaded a similar file that is derived from a German Mirrorsite - this one here: http://ftp5.gwdg.de/pub/misc/openstreetmap/download.geofabrik.de/ i downloaded the file called germany.osm,bz2 according to this instruction > type bzip2 -d planet.osm.bz2; or your OS may support double-click unpacking. > See Wikipedia's list of compression programs. i did the following - but without success-,               martin@linux-wyee:~/poi_abzug/geofabrik> bzip2 -d germany.osm,bz2;     bzip2: Can't open input file germany.osm,bz2: No such file or directory.      why - can anybody explain why i did not succeed? **end of update;** further **questions:** Is this a XML-File ? BTW; i only need the data for the Pois in Germany: can i use this file here - http://ftp5.gwdg.de/pub/misc/openstreetmap/ As this might be of general interest i think that many users can benefit from an answer to this question. So my **question** is: How to tread the above mentioned dataset - how to read it - on a OpenSuse-Linux-machine!? Should i do it in the terminal? Note - i want to store all in a csv file - and i look for all the restaurants... What is the easiest way to do this!? As i am a newbie - thx for your patience!


title: How to convert Garmin Waypoints to Google Mercator?

text: I have georef some topocards with QGIS and the Google Hybrid Layer. So, this cards are connected with the Google Mercator System. Now I have some waypoints taken with Garmin GPSmap62 (in WGS84). When I try to load both datasets in QGIS, I get a mismatch of around 2000 km (don't know, the topocards are in bulgaria and the waypoints somewhere in the southern atlantic). I tried to convert the topos to WGS84 and/or the waypoints to Google Mercator. Either qgis crashed down or nothing happens. I'm working with Win7, and QGIS 1.7.4 \edit: First the Topocard. And I try to load up the gpx File, but I'm not sure how _gg_. Ah okay, it isn't allowed me to upload pictures, cause I'm a new user. sry


title: How to connect the selectionChanged signal with pyqt?

text: Hi I am developing a plugin for Qgis in python, I need to connect the sselectionChanged signal emitted when a feature of the layer is selected, I could not find any examples on internet, here is what I have done as far                   QObject.connect(self.iface.mapCanvas(),SIGNAL("selectionChanged(QgsMapLayer)"), self.test)       on the log of qgis I get this warning:                   Warning: Object::connect: (receiver name: 'MainWindow')      can you help me? please, thanks in advance


title: Save multiple selection sets in QGIS

text: I'm looking for a feature in QGIS but I am not sure if it exists.   I don't know the name of such feature, but I call it multiple selection sets. It exists in TransCAD and I find it very useful. It basically allows one to save the selected features under a name, and have multiple sets that one can compose (union, intersection, difference, etc).   I am under the impression that QGIS allows only to save the selected set as a new file, which is not what I am looking for.   Does it exists in QGIS? How can one use it? If not, is there a work-around?


title: How to make two feature layers congruent using ArcGIS for Desktop?

text: I have two feature layers (polyline/polygon) which should be congruent but there's an offset of about 220 meters. The Georeferencing tool only works with raster layers. How do I do the equivalent with a feature layer?


title: Is it possible to georeference a shapefile?

text: I need to align a polyline shapefile with a pre-existing one. This new one has a couple of edits, but for some reason shows up in a different location (I believe the person who made this edit did so in a different program other than ArcMap). I've already had someone take a look at the projection to see if it was an issue and so far it seems to be okay. How I can change over to align the two lines?


title: GeoJson layers with leaflet in IE8

text: I have built a website displaying a lot of line shapes, using geojson, on a leaflet map. This is working relatively well in most browsers, though speed could be improved. As usual IE8 is being the thorn in my side, originally it was not known that IE8 would be a requirement of this project. I currently cannot get the geojson lines to display on my large searchable map. The lines should display after a couple zooms. Leaflet is not throwing an error, the lines just aren't visible. Has anyone else experienced such an issue, or have an alternative to geojson layers for IE8? My current setup is a MySQL database with spatial data from a shape file, which is then converted into geojson, for leaflet to render.


title: City density in GIS

text: There's a cite calles http://map.vmr.gov.ua/ in which each house has a type (for living, offices etc) and number of floors. I want to find out city density - i.e. to assign certain amount of people to each house, avaliable for living (using average family statistics and typical amount of families for each floor) to find out general density, make a sort of heat map and to find out most and least dense parts of city. I'm mostly new to GIS systems so I would be very greatfull for any suggestions how to solve that problem. Thank you in advance.


title: Arcmap classification data exclusion

text: I am presenting some percentage change data, and I created a color ramp such that negative values are shades of red and the positive values are shades are blue. I want to use a custom classification with data exclusion so that I create a classification scale for negative values and scale for positive values. However, data exclusion doesn't seem to be working. Is there a different way to achieve my goal? I am using Arcmap 10.1.


title: Google Earth tour LookAt not looking at same place as Placemark with same coordinates

text: I have tried to create a tour from a csv file. I want the tour to look at each Placemark in turn but the view is of a different set of locations when run in Google Earth. The Placemarks are plotting correctly (as far as I can tell) and the LookAt seems to be off. Please could someone tell me what I'm doing wrong?               <Document>        <name>Test</name>        <open>1</open>        <gx:Tour>        <name>Play me!</name>        <gx:Playlist>        <gx:FlyTo>        <gx:duration>10</gx:duration>        <LookAt>            <longtitude>-0.0441561</longtitude>            <latitude>51.4354051</latitude>            <altitude>400</altitude>            <heading>1</heading>            <tilt>0</tilt>            <range>1</range>            <altitudeMode>relativeToGround</altitudeMode>        </LookAt>        </gx:FlyTo>                <gx:wait>                    <gx:duration>10</gx:duration>                </gx:wait>        <gx:FlyTo>        <gx:duration>10</gx:duration>        <LookAt>            <longtitude>-0.0382701</longtitude>            <latitude>51.452901</latitude>            <altitude>400</altitude>            <heading>1</heading>            <tilt>0</tilt>            <range>1</range>            <altitudeMode>relativeToGround</altitudeMode>        </LookAt>        </gx:FlyTo>                <gx:wait>                    <gx:duration>10</gx:duration>                </gx:wait>        <gx:FlyTo>        <gx:duration>10</gx:duration>        <LookAt>            <longtitude>-0.0242591</longtitude>            <latitude>51.4303771</latitude>            <altitude>400</altitude>            <heading>1</heading>            <tilt>0</tilt>            <range>1</range>            <altitudeMode>relativeToGround</altitudeMode>        </LookAt>        </gx:FlyTo>                <gx:wait>                    <gx:duration>10</gx:duration>                </gx:wait>              <gx:TourControl>         </gx:TourControl>        </gx:Playlist>        </gx:Tour>      <Placemark>            <name>T1</name>            <description></description>            <Point>                <coordinates>-0.0441561,51.4354051</coordinates>            </Point>        </Placemark>        <Placemark>            <name>T2</name>            <description></description>            <Point>                <coordinates>-0.0382701,51.452901</coordinates>            </Point>        </Placemark>        <Placemark>            <name>T3</name>            <description></description>            <Point>                <coordinates>-0.0242591,51.4303771</coordinates>            </Point>        </Placemark>     </Document>


title: How I can insert feature attribute after feature drawing complete

text: I am new to OpenLayers. I am trying to insert wfs feature information to my map after drawing like this.   I can draw a WFS feature like (polygon, line, point) on map and after drawing the feature complete I can insert attribute information of feature (like polygon name, area, etc) into permanently. That is after refresh the map the feature (polygon) visible on the map and on click on the feature (polygon) the feature attribute information is display by popup.   I have done drawing feature on the map by finding the some demo on Google.   But I do not understand how to store the feature attribute. Currently I am using Geoserver and shapefile on Geoserver for layer and fetching this shapefile as layer in open layer by WFS service of Geoserver. Now my question is: How to store feature attribute by using Shapefile or Database (PostgresSQL)? If we use shapefile as WFS layer then, it is possible to store feature attribute after drawing the feature in shapefile by fetching it from Geoserver in OpenLayers. or For storing the feature attribute we should use Database (PostgresSQL)? I am totally new to OpenLayers. I don’t know which way is better. I want that I am able to draw feature on OpenLayers map and I can insert attribute information of feature by popup and save this feature on the map. My website Url is: http://mindcrewgis.com/


title: Using OS Mastermap SLDs in QGIS

text: In an answer to a question asked by @Canningmister, @tjmgis posted an excellent answer Styling OS Vectormap Local in QGIS? about making OS SLD files work in QGIS.   As an amateur user who uses QGIS to produce maps for my Parish Council I have struggled with styles for OS data in QGIS.   I can see how to edit the sld files using a text editor such as TextPad which I assume will work with Mastermap area and line layers but I am confused with making the text layer SLD work and have several questions:   1. Are the layers imported into QGIS from ESRI shapefiles or directly from gml?   2. How do you post process the text layer data to add extra columns without having the data in a database and if so which one?  I am familiar with executing SQL scripts in a database. Am I missing something - I can't get my head round this.


title: addins vs python tools, which makes more sense?

text: With the upcoming release of ArcGIS Pro and the eventual retirement of Desktop and ArcObjects does it make more sense to start new customizations in python or continue to use ArcObjects/addins for the time being? I have a new project about to start that will involve a custom toolbar for desktop with about 10 custom editing tools. Previously I have always carried out such a project in .NET and ArcObjects. But now I wonder:   1. Is the functionality now there to do this in python (At arcgis 10.1)?   2. If I continue to use ArcObjects what is a resonable expectation on the toolbar life expectancy? The current editing environment is ArcGIS desktop 10.1 and is likely to stay that way for the next 1 - 2 years. Any advice or opinions greatly appreciated


title: Associate points with polygon QGIS 2.2

text: I want to associate a point layer with a polygon layer so that I can do zonal statistics. In the image below, the polygon layer represents roughly the metro area of a certain city. The point represents the city. The polygon is a result of a polygonized process from a raster image. I want to use these new polygons as buffers to analyze the original raster. However, I need a way to associate the city points with the surrounding polygons. In other words, I want to find zonal statistics for each city ("oecdminusa") using the corresponding metro polyon ("polygontest"). When I ran zonal stats the first time, it simply named each of the polygons in "polygontest" as 1,2,3,4 etc. How do I make it correspond with the appropriate city? ![enter image description here](http://i.stack.imgur.com/SSsMj.jpg)


title: Filling Field2 to icorresponding values in Field1 of shapefile using Field Calculator of ArcGIS for Desktop?

text: I have to give explanation in a new field of the existing field which contains abbreviation. Example: Field "Geology" has values Q, Qs,N....(abbreviation), I have to give explanation of those in a new field "Geologic Age", where Q is Quaternary Sediments and Qs is Quaternary sand & dunes.. and it need to be done for more than 100 entries, so manual is practically time taking, something like "IF[GLG]=Q;THEN[Geologic_Age]=Quaternary sediments"


title: Looking for QGIS 2.0.1

text: So I deleted my old version of QGIS Dufour (2.0.1) with the intention of replacing it with the new version (2.2). It wasn't until I deleted my old version that I realized my Mac OS 10.6.8 is not compatible with the new QGIS. Where can I go to re-download my old version of QGIS (2.0.1)?


title: How to do geographic operations (intersect) on OSM XML?

text: I would like to create bicycle router which routes more through shade. For that I would like to use OSM. Most routers use OSM xml so my question is: How to add tags to a street that goes through a forest and write that back to OSM xml? Getting part of a street through forest is easy with POSTGIS. Some intersect operation probably. I can insert OSM to Postgresql with osm2psql but I have found no way to get it back in usable xml again. My current idea is parsing OSM XML and saving forest polygons somewhere. Then in another pass intersect each road with a forest polygon if bounding boxes intersects. If road is fully in forest tag it as in_forest. If not split road and tag forest part as in forest and other part as is. I would use python libraries (Shapely) or C++ if Python isn't fast enough. But this would probably take ages. My country is only 350MB unzipped OSM XML. Is there any better way? I looked into osmosis tag-transform but this is only for regex changes.


title: How to flip line based on a distance with another feature class

text: I have two feature classes Line and point. So what I need is to change the flow direction based on a distance with a point feature class. If the start vertices is far from the point, the system should flip the line. So I want the starting for all line features to be near the point feature class. what shall i do? Shall I do a model builder or python script?


title: How to transfer attribute of nearest Point to Line feature?

text: **Feature class 1: Addresses (Point) Feature class 2: Road centerlines (Polyline)** _Which is the most appropriate tool for transferring the Address attribute (e.g. Road No. or Bldg) to the nearest Road line segment?_ Tools reviewed and rejected **Attribute transfer ( Spatial adjustment) Near/Near table**


title: Need assistance setting center coordinate with GDALWARP

text: I'm having some issues converting my project into an MBTiles format. I'm trying to use gdalwarp to convert my main file to geotiff format, but I can't seem to get the center coordinate to how it is with my google maps custom map. My center coordinate is latitude:-53.722541 and longitude:-100.625000 using Google Maps. I'm not sure how to use gdalwarp to ensure my center coordinate to remain at the same position. My files dimensions are 8192x8192. The command I've been running is:               gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3785 -r bilinear -te -4096 -4096 4096 4096 input.png output.tif      The SRS projection I'm using with TileMill is:               +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over      I'm at a loss at which of those two I would need to change in order to keep my markers properly placed, and not need to redo each and every one of them. Thanks in advance for any support with this,


title: How to display a raster using Javascript API

text: It's all a proof of concept for myself and my learning but I have definitely hit a brick wall here! Here's what I've done and can get to work. I have written a Python script that:   1. Sets its workspace to and SDE instance.    2. Checks out a Spatial licencse.    3. Queries raster layers.    4. Takes a raster layer based on some string match - and then adds 100 to each pixel to create an output raster. Raster Algebra stuff....    5. The output raster is saved in the scratch workspace as "wibble". Hmmm, nice. **Next....** I can Share As A Geoprocessing Service and it does all the same stuff and saves an output raster to the ArcGIS Servers output directory with a Scratch FGDB there and a raster called "wibble" - hooray!! **Next...** I want to run my Geoprocessing task from the browser using JS API and that works too, and it puts the output raster in a jobs folder - looking something like this:               d:\arcgisserver\directories\arcgisjobs\testraster_gpserver\jccffbe1165404cf4b930b971c86be46a\scratch\scratch.gdb\wibble      My code for the geoprocessor looks like:               var gp = new Geoprocessor("http://ls8-webd-w8v:6080/arcgis/rest/services/TestRaster/GPServer/TestRaster");     var params = { "inputMessage": dom.byId("testMessage").value };             gp.submitJob(params, completeCallback, statusCallback);      So how do I now hook up my JavaScript code to get hold of this "wibble" raster and display it on a map in the browser. Now, I know you are not able to directly draw out rasters to a map service, it has to be an image, so you use getResultImageLayer from the geoprocessing object, but it still doesn't display, I think the code is wrong:               function completeCallback(jobInfo) {         imageParams = new ImageParameters();         imageParams.imageSpatialReference = map.spatialReference;         gp.getResultImageLayer(jobInfo.jobId, "wibble", imageParams, function (gpLayer) {             gpLayer.setOpacity(0.5);             map.addLayer(gpLayer);         });     }      I'm not even doing anything with ImageParameters - as I'm not sure what to do. Then I tried a different way, still nothing.               function completeCallback(jobInfo) {         mapserviceurl = "http://ls8-webd-w8v:6080/arcgis/rest/services/TestRaster/GPServer/TestRaster/jobs";         var mapurl = mapserviceurl + "/" + jobinfo.jobId;              var outputLayer = new ArcGISDynamicMapServiceLayer(mapurl, {             "id": "TestRaster",             "opacity": 0.7         });              map.addLayer(outputLayer);     };      Oh, by the way - here is the python that makes the geoprocessing script - it's nothing major!               import arcpy, os, sys     from arcpy.sa import *          inputSde = r"Database Connections\ls8-dbd-w8v-2_geodb_GEOLOGY.sde"     arcpy.env.overwriteOutput = True          inputString = arcpy.GetParameterAsText(0)          def main():         try:             if arcpy.CheckExtension("Spatial") == "Available":                          arcpy.CheckOutExtension("Spatial")             else:                 raise LicenseError                  arcpy.env.workspace = inputSde             outRas = None                  rastersStrings = arcpy.ListRasters("*", "All");                  for strRaster in rastersStrings:                 if strRaster == "GEOLOGY.MAAS_BATHYMETRY":                     path = os.path.join(inputSde, strRaster)                     raster = Raster(path)                     outRas = raster + 100                     if outRas is not None:                         tempFile = os.path.join(arcpy.env.scratchGDB, "wibble")                         outRas.save(tempFile)                     else:                         arcpy.AddWarning("Failed to make raster")              except LicenseError:             print "Spatial license is unavailable"           except:              arcpy.AddWarning("Error: {0}".format(sys.exc_info()[0]))         finally:             arcpy.CheckInExtension("Spatial")          if __name__ == '__main__':         main()      Any help - much appreciated guys! ## UPDATE I get the following Javascript error message for using the ImageParameters code (as above) - saying its "400 a bad request" - it looks like its looking in a folder "results" but I don't have that on the server - could it be to do with how I am saving the raster in the Python script?               http://ls8-webd-w8v:6080/arcgis/rest/services/TestRaster/GPServer/TestRaster/jobs/j7613d1e015a74976986c6f7447ea5ba5/results/wibble?dpi=96&transparent=true&format=png8&imageSR=4326&f=image&bbox=%7B%22xmin%22%3A-220.57570097228788%2C%22ymin%22%3A-48.21623734392513%2C%22xmax%22%3A220.57570130793587%2C%22ymax%22%3A135.59684693950146%2C%22spatialReference%22%3A%7B%22wkid%22%3A4326%2C%22latestWkid%22%3A4326%7D%7D&bboxSR=4326&size=960%2C400


title: CreateRasterCatalog on the fly using ArcObject 10 .NET SDK

text: My debugger crashes without showing any error. It crashes not even catches exception on the line:               rasterWorkspace.CreateRasterCatalog(catalogName, fields, "Shape", "Raster", "defaults");      Code snippet below:               CreateRasterCatalog(@"D:\TestfGDB\TestfGDB.gdb", "Test2");              public void CreateRasterCatalog(string fGDBPath, string catalogName)         {             try             {                 IRasterWorkspaceEx rasterWorkspace = OpenFileGDBWorkspace(fGDBPath);                 IFields fields = CreateField(null, null);                 IRasterCatalog rasterCatalog = rasterWorkspace.CreateRasterCatalog(catalogName, fields, "Shape", "Raster", "defaults");             }             catch (Exception e)             {                 string msg = e.Message;             }         }         /// <summary>         /// Opens a File geodatabase workspace.          /// </summary>         /// <param name="filePath">FGDB file path</param>         protected IRasterWorkspaceEx OpenFileGDBWorkspace(string filePath)         {             IWorkspaceFactory workspaceFactory = new FileGDBWorkspaceFactoryClass();             IRasterWorkspaceEx workspaceEx = (IRasterWorkspaceEx)workspaceFactory.OpenFromFile(filePath, 0);             return workspaceEx;         }              public IFields CreateField(ISpatialReference rasterSpatRef, ISpatialReference geoSpatRef)         {             IFields fields = new Fields();             IFieldsEdit fieldsEdit = (ESRI.ArcGIS.Geodatabase.IFieldsEdit)fields;             ESRI.ArcGIS.Geodatabase.IField field = new ESRI.ArcGIS.Geodatabase.FieldClass();             // create a user defined text field             ESRI.ArcGIS.Geodatabase.IFieldEdit fieldEdit = (ESRI.ArcGIS.Geodatabase.IFieldEdit)field;                  // setup field properties             fieldEdit.Name_2 = "OBJECTID";             fieldEdit.Type_2 = ESRI.ArcGIS.Geodatabase.esriFieldType.esriFieldTypeOID;             //fieldEdit.IsNullable_2 = true;             //fieldEdit.AliasName_2 = "Sample Field Column";             //fieldEdit.DefaultValue_2 = "test";             //fieldEdit.Editable_2 = true;             //fieldEdit.Length_2 = 100;             fieldsEdit.AddField(field);                       // setup field properties             fieldEdit.Name_2 = "Shape";             fieldEdit.Type_2 = ESRI.ArcGIS.Geodatabase.esriFieldType.esriFieldTypeGeometry;             //fieldEdit.IsNullable_2 = true;             //fieldEdit.Editable_2 = true;             fieldsEdit.AddField(field);                  // setup field properties             fieldEdit.Name_2 = "Raster";             fieldEdit.Type_2 = ESRI.ArcGIS.Geodatabase.esriFieldType.esriFieldTypeRaster;             //fieldEdit.IsNullable_2 = true;             //fieldEdit.Editable_2 = true;             fieldsEdit.AddField(field);                  // setup field properties             fieldEdit.Name_2 = "Name";             fieldEdit.Type_2 = ESRI.ArcGIS.Geodatabase.esriFieldType.esriFieldTypeString;             //fieldEdit.IsNullable_2 = true;             //fieldEdit.Editable_2 = true;             fieldsEdit.AddField(field);                  // setup field properties             fieldEdit.Name_2 = "Metadata";             fieldEdit.Type_2 = ESRI.ArcGIS.Geodatabase.esriFieldType.esriFieldTypeBlob;             //fieldEdit.IsNullable_2 = true;             //fieldEdit.Editable_2 = true;             fieldsEdit.AddField(field);                  // setup field properties             fieldEdit.Name_2 = "Shape_Length";             fieldEdit.Type_2 = ESRI.ArcGIS.Geodatabase.esriFieldType.esriFieldTypeDouble;             //fieldEdit.IsNullable_2 = true;             //fieldEdit.Editable_2 = true;             fieldsEdit.AddField(field);                  // setup field properties             fieldEdit.Name_2 = "Shape_Area";             fieldEdit.Type_2 = ESRI.ArcGIS.Geodatabase.esriFieldType.esriFieldTypeDouble;             //fieldEdit.IsNullable_2 = true;             //fieldEdit.Editable_2 = true;             fieldsEdit.AddField(field);                  fields = (ESRI.ArcGIS.Geodatabase.IFields)fieldsEdit;                        //// add OID field                  //fields.addField(createOIDField("ObjectID"));             //// add NAME field             //fields.addField(createNameField("name"));             //// add RASTER field             //fields.addField(createRasterField("raster", isByRef, rasterSpatRef));             //// add SHAPE field             //fields.addField(createShapeField("shape", geoSpatRef));                  return fields;         }      Any help would be appreciated.


title: How to use QGIS Accuracy Assessment plugin?

text: I am new user and i want perform accuracy assessment in qgis. I have 1 classified map and another for validation. I have installed the plugin. what am i doing after? Plugin: http://plugins.qgis.org/plugins/accassess/ Docu: https://bitbucket.org/jkibele/accassess/overview


title: How to hide features based selection using API for Silverlight?

text: my challenge is to, when i want to select feature on the map by the mouse the action must be zoom to the feature and all other features hidden in the same layer my application created by API for silverlight with the latest version. more details; i have 13 region in one polygon layer, our vendor required, when he want to select any region on the map by mouse, the action on the map must be zoom to the select feature and other features which "12" region must be hidden.


title: issues with updating data source, script runs but nothing happens, arcpy

text: I have seen many questions posted that are similar to mine but still havent found the answer so figured I would post my script. I want to change all the data sources found in all mxds in a location to point to one geodatabase. The mxds are currently referencing feature classes, shapefiles and rasters. I am wondering if there is an issue with the shapefiles that have now been moved to the gdb and no longer have the file extension.. even if so why wouldn't the feature classes sources change? Here is what I have, it runs, but there are no changes to the data sources that I can see. UPDATE: I've now modified it to use lyr.ReplaceDataSource, which works it changing the data source for the first shapefile it encounters in the mxds, but it doesn't change the source for the other feature classes and rasters. UPDATE 2: Im still struggling to find the answer here. I now am trying to use findandreplaceworkspacepaths on each mxd, and again the script runs, but the data sources don't change. I am printng out the old and new paths in the script and they are printing correctly.. so what is happening?? updated code below:               import arcpy, os     mxd_path = r'C:\GIS\Assimilate_Mxd_Data\Assimilate_Mxd_Data_TESTING\Input Mxd'     new_path = r'C:\GIS\Assimilate_Mxd_Data\Assimilate_Mxd_Data_TESTING\Output Data\test_output.gdb'     mxd_list = []          for dirpath, dirnames, filenames in os.walk(mxd_path):         for filename in filenames:             fullPath = os.path.join(dirpath, filename)             basename, extension = os.path.splitext(filename)             if extension == ".mxd":                 each_mxd = fullPath                 mxd_list.append(each_mxd)     for m in mxd_list:         mxd = arcpy.mapping.MapDocument(m)         for lyr in arcpy.mapping.ListLayers(mxd):             if lyr.supports("DATASOURCE"):                 old_path = lyr.dataSource                 mxd.findAndReplaceWorkspacePaths(old_path,new_path,validate=True)                 print old_path,new_path     mxd.save()     del mxd


title: Is it easy to digitize in ArcGIS using the Wacom tablet?

text: I do just enough digitizing (and touching up/ drawing when post-processing) to perhaps warrant the purchase of a Wacom tablet. Just wondering how well it plays with ArcGIS. What are your experiences? Is there a learning curve. Can the tablet and pen replace the mouse all together?


title: Help with a Simple Code to Read a Metadata from a Raster File

text: Please help me with a python code to read metadata from a Raster file. Thanks


title: Select by attributes - records with the largest attribute value

text: Using 'Select By Attributes' how can I select records from a shapefile database table where a numeric field contains the largest value. Need to use MAX function or something like that but can't get the syntax right.


title: What URL should i use?

text: By using the code below, we can query featureSet to ArcGIS Sample Server.               queryTask = new esri.arcgis.gmaps.QueryTask("http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/3");      I have published data on Geoserver by REST service.   1. Can I make this type of query (QueryTask) to Geoserver ?   2. What URL string should I use?  As I understood, the URL should be ("http://localhost:8080/geoserver/rest/"), but I don't know the complete URL. I mean what should I append with the URL after /rest/..... Can anybody explain the complete format of the URL?


title: Is it possible to draw direction arrows using SLD of GeoServer

text: I have one record in MySQL with LINESTRING field which contains list of points. Is it possible to draw direction arrows over whole line using SLD of GeoServer? Or i need split my one line?


title: How to export raster from PostGIS enabled DB as GTiff with gdal_translate

text: When I want to export a GTiff from my PostGIS DB with this code:               gdal_translate -of GTiff PG":host='localhost' port:'5432' dbname='landslide' user='postgres' password='pass' schema='public' table='dgm' mode='2'" test.tiff      I get the following error:               Illegal filename. character <'> not allowed.Illegal filename. character <'> not allowed.ERROR 4: `PG:host='localhost' port:'5432' dbname='landslide' user='postgres' password='pass' schema='public' table='dgm' mode='2'' does not exist in the file system, and is not recognised as a supported dataset name.          GDALOpen failed - 4 `PG:host='localhost' port:'5432' dbname='landslide' user='postgres' password='pass' schema='public' table='dgm' mode='2'' does not exist in the file system, and is not recognised as a supported dataset name.      I followed the instructions from gdal. Is the syntax incorrect? I tried to change the quotes, but no success.


title: Does Earth's magnetic field affect any datum definition?

text: I had an interview last week. One very interesting question was "does earth's magnetic field affect any datum definition?" I was sure it does not for geocentric datums. But I thought maybe **the geodetic datums** are affected by earth's magnetic field. But then, none of the physical parameters of a geodetic datum are related to Earth's magnetic field! And then I got confused! Can anybody give the correct answer?


title: Mapnik TextSymbolizer showing nothing above OpenStreetMap base layer

text: I'm trying to use Python to create a TextSymbolizer to display a trail name above the trail LineSymbolizer. Below this layer is the OSM base layer which contains its own data for that trail. I can display the line but I cannot display any text. I've run out of ideas, Google is proving fruitless so I'm looking for new ideas, if I may be so bold. I've compiled Mapnik with debug turned on only to find that there's no logging near the text symbolizer code. I had a go with gdb but I haven't spent enough time to workout how to debug from python into the Mapnik shared library. If anybody has any guidance for how to do this then I'd be grateful. I'm happy enough to debug but I'm not familiar with gdb. I also edited the osm.xml and took away the trails and trail names which might be interfering with the TextSymbolizer but I didn't see any difference. Relevant code:                   imgx = 400         imgy = 400              m = mapnik.Map(imgx,imgy)         m.buffer_size = 1000              workoutParams = dict(host='localhost',port='5432',user='postgres',password='password',dbname='postgis')         workoutParams['estimate_extent'] = False         workoutParams['extent']='-180,-89.99,180,89.99'         workoutTable = "(select id, geometry, 'really long trail name' as trail_name from pub_trail_segment where pub_trail_id = %s) as route" % idNumber         workoutParams['table']= workoutTable         workoutDatasource = PostGIS(**workoutParams)         workoutLayer = Layer('route','+proj=latlong +datum=WGS84')         workoutLayer.datasource = workoutDatasource         workoutStyle = mapnik.Style()         segmentStyleRule = mapnik.Rule()         segmentStyleRule.symbols.append(mapnik.LineSymbolizer(mapnik.Color('blue'),1)) #50 for print         segmentTextSymbolizer = mapnik.TextSymbolizer('trail_name', 'book-fonts', 20, mapnik.Color('blue'))         segmentTextSymbolizer.halo_radius = 1         segmentTextSymbolizer.label_placement = mapnik.label_placement.LINE_PLACEMENT         segmentTextSymbolizer.allow_overlap = True         segmentTextSymbolizer.avoid_edges = False              segmentStyleRule.symbols.append(segmentTextSymbolizer)         workoutStyle.rules.append(segmentStyleRule)         m.append_style('route',workoutStyle)         workoutLayer.styles.append('route')         m.layers.append(workoutLayer)         m.zoom_to_box(bbox)         im = mapnik.Image(imgx,imgy)         mapnik.render(m, im)      What am I doing wrong? Many thanks for any ideas Cros-posted to StackOverflow


title: Why is select by drawing a rectangle not working in ArcGIS?

text: I click on the select features tool on the toolbar and when I make a rectangle around the points I want to select it doesn't select them. Does anyone know what I may be doing wrong or how I may select points by making a shape around them? Thanks


title: How to determine the location that has the most intersections with other polygons in a certain distance/buffer?

text: I have a very scattered map with a lot of very differently shaped polygons. Now I want to determine the location where a circle polygon of a certain size intersects the most surface of these scattered polygons. It would be even better to be able to determine not only one best location but maybe even some more locations that are the most intersecting to get something like a heatmap. Here is an exemplary picture: ![enter image description here](http://i.stack.imgur.com/W8kbw.jpg)


title: Failing to find a location in Google Earth

text: I type a location into the Google Earth search bar and click search. Most of the time I get zoomed to the correct location, but on about 30% of occasions nothing happens. GE fails to take me to that location (fly to location) and simply refuses to zoom to the keyed in location. This occurs both when the place name is recognised by GE from the suggestions box: ![enter image description here](http://i.stack.imgur.com/IDxB7.jpg) And also when it is keyed in without a suggestion. I don't mind too much, as most of the time I do not have this problem. However, when it fails to work it can stall my work and I have to wait until GE decides to work again. There must be some kind of quick fix for this problem **Version 7.1.2.2041**


title: How to get the path of enterprise geodatabase stored in sql server?

text: I want to enter a geodatabase in a ListVersions tool.Output will be list of versions in a table. Do i need to give a connection file path for Enterprise_GDB: ___ _ or path of actual GDB. and how???? I am getting an error. Tool is not working successfully...


title: Road Miles Calculation without Google Maps

text: We have two ZIP codes we want to calculate the road mile difference. Google maps states you cannot use their API if you do not show a map. Is there another service that provides this type of calculation? I have searched Google but nothing has turned up useful. Currently we use lat/lon to get as the crow flies calculations.


title: How can I delete from attribute table in a dxf file in ArcGIS?

text: Hy everybody! I have tried delete one row from an attribute table in a dxf file, but I can't. ArcMap said I need to export that to a geodatabase or something similar. What exactly should I do? From step to step, please.


title: How to extract mountain peak points from a given digital elevation model?

text: For past couple of days, i m trying to determine how many mountain peaks are visible (line of sight) from a given location using POSTGIS. The data i have is DEM for Swiss and address locations from Swiss which are stored in a table in postgis database. I downloaded the mountain peak points from OSM(Open street Map), and then wrote a function in postgis to find if the peak is visible from a given address location using address data from the table, elevation raster(imported in the postgis db), and mountain peak data imported from OSM to POSTGIS... The Funtion is giving me output but the quality of the results are very poor as OSM mountain peaks are not proper and are not aligned to the mountain peak points in the elevation raster(problem with OSM data Quality). projections are matching of both the layers.. so now to solve this problem i need to get the exact mountain peaks. So i was just wondering if i can extract mountain peaks from elevation raster. with it co-ordinates n elevation.. which i can then used to calculate Line of sight... can anyone here please help me?? i m not sure if this is possible Thanks in advance


title: How to change the size of Hex bins?

text: I am creating a map in QGIS and I am using a hex bin layer. To create this layer I used the MMQGIS plugin and created a grid layer. However, the size of my hexagons that I created are too small, and I was wondering how I can change their size to make them bigger.


title: Is it possible to create DEMs from ORT format in QGIS?

text: I recently requested DEMs for my area of interest. The data has however arrived in ORT file format - how do I create DEMs from this data in QGIS? When I open the files in Notepad there are three columns which respectively contain x and y co-ordinates in m (Hartebeeshoek Lo31) and a third column containing the height value.


title: Displaying CartoDB visualization on mobile with "Where Am I" button

text: I want to display a CartoDb visualization on my mobile device, but also have an icon that shows where the device is and centers/zooms to it when I press it. I'm trying to find nearby features. If this isn't available via a web page, then perhaps via a PhoneGap plugin?


title: Populating combobox and zoom to selected field in box?

text:
I am very new to arcgis javascript .I am not able to populate combobox and zoom to the selected field in the box               <!DOCTYPE html>                        <meta http-equiv="Content-Type" content="text/html; charset=utf-8">         <meta name="viewport" content="initial-scale=1, maximum-scale=1,user-scalable=no">         <title>Feature Layer Only Map</title>         <link rel="stylesheet" href="http://desktop101.ksrsac.com/KSRSAC/arcgis_js_api/library/3.9/3.9/js/dojo/dijit/themes/claro/claro.css">         <link rel="stylesheet" href="http://js.arcgis.com/3.9/js/esri/css/esri.css">

"Tank Information System"

                         Content for the first tab

                         Content for the second tab                         <select id="mySelect" dojoType="ComboBox" style="width:200px;font-size:16px;"               autoComplete="true"                  forceValidOption="false"                  value="Select Zoning Type"                  onchange="applyLayerDef(this)" >                         </select>


title: Alternative backup File or Personal Geodatabase

text: I'm wondering if there is a way around or an alternative way to backup a File or Personal Geodatabase. For example: If I have a very large dataset within the File or Personal Geodatabase and I want to backup it up to several media storage devices such as External Hard drive, flash drive, DVD etc. Can this process be automated instead of performing it manually. I work under the Government, so I may be limited to what I can do here in the office.


title: GeoReference tiff image into the Google map?

text: I want to Geo-reference a tiff image contain a Map for a specific region, I need to add this region(image) and Geo-reference it on the Google map, SO I can see on the map for the region I select the image which I have inserted. the image also need information such as the Latitude and Longitude of corners and size, and the pixel per meter measurement. How to do the correct Geo-reference for that ?


title: Geo locate own view pic on google map

text: I need to use Google map( GMap API) on C#, and I need to import my own picture for some place I want (with high resolution), plus it has its information on map such as : width,high, size of pixel , meter per pixel and its latitude and longitude. then I will see my view of this region on the map.


title: How do I add Garmin .img files to QGIS?

text: I downloaded 10 .img files from garmin.openstreetmap.nl which are supposed to contain the topo tiles for the area I would like to add in qgis. I've been unable to find any instructions on how to use the files in qgis. Can someone explain or point me to where instructions are available to accomplish this? I need to apply contour lines to a map of Thailand in qgis. It appears the img tiles I currently have are unusable to achieve this. Where can I download whatever file type is necessary containing the contour lines for Thailand and by what steps, using what application can that file be converted to one which can be added as a layer to qgis?


title: How to get Coordinates values from the latitude and lontitude

text: I having the form, i have two fields, Latitude and longtitute, if i put the value from here, i need automatically mapping on the map. I want to get the value from Lat and Long, it must be generate this format `{ "type": "Point", "coordinates": [3427927.775 5793243.975] },` How to do this, Help me i'm new for this.


title: How to set and configure a value relation edittype in QGIS 2.0?

text: I am trying to set the edit type of a vector layer to use a value relation. The Python commands below are processed without errors or warnings, but the when I check the properties in QGIS, only the EditType is set, but the data settings are empty. I use:               vlayer.setEditType(12, QgsVectorLayer.ValueRelation)     vlayer.ValueRelationData( 'releatedtable20131027113855772', 'key', 'value',False, True)      I got the name of the related table by inspecting:               ml=QgsMapLayerRegistry.instance().mapLayers()     dict(ml)      My guess is that I use the wrong string for referencing the layername. I have two questions: How to get the layername string to be used in ValueRelationData? Should I add further code to get it all working?


title: converting a distance to a map distance

text: I have a point and a polyline. I need to find the perpendicular distance between each segment of the polyline and the point. By each segemnt of the polyline I mean that, each polyline may contain more than one line and I need to find every line of the polyline and find the perpendicular distance. I know its really easy using ICurve.QueryPointAndDistance Method. But to use this method, I need to create polylines from each segment of the actual polyline. But it is really time consuming for more than 13k polyline features. So I tried to create line from every segment of the polyline. And calculate perpendicular distance from a point to the line. but this distance is not same as the distance I get from _ICurve.QueryPointAndDistance_ Method's _distanceFromCurve_. ![image is attached for better understanding](http://i.stack.imgur.com/g2BZX.png) _**Question:_** Is there any better option to calculate perpendicular distance that will be same as the `ICurve.QueryPointAndDistance` Method's _distanceFromCurve_? Or what is the reason I am not getting the same result as `ICurve.QueryPointAndDistance` Method's _distanceFromCurve_? I am using arcObjects and c#. the input feature layer is always projected. So I am getting `X`, `Y` value in meter/feet. For calculating perpendicular distance I am using following method:               private double GetPerpendicularDistance(Line line, Point point)     {         var a = line.Point1.Y - line.Point2.Y;         var b = line.Point2.X - line.Point1.X;         var c = (line.Point1.X * line.Point2.Y) - (line.Point2.X * line.Point1.Y);              var x = point.X;         var y = point.Y;              var d1 = Math.Abs((a * x) + (b * y) + c);         var d2 = Math.Sqrt((a * a) + (b * b));         var d = d1 / d2;         return d;     }      Line is an internal class.               internal class Line     {         internal Point Point1 { set; get; }         internal Point Point2 { set; get; }     }      Point is another internal class.               internal class Point      {         internal double X { set; get; }         internal double Y { set; get; }     }


title: Invalid references created by pgr_createTopology

text: I am trying pgr_createTopology on a simple network and it is leaving invalid node references in the `source` and `target` columns. Apparently the segments table is not being updated with the new identifiers created in the `_vertices_pgr` table. An example:               # \d esch.roads                    Table "esch.roads"      Column |           Type            | Modifiers      --------+---------------------------+-----------      id     | numeric(10,0)             | not null      source | smallint                  |       target | smallint                  |       geom   | geometry(MultiLineString) |      Indexes:         "roads_pkey" PRIMARY KEY, btree (id)          # SELECT pgr_createTopology('esch.roads', 0.0001, 'geom', 'id', 'source', 'target', 'true');     NOTICE:  PROCESSING:     NOTICE:  pgr_createTopology('esch.roads',0.0001,'geom','id','source','target','true')     NOTICE:  Performing checks, pelase wait .....     NOTICE:  Creating Topology, Please wait...     NOTICE:  -------------> TOPOLOGY CREATED FOR  8 edges     NOTICE:  Rows with NULL geometry or NULL id: 0     NOTICE:  Vertices table for table esch.roads is: esch.roads_vertices_pgr     NOTICE:  ----------------------------------------------      pgr_createtopology      --------------------      OK     (1 row)          # SELECT id FROM esch.roads_vertices_pgr WHERE id NOT IN (SELECT source FROM esch.roads);      id      ----       2     (1 row)          # SELECT id FROM esch.roads_vertices_pgr WHERE id NOT IN (SELECT target FROM esch.roads);      id      ----       6     (1 row)      This seems wrong to me, or am I missing something?


title: Infowindow selection change

text: I am using infowindow to display the results of identifyTask that has two layers from my map service so it returns multiple features.   The infowindow has more button in each feature it identifies so that when I click the more button I get the data of current feature. How can I get the current selected feature? I am using this example.   I want to get incase the popup identifies more than one feature to get current feature How can this be done.


title: How to get information about only one data point?

text:
I have a vector-point shapefile. The data are not dense (about 1 x 1 degree). Lets say that every point includes value of average wind speed. I would like to create a interactive web-map so that every time a user click on it, he/she get information about the point of interest (POI). ![enter image description here](http://i.stack.imgur.com/kYzME.png) This map includes 3 layers   1. `bluemarble` base layer   2. WMS layer of points with data   3. TIFF layer created from wind speed data (red - high, blue - low) I have two problems. **1) to many points clicked** when the map is zoomed-out it happens that user's click covers more then 1 point of the vector so that GetFeatureInfo returns too much info. In the `event.text` I got info of the TIFF raster clicked pixel, GREEN, RED, BLUE, ALPHA_BAND and values of all the clicked points:                fid	       long	       lat	       codetec	       modelno	       tecdiff

       allinfo.122239	        -4.000000000000000	       4.000000000000000	        27.210000000000001	             2	       2.150000000000000

       allinfo.122240	       -4.000000000000000	       3.700000000000000	       27.210000000000001	       2	       2.150000000000000

       allinfo.122241	      -4.000000000000000	      3.300000000000000	      27.030000000000001	      2	      2.250000000000000

       allinfo.122242	       -4.000000000000000	       3.000000000000000	       27.030000000000001	       2	       2.250000000000000
       and some more, but you got the point. **2) no point clicked when zoomed-in** When it is zoomed-in too much the user obviously does not click any point of the vector layer. So there is no info about the vector layer only the TIFF raster colour band info. * * * I would say that the 1) problem can be solved by formatting the GetFeatureInfo response but I am lost in the problem 2). Is there a way how to force to get info from the nearest point? Or show info always from one point only? Should I change the vector layer? I tried to use create Voronoi polygons but I did not succeed to produce any layer (and there two links : Q1 and Q2 suggest that that does not work that well in QGIS anyway). Is it possible to retrieve the value from the raster colour band?


title: Can we create relationship classes based on a geometric network?

text: I have a point Feature Class and a polyline Feature class which I connected in a geometric network. I also need to connect them logically. In my theoretical model, an edge is defined by the two points at its end. So I need a Many-to Many relationship. An edge can have two points and a point can have two edges. Connecting them in a relationship requires the maintenance of unique ID's etc. So I was wondering is there a possibility that based on their geographical relationship of points and line to create as well a logical relationship in the database, and save me from the whole process of giving them ID's, connecting them, etc? It's probably a rather general and abstract question, but I guess some of you might understand what I mean. Thank you! I use ArcGIS 10.0.


title: How to use centroid option in Field Calculator in QGis

text: I want to use Centroid option in field calculator in QGis .How to access it?


title: How to make a field "Unique ID field" in ArcGIS?

text: I need to perform an OLS (Ordinary Least Squares) but ArcGIS asks for a unique ID field. I have a field which would be the right candidate, it's a FIPS number (12 chars), which I have made a number, and text in Excel to no avail. OLS does not work without it... How can I make it a legal Unique ID field?


title: Conceptually, what's the ideal way to develop detour maps for a large number of nodes

text: I have a road layer of roughly 30,000 links. On these links, I have a total of 15,000 nodes that are sorted in a point layer. I'm trying to build a detour distance and routing around each asset in case that asset fails. I have access to ArcGIS (not with network analyst) and QGIS latest build. I have access to PostGIS and pgRouting but I'd rather stay perhaps in network tools whether it's ArcPy, QGIS python or even straight python with shapely and netowkrX Currently, I'm in the solution design stage. I'm curious as to how you'd recommend I build a "rerouting" engine that can efficiently output `detour_distance` and `list_of_routes_of_detour` **How would you recommend I build such a system?**


title: Routing using .MP file format

text: I am having .MP file for a particular map location, the .MP file consists of points and polylines. Now I want to find the shortest path between two points. What is perfect and simple method to find shortest path. Also please provide me good tutorial, tools.. on this... Thank in advance.


title: Geoserver transparency of overlapping features

text: I'm serving a layer via WMS that has multiple Polygons with transparency set to 30%. The polygons can overlap each other resulting in darker areas where the transparencies combine. Is it possible to prevent this? UPDATE: Transparency is set via `<CssParameter name="fill- opacity">0.5</CssParameter>` of `<Fill>` element in SLD definition. The whole SLD looks like this:               <?xml version="1.0" encoding="UTF-8"?>     <StyledLayerDescriptor version="1.0.0" xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc"       xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"       xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.0.0/StyledLayerDescriptor.xsd">       <NamedLayer>         <Name>plochaQ5ZU</Name>         <UserStyle>           <Name>Q5ZU</Name>           <FeatureTypeStyle>                  <Rule>               <Name>Q5ZU</Name>               <Title>Q5ZU</Title>               <PolygonSymbolizer>                 <Fill>                   <CssParameter name="fill">#0070D4</CssParameter>                   <CssParameter name="fill-opacity">0.5</CssParameter>                 </Fill>               </PolygonSymbolizer>             </Rule>           </FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>      I'm using OpenLayers to display the layer.


title: Can the software be installed silently?

text: Can QGIS be silently installed and uninstalled? We are looking to deploy this to several users for evaluation.


title: Test of ORACLE DB functions converting lat-long to utm and mgrs with real data known to be correct?

text: I have rewritten some Delphi function to ORACLE DB functions for converting from lat-long to utm and mgrs. Anyone care to validate the output with real data they know is correct? package specification:               create or replace package gedaco as          function MGRS(lat in number, Lon in number,  a in number, InverseFlattening in number,                         Coding in number,  Digits in number) return varchar2;         Function MGRSLatZone(lat in number) return varchar2;         function SquareID(UTMzn in number, Northing in number, Easting in number,                            Coding in number) return varchar2;         function UTM(lat in number, Lon in number, a in number, InverseFlattening in number)                                     return varchar2;         Function UTMX(UTMs in varchar2) return number;         Function UTMY(UTMs in varchar2) return number;         function UTMZone(lat in number, Lon in number) return number;     end gedaco;      package body:               create or replace package body gedaco as           function MGRS(lat in number, Lon in number,  a in number, InverseFlattening in number,                   Coding in number,  Digits in number) return varchar2 is        result varchar2(32);       UTMs1 varchar2(32);       E1 number;       N1 number;       Zn number;       Lzn varchar2(32);       Sq varchar2(32);     begin       UTMs1 := UTM(lat, Lon, a, InverseFlattening) ;       E1 := UTMX(UTMs1);       N1 := UTMY(UTMs1);       Zn := UTMZone(lat, Lon);       Lzn := MGRSLatZone(lat);       Sq := SquareID(Zn, N1, E1, Coding);       result := replace(           to_char(Zn,'00') || LZn || Sq ||           to_char(round(E1 - 100000 * trunc(E1/100000)),'00000') ||           to_char(round(N1 - 100000 * trunc(N1/100000)),'00000')                ,' ', '');       return result;     end MGRS;               Function MGRSLatZone(lat in number) return varchar2 is          result varchar2(1);           GridZones CONSTANT varchar2(20) := 'CDEFGHJKLMNPQRSTUVW';     begin       If (lat >= 72) Then Result := 'X';       Else Result := substr(GridZones, Trunc((lat + 88) / 8), 1);       End If;       return result;     end MGRSLatZone;          function SquareID(UTMzn in number, Northing in number, Easting in number,                        Coding in number) return varchar2 is         result varchar2(32);         N number;         E number;         ZoneSet number;         Col varchar2(32);         Rov varchar2(32);         Col1 CONSTANT varchar2(20) := 'ABCDEFGH';         Col2 CONSTANT varchar2(20) := 'JKLMNPQR';         Col3 CONSTANT varchar2(20) := 'STUVWXYZ' ;         Row1 CONSTANT varchar2(20) := 'ABCDEFGHJKLMNPQRSTUV';         Row2 CONSTANT varchar2(20) := 'FGHJKLMNPQRSTUVABCDE';         Row3 CONSTANT varchar2(20) := 'LMNPQRSTUVABCDEFGHJK';         Row4 CONSTANT varchar2(20) := 'RSTUVABCDEFGHJKLMNPQ';        begin       N := Trunc(Northing / 100000);       N := N - 20 * Trunc(N / 20);       E := Trunc(Easting / 100000);       ZoneSet := UTMzn - 6 * Trunc(UTMzn / 6);       If ((ZoneSet = 1) Or (ZoneSet = 4)) Then          Col := SubStr(Col1, E, 1);       End If;       If ((ZoneSet = 2) Or (ZoneSet = 5)) Then          Col := SubStr(Col2, E, 1);       End If;       If ((ZoneSet = 3) Or (ZoneSet = 0)) Then          Col := SubStr(Col3, E, 1);       End If;            ZoneSet := ZoneSet - 2 * Trunc(ZoneSet / 2);       If ((Coding = 1) And (ZoneSet = 1)) Then         Rov := SubStr(Row1, N + 1, 1);       End If;       If ((Coding = 1) And (ZoneSet = 0)) Then         Rov := SubStr(Row2, N + 1, 1);       End If;       If ((Coding = 2) And (ZoneSet = 1)) Then         Rov := SubStr(Row3, N + 1, 1);       End If;       If ((Coding = 2) And (ZoneSet = 0)) Then         Rov := SubStr(Row4, N + 1, 1);       End If;            Result:= Col || Rov;       return result;     end SquareId;               function UTM(lat in number, Lon in number, a in number, InverseFlattening in number)       return varchar2       is result varchar2(320);       ZoneWidth CONSTANT number := 6;       CentralScaleFactor CONSTANT number := 0.9996;       Zone1CentralMeridian CONSTANT number := -177;       Zone0WestMeridian number;       Zone0CentralMeridian number;       FalseEasting CONSTANT number := 500000;            Pi number;       SemiMajorAxis number;       Flattening number; Eccent2 number; Eccent4 number; Eccent6 number;       A0 number; A2 number; A4 number; A6 number;       LatRad  number;       LonRad  number;       Sin1Lat  number; Sin2Lat  number; Sin4Lat  number; Sin6Lat  number;       Rho  number;       Nu  number;       Psi  number; Psi2  number; Psi3  number; Psi4   number;       CosLat  number; CosLat2  number; CosLat3  number; CosLat4  number; CosLat5   number;       CosLat6  number; CosLat7    number;       TanLat  number; TanLat2  number; TanLat4  number; TanLat6    number;       DifLon  number; DifLon2  number; DifLon3  number; DifLon4  number; DifLon5   number;       DifLon6 number; DifLon7  number; DifLon8   number;       DistOverMeridian    number;       Zone  number;       CentralMeridian  Integer;       East1  number; East2  number; East3  number; East4  number;       North1  number; North2  number; North3  number; North4  number;       X  number;       Y  number;       Hemi  varchar2(1);       FalseNorthing  number;          begin     Zone0WestMeridian := Zone1CentralMeridian - (1.5 * ZoneWidth);     Zone0CentralMeridian := Zone0WestMeridian + ZoneWidth / 2;     Pi := 3.141592653589793238462643383279502884197169399375105820974944592307816406;               SemiMajorAxis := 1000 * a  ;          Flattening := 1.0 / InverseFlattening   ;     Eccent2 := 2.0 * Flattening - (Flattening * Flattening);     Eccent4 := Eccent2 * Eccent2   ;     Eccent6 := Eccent2 * Eccent4 ;     A0 := 1 - (Eccent2 / 4.0) - ((3 * Eccent4) / 64.0) - ((5.0 * Eccent6) / 256.0);     A2 := (3.0 / 8.0) * (Eccent2 + (Eccent4 / 4.0) + ((15.0 * Eccent6) / 128.0)) ;     A4 := (15 / 256) * (Eccent4 + ((3.0 * Eccent6) / 4.0));     A6 := (35.0 * Eccent6) / 3072.0 ;       --  ' Parameters to radians         LatRad := lat / 180 * Pi;         LonRad := Lon / 180 * Pi ;                 --  'Sin of latitude and its multiples         Sin1Lat := sIn(LatRad) ;         Sin2Lat := sIn(2 * LatRad) ;         Sin4Lat := sIn(4 * LatRad);         Sin6Lat := sIn(6 * LatRad);            --  'Meridian Distance         DistOverMeridian := SemiMajorAxis *                              (A0 * LatRad - A2 * Sin2Lat + A4 * Sin4Lat - A6 * Sin6Lat);                 --  'Radii of Curvature         Rho := SemiMajorAxis * (1 - Eccent2) /Power( (1 -                 (Eccent2 * Sin1Lat * Sin1Lat)) , 1.5);         Nu := SemiMajorAxis /Power( (1 - (Eccent2 * Sin1Lat * Sin1Lat)) , 0.5);         Psi := Nu / Rho  ;         Psi2 := Psi * Psi ;         Psi3 := Psi * Psi2;         Psi4 := Psi * Psi3  ;            --  'Powers of cos latitude         CosLat := Cos(LatRad);         CosLat2 := CosLat * CosLat  ;         CosLat3 := CosLat * CosLat2 ;         CosLat4 := CosLat * CosLat3 ;         CosLat5 := CosLat * CosLat4 ;         CosLat6 := CosLat * CosLat5 ;         CosLat7 := CosLat * CosLat6 ;               --    'Powers of tan latitude         TanLat := Tan(LatRad) ;         TanLat2 := TanLat * TanLat ;         TanLat4 := TanLat2 * TanLat2  ;         TanLat6 := TanLat2 * TanLat4  ;           --   'Zone      --   'Zone := Int((Lon - Zone0WestMeridian) / ZoneWidth)         Zone := UTMZone(lat, Lon)   ;                        CentralMeridian := Trunc((Zone * ZoneWidth) + Zone0CentralMeridian ) ;         DifLon := (Lon - CentralMeridian) / 180 * Pi    ;         DifLon2 := DifLon * DifLon  ;         DifLon3 := DifLon * DifLon2 ;         DifLon4 := DifLon * DifLon3 ;         DifLon5 := DifLon * DifLon4 ;         DifLon6 := DifLon * DifLon5 ;         DifLon7 := DifLon * DifLon6 ;         DifLon8 := DifLon * DifLon7 ;              East1 := DifLon * CosLat  ;         East2 := DifLon3 * CosLat3 * (Psi - TanLat2) / 6.0;         East3 := DifLon5 * CosLat5 * (4.0 * Psi3 * (1.0 - 6.0 * TanLat2) + Psi2 *          (1.0 + 8.0 * TanLat2) -Psi * (2.0 * TanLat2) + TanLat4) / 120.0;         East4 := DifLon7 * CosLat7 * (61.0 - 479.0 * TanLat2 + 179.0 * TanLat4 - TanLat6)          / 5040.0  ;         X := CentralScaleFactor * Nu * (East1 + East2 + East3 + East4) + FalseEasting  ;              If (lat >= 0) Then           Hemi := 'N';           FalseNorthing := 0;         Else           Hemi := 'S';           FalseNorthing := 10000000;         end if;              North1 := Sin1Lat * DifLon2 * CosLat / 2.0 ;         North2 := Sin1Lat * DifLon4 * CosLat3 * (4.0 * Psi2 + Psi - TanLat2) / 24.0 ;         North3 := Sin1Lat * DifLon6 * CosLat5 * (8.0 * Psi4 * (11.0 - 24.0 * TanLat2)                   - 28.0 * Psi3 * (1.0 - 6.0 * TanLat2) +          Psi2 * (1.0 - 32.0 * TanLat2) - Psi * (2.0 * TanLat2) + TanLat4) / 720;         North4 := Sin1Lat * DifLon8 * CosLat7 * (1385 - 3111 * TanLat2 + 543 *                    TanLat4 - TanLat6) / 40320.0 ;         Y := CentralScaleFactor * (DistOverMeridian + Nu *              (North1 + North2 + North3 + North4)) + FalseNorthing;            Result := Zone || Hemi || ' ' ||                  to_char(round(X, 3),'0000000.000') ||                  to_char(round(Y, 3),'0000000.000');       return result;     End UTM;               Function UTMX(UTMs in varchar2) return number         is result number;     begin       Result := to_number(substr(UTMs, 6, 11));       return result;     End UTMX;               Function UTMY(UTMs in varchar2) return number         is result number;     begin       Result := to_number(substr(UTMs, 18, 11));       return result;     End UTMY;               function UTMZone(lat in number, Lon in number) return number is       result number;       UTMZone number;       e number;       d number;       ZoneWidth CONSTANT number := 6;       Zone1CentralMeridian CONSTANT number := -177;       Zone0WestMeridian number;     begin       Zone0WestMeridian := Zone1CentralMeridian - (1.5 * ZoneWidth);       d:=ZoneWidth;             UTMZone := Trunc((lon - Zone0WestMeridian) / d);         --Special Cases for Norway & Svalbard         CASE          WHEN (lat > 55) AND (UTMZone = 31) AND (lat < 64) AND (lon >  2) THEN UTMZone := 32;         WHEN (lat > 71) AND (UTMZone = 32) AND (lon <  9) THEN UTMZone := 31;         WHEN (lat > 71) AND (UTMZone = 32) AND (lon >  8) THEN UTMZone := 33;         WHEN (lat > 71) AND (UTMZone = 34) AND (lon < 21) THEN UTMZone := 33;         WHEN (lat > 71) AND (UTMZone = 34) AND (lon > 20) THEN UTMZone := 35;          WHEN (lat > 71) AND (UTMZone = 36) AND (lon < 33) THEN UTMZone := 35;         WHEN (lat > 71) AND (UTMZone = 36) AND (lon > 32) THEN UTMZone := 37;         ELSE UTMZone := UTMZone;           END CASE;            Result := UTMZone;       return result;     end UTMZone;          end gedaco;      Function is used MGRS(:latitude, :longitude, 6378.137, 298.2572236, 1, 5) for WGS84 with 5 digits precision.   list of datums:               Datum           Radius          InverseFlattening     'WGS84',        6378.137,       298.2572236     'NAD27',        6378.2064,      294.9786982     'NAD83',        6378.137,       298.2572221     'WGS66',        6378.145,       298.25     'GRS67',        6378.16,        298.2472     'IAU68',        6378.16,        298.2472     'WGS72',        6378.135,       298.26     'Clarke66',     6378.2064,      294.9786982     'GRS80',        6378.137,       298.2572221     'Krasovsky',    6378.2064,      298.3     'Bessel',       6377.397155,    299.1528128      Just to be completely specific, I have of course done some random checks with Earth Point, but I am asking for someone to control my functions with a sustainable amount of data. A proper reply to my question can simply be test with 10 000 records without errors found, and if errors are found I am of course interested in the location not calculated correctly if that is possible.


title: How to change the color for Openlayer Layerswitcher?

text: I tried to change the layerswitcher color. But its not changed. Here is my code. I tried in javascript also osMap.addControl(new OpenLayers.Control.LayerSwitcher({ 'activeColor': "white", 'fontColor': "black" })); but there is no effect.               .olControlLayerSwitcher     {         position: absolute;         top: 25px;         right: 0;         width: 20em;         font-family: sans-serif;         font-weight: bold;         margin-top: 3px;         margin-left: 3px;         margin-bottom: 3px;         font-size: smaller;         color: white;         background-color: transparent;         z-index: 10000;     }     .olControlLayerSwitcher .layersDiv     {         padding-top: 5px;         padding-left: 10px;         padding-bottom: 5px;         padding-right: 10px;         background-color: #CCCCCC;     }


title: Table to dBase tool failure when converting .xls files

text: Whenever I load my Excel 2010 worksheets into ArcMap 10 Table to dBASE tool, it fails, returning the error: "Failed to Copy Rows" The Copy Rows tool works just fine for me on a single document. There are no special characters or spaces in my file paths, workbook or worksheet names, or in column headers. Has anybody else out there had an issue with this tool? Or is there a simple fix for this that I am missing?


title: How can I fix a broken gdal-tools plugin in QGIS?

text: I have QGIS installed on my 2011 macbook pro and I recently upgraded my OS to 10.9. After the upgrade, I noticed that QGIS gives me this error (see attachment 1) saying that the plugin is not installed (but it is; the second screenshot shows that it is "local only"). I suppose I can re-install the plugin but I'm hoping that it isn't necessary. Can someone provide some guidance as to how to get the plugin working again? Is it going to require that I completely remove the plugin and re-install? ![Attachment 1 showing error message at startup](http://i.stack.imgur.com/mN7s7.jpg) ![Attachment 2 showing dialog from plugins menu](http://i.stack.imgur.com/uyQTZ.jpg)


title: arcgis10: export eventlayer (from querylayer) to ArcSDE

text: I am not using arcmap but I am creating a arcgis desktop console app. I am pulling an external SQL table in as a querylayer. This table has lat & long so I created an eventlayer (or IXYEventSource xyEventSource). Now I want to export it to SDE as a featureclass but now I am stuck which route to take. Eventually, I want to put this exe in windows scheduler. Any ideas on how to export to SDE? Thanks, saim