Changeset 2922 for production

Timestamp:

May 19, 2006, 2:47:30 PM (19 years ago)

Author:

sexton

Message:

trying to include location of indigeneous communities (need to check projection)

Location:

production

Files:

• onslow_2006/damage.tex (modified) (1 diff)
• onslow_2006/onslow_communities.png (added)
• onslow_2006/plot_data_extent.py (modified) (2 diffs)
• onslow_2006/project.py (modified) (1 diff)
• pt_hedland_2006/damage.tex (modified) (1 diff)
• pt_hedland_2006/project.py (modified) (2 diffs)
• pt_hedland_2006/pt_hedland_communities.png (added)
• pt_hedland_2006/results.tex (modified) (2 diffs)

production/onslow_2006/damage.tex

@@ -1 +1 @@
 
 this will be a discussion of the damage modelling
+
+need to talk about potential impact on indigeneous
+communities here
+
+
+\begin{figure}[hbt]
+
+  \centerline{ \includegraphics[width=100mm, height=75mm]{onslow_communities.png}}
+
+  \caption{Location of indigeneous communities in study area.}
+  \label{fig:communities}
+\end{figure}

production/onslow_2006/plot_data_extent.py

@@ -11 +11 @@
 ion()
 
-plot_data = True
+plot_data = False
 plot_with_poly = True
+
+def read_community_file(filename1, filename2, polygon):
+    from coordinate_transforms.redfearn import redfearn
+    from utilities.polygon import inside_polygon
+    fid1 = open(filename1)
+    lines = fid1.readlines()
+    fid1.close()
+    xlocs = []
+    ylocs = []
+    community = []
+    population = []
+    line1 = lines[0]
+    line11 = line1.split(',')
+    fid2 = open(filename2,'w')
+    s = 'Easting, Northing, Community, Population\n'
+    fid2.write(s)
+    for i in range(len(line11)):
+        if line11[i].strip('\n').strip(' ') == 'LATITUDE': lat_index = i
+        if line11[i].strip('\n').strip(' ') == 'LONGITUDE': lon_index = i
+        if line11[i].strip('\n').strip(' ') == 'COMMUNITY': name_index = i
+        if line11[i].strip('\n').strip(' ') == 'POPULATION': pop_index = i
+    for line in lines[1:]:
+        fields = line.split(',')
+        lat = float(fields[lat_index])
+        lon = float(fields[lon_index])
+        z, easting, northing = redfearn(lat,lon)
+        v = inside_polygon([easting,northing],polygon)
+        if v == True:
+            xlocs.append(easting)
+            ylocs.append(northing)
+            community = fields[name_index]
+            population = fields[pop_index]
+            s = '%f, %f, %s, %s\n' %(easting, northing, community, population)
+            fid2.write(s)
+            print community, easting, northing
+    
+    return xlocs, ylocs, community, population
+
+def read_file(filename1):
+    fid1 = open(filename1)
+    lines = fid1.readlines()
+    fid1.close()
+    x = []
+    y = []
+    line1 = lines[0]
+    for line in lines[1:]:
+        fields = line.split(',')
+        x.append(float(fields[0]))
+        y.append(float(fields[1]))
+    
+    return x, y
 
 if plot_data == True:
@@ -52 +103 @@
                         figname,
                         verbose = True)
+    figure(4)
+    x1, y1 = poly_xy(project.polyAll)
+    figname = 'onslow_communities'
+    x, y, community, population= \
+           read_community_file(project.community_filename, \
+                               project.community_scenario, project.polyAll)
+    #coastx, coasty = read_file(project.coast_filename)
+    #plot(x1,y1,'r-',x,y,'r+', coastx,coasty,'g.')
+    plot(x1,y1,'r-',x,y,'r+')
+    xlabel('Eastings')
+    ylabel('Northings')
+    savefig(figname)
 
 close('all')

production/onslow_2006/project.py

@@ -56 +56 @@
 #for MOST
 gauge_filename = gaugedir + 'onslow_gauges1.xya'
-
-
+buildings_filename = gaugedir + 'all_bld_ind.csv'
+community_filename = gaugedir + 'CHINS_v2.csv'
+community_scenario = gaugedir + 'community_onslow.csv'
 gaugetimeseries = gaugedir + 'onslow'
 

production/pt_hedland_2006/damage.tex

@@ -1 +1 @@
 
 this will be a discussion of the damage modelling
+
+need to talk about potential impact on indigeneous
+communities here
+
+
+\begin{figure}[hbt]
+
+  \centerline{ \includegraphics[width=100mm, height=75mm]{pt_hedland_communities.png}}
+
+  \caption{Location of indigeneous communities in study area.}
+  \label{fig:communities}
+\end{figure}

production/pt_hedland_2006/project.py

@@ -21 +21 @@
 offshore_name2 = 'pt_hedland_offshore_points_fairsheet'
 # coastline developed from 30m DTED
-coast_name = 'pt_hedland_coastline_points'
+coast_name = 'pt_hedland_coastline_points.xya'
 
 boundary_basename = 'SU-AU_clip'
@@ -61 +61 @@
 #for MOST
 gauge_filename = gaugedir + 'pt_hedland_gauges.xya'
-
+community_filename = gaugedir + 'CHINS_v2.csv'
+community_scenario = gaugedir + 'community_pt_hedland.csv'
 tidal_filename = tidedir + 'pt_hedland_tide.txt'
+coast_filename = datadir + coast_name
 
 # boundary source data

production/pt_hedland_2006/results.tex

@@ -12 +12 @@
 \begin{figure}[hbt]
 
-  \centerline{ \includegraphics[width=75mm, height=75mm]{pt_hedland_data_poly.png}}
+  \centerline{ \includegraphics[width=100mm, height=75mm]{pt_hedland_data_poly.png}}
 
   \caption{Study area for Pt Hedland scenario}
@@ -21 +21 @@
 \begin{figure}[hbt]
 
-  %\centerline{ \includegraphics[width=75mm, height=75mm]{}}
+  %\centerline{ \includegraphics[width=100mm, height=75mm]{}}
 
   \caption{Computational mesh for Pt Hedland study area}