I have a group of geographical polygons in lat long (epsg:4326)
coordinates.
I would like to lay this polygon out at an angle.
I would like to put two 3D bar charts at the midpoint or close inside
each polygon.
Can anyone provide an example of this? Or a web site with an example?
Thanks,
Jim
28480 Levels: 481130001001000 481130001001001 481130001001002 ...
481130199004013
Thanks,
Jim Burke
PS My thanks to Gledson Luiz Picharski for his help initially with this
"outside of %IN%" logic.
___
R-sig-Geo mailing list
R-sig-Geo@s
",
package="maptools")[1], proj4string=CRS("+init=epsg:2166"))
Error in getinfo.shape(filen) : Error opening SHP file
Good luck,
Jim Burke
Jim Burke wrote:
Jian,
First realize where your R code was trying to
find your shapefile. In your code below you
are searching
<- readShapePoly(system.file("nro2.shp", package="maptools")[1],
proj4string=CRS("+init=epsg:2166"))
Then in R, make sure R is pointing to the
directory where your nro2.shp file exists.
This can be done by selecting in R, File,
Change Directory.
Good luck,
J
Hi All,
Everyone's suggestions work. There are lots
of roads to solutions in R. This is wonderful!
It seems that Dan's and Pedro's solutions
work faster for my largish 1,897 spatial
polygon data frame.
Thanks,
Jim Burke
All re
t %in%?
I am processing two different sets of polygons (A,B) in
a loop. Operation is seeing in center points in polys B
fall in A. I would like to remove the spatial polygons
that I found (B) as I traverse the geographical area.
Thanks,
Jim Burke
___
R-sig-Geo m
GridDataFrame", CRSobj="CRS"
x="SpatialLines", CRSobj="CRS"
x="SpatialLinesDataFrame", CRSobj="CRS"
x="SpatialPixelsDataFrame", CRSobj="CRS"
x="SpatialPoints", CRSobj="CRS"
x="SpatialPointsDataFrame", CRSobj=&
a KML file just fine for example.
census_sp <- readShapePoly("HD_102_census_blocks_sp.shp",
IDvar="BLKIDFP00", proj4string=CRS("+init=epsg:4326"))
Thanks,
Jim Burke
Jim Burke wrote:
How can I convert a "SpatialPolygons" object to a "polylist"
orks according to the doc with polygon data sets.
An example Map2Poly which works with "inout" lists itself as a class
"polylist".
Thanks,
Jim Burke
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
d. c_df output data frame example.
BLK PCT
"480019510006079",1811
"480019510006080",1811
"480019510006998",1811
"480019510006999",1811
"480039501003000",1120
"480039501003001",1120
"480039501003002&
"480039501003000",1120
"480039501003001",1120
"480039501003002",1120
"480039501003003",1120
Thanks,
Jim Burke
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
ER["Central_Meridian",-98.5],
PARAMETER["Standard_Parallel_1",32.1333333333],
PARAMETER["Standard_Parallel_2",33.97],
PARAMETER["Latitude_Of_Origin",31.67],
UNIT["Foot_US",0.3048006096012192]]
Thanks,
Jim Burke
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_sp$PCT, cex=.02),
as.table = TRUE)
)
plot(vote_percent, position=c(0, 0.50, 1, 1 ))
plot(total_votes, position=c(0, 0,1, .50), newpage = FALSE)
dev.off()
Good luck,
Jim Burke
Jim Burke wrote:
I would like to print two polygon maps per
PDF page. The following code prints each
map on two PD
;),
colorkey=list(height=1.0, width=0.6),
main = list(label="My second title",cex=.5),
line=0,
sp.layout = list("sp.text",coordinates(tx3_sp),
as.character(tx3_sp$PCT, cex=.02),
as.table = TRUE),
newpage = FALSE
)
dev.off()
Thanks in advance for your sugg
6286 8954
1931-40 347566 16595 160037 1750 2483
1941-50 621147 37028 354804 7367 14693
1951-60 1325727 153249 996944 14092 25467
1961-70 1123492 427642 1716374 28954 25215
1971-80 800368 1588178 1982735 80779 41254
1981-90 761550 2738157 3615225 176893 46237
1991-00 1359737 2795672 4486806 354939 9
A partial answer. plot() these separately shows
they plot at vastly different dimensions in order
to fill the plot frame.
I will try a merge in R and see what happens.
Thanks,
Jim
Jim Burke wrote:
Thanks in advance everyone. I need assistance
resolving the above error. Trying to spplot join
ight=1.0, width=0.6),
main = list(label="HD-108 - Barack Obama \n Votes by Precinct
11/04/08",cex=1),
sub = list(label="by Jim Burke, j.bu...@earthlink.net",cex=.3),
line=0,
sp.layout = list("sp.text",coordinates(tx3_sp),
as.character(tx
+datum=WGS84"),
verbose=TRUE,repair=FALSE)
Shapefile type: PolyLine, (3), # of Shapes: 29161
> class(tx3_sp)
[1] "SpatialPolygonsDataFrame"
attr(,"package")
[1] "sp"
> class(main_roads)
[1] "SpatialLinesDataFrame"
attr(,"package")
[1] "sp&q
://mobilab.mit.edu/mashup/data_format.html
MAS GIS
http://www.mass.gov/mgis/
Good luck,
Jim Burke
Alina Sheyman wrote:
> Does anyone know where I can find a shapefile with zipcodes for Boston
>
> [[alternative HTML version deleted]]
>
> ___
>
opened a bug with Thunderbird.
Back to R, I really hope to get going on the
problem in that note. I have spent some
thought and time trying to get it fixed but
no success.
Thanks,
Jim Burke
Jim Burke wrote:
Hi Everyone,
QUESTION: How can I plot a spatial polygon (SpatialPolygonsDataFrame)
as 3-d
onsDataFrame
Coordinates:
min max
r1 2488376 2542115
r2 7003736 7046215
Is projected: TRUE
proj4string :
[+proj=aea +ellps=GRS80 +datum=WGS84 +lat_1=29.5 +lat_2=45.5
+towgs84=0,0,0]
Data attributes: (I did not list these)
III.
plot(zip2_sp)
invisible(text(getSpPPolygonsLabptSlots(zip2_sp),
labels=as.character(zip2_sp$CUST_CNT), cex=0.4))
title("Wyoming Customer count", sub = "Something at the bottom of the
graph",
cex.main = 2, font.main= 4, col.main= "blue",
cex.sub = 0.75, font.su
he polygon $z
value? Assuming lat and long are in the sp.
http://r-spatial.sourceforge.net/gallery/#fig21.R
Then if its not as straightforward as above, what would be the steps to
convert a spatial polygon to points or something to graph that way?
Then also I am looking forward to a
good read of "Lattice Multivariate Visualization".
Thanks,
Jim Burke
Roger Bivand wrote:
On Mon, 16 Feb 2009, Jim Burke wrote:
(Just stick to plain text, but don't use fonts or bold or stuff - keep
it simple).
I have a nice solid SpatialPol
s. It does not behave like a Shapfile.
>main_roads <-readOGR("main roads.shp", layer = "main roads")
OGR data source with driver: ESRI Shapefile
Source: "main roads.shp", layer: "main roads"
with 29161 rows and 33 columns
Feature type: wkbL
ata source with driver: ESRI Shapefile
Source: "main roads.shp", layer: "main roads"
with 29161 rows and 33 columns
Feature type: wkbLineString with 2 dimensions
Thanks,
Jim Burke
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
#shows representation of
geometries too.
names(tx3_sp) #nice but lengthy
Let us know if this helps you any.
Good luck,
Jim Burke
Peter S. Hayes wrote:
Hi All,
I'm trying to teach a spatial analysis class using R as a means of
learning some details..
]
+ }
+ }
> save(wrld_simpl,'SIMPLEMAP.Rdata')
Error in save(wrld_simpl, "SIMPLEMAP.Rdata") :
object ‘SIMPLEMAP.Rdata’ not found
Thanks,
Jim Burke
Jim Burke wrote:
Thanks Hans,
Its been a day or so since I have used R graphics.
Now sp seems to lack grid. I cannot run
e 'application/zip' length 845890 bytes (826 Kb)
opened URL
downloaded 826 Kb
package 'sp' successfully unpacked and MD5 sums checked
The downloaded packages are in
C:\Documents and Settings\Jim Burke\Local
Settings\Temp\RtmpqCMAZS\downloaded_packages
updating HTML pac
29 matches
Mail list logo
