3.30.2011

Climate change in the Philippines

Wondering about the recent temperature changes in the Philippines for the whole year round.
Here's a climate graph of the Philippines.



Philippines Climate graph contributed by climatetemp.info

How to view shapefiles in Google Earth

I have tried it in ArcGIS 9.3.

In Arc Toolbox choose Conversion tools > To KML then choose the layer that you want to convert then choose scale as 1 if your data is currently in WGS84 projection. It will be saved as a zipped kml file (.kmz) because you chose to save a layer of your map document. Just double click the output kmz file and there you go.

Here's a screenshot of my sample shapefile (Isabela, Philippines) as viewed from Google Earth.


In one of my graduate course project, we have studied dengue and identify high, medium and low risk areas in  Quezon City, Philippines. We have also included the Google Earth upload of our output layers as part of our methodology. Here's a sample shot showing how the fields and attributes are preserved after exporting the shapefile to kmz.

Dengue Hotspots in Quezon City, Philippines
So there you go. I'll be posting too on how to work the other way around; from KML to SHP next time. Enjoy googling (^_^).

3.17.2011

Instant information experience for GPS using Qwiki

A "Qwiki" is a short, interactive story: a drastically improved information experience provided via interactive video.

This is a great tool for research. It has an interactive interface highly visual in nature with short notes that is really comprehensive and smart.

I have tried querying on GPS and this is the result:


View Global Positioning System and over 3,000,000 other topics on Qwiki.


Qwiki is really super nice with the instant upload on my favorite social networking
sites. You must try it too!

3.16.2011

PYTHON CODE THAT ADDS A NEW FIELD TO A SHAPEFILE AND SPECIFIES ITS CONTENT

Arcpy class is new to ArcGIS 10.0. ArcPy is a site-package that builds on (and is a successor to) the successful arcgisscripting module. Its goal is to create the cornerstone for a useful and productive way to perform geographic data analysis, data conversion, data management, and map automation with Python.

Here's a code that will add a new field to a shapefile and then specify the contents of the added field. In this case, all areas greater than 1000 square meters will be labelled "YES" meaning it is subdividable otherwise it will be "NOO".

**********************************************************************************
import arcpy
def cursor():
    inFeatures = GetParameterAsText(0)
    subdvdable = GetParameterAsText(1)
    #inFeatures = "C:\\Users\\RS Lab\\Desktop\\222\\vector\\parcel_area.shp"
    arcpy.AddField_management(inFeatures, subdvdable, "TEXT")
    rows = arcpy.UpdateCursor(inFeatures)
    for row in rows:
        if row.Shape_Area > 1000:
            row.subdvdable = "YES"
        else: row.subdvdable = "NOO"
     rows.updateRow(row)

    del row,rows

**********************************************************************************

To implement this, a new toolbox should be created then the parameters will be set using a GUI. For the inFeatures parameter, it can be specified as Feature Class. The next parameter subdvdable will be of String type which will specify the name of the field to be added.

When run, the specified shapefile will already contain the added field and the specified content based on the logic set.