Hot answers tagged

40

The following solution is based on a post by Roger Bivand on R-sig-Geo. I took his example replacing the German shapefile with some census data from Oregon you can download from here (take all shapefile components from 'Oregon counties and census data'). Let's start with loading the required packages and importing the shapefile into R. # Required packages ...


18

The first step to determining the correct projection of any layer, is to find the projection information, if any, that came with your layer. In the case of a Shapefile, like what you downloaded from Census.gov, that information is contained in a .prj file, short for Projection. Here are the contents of the projection file from the census data: GEOGCS["...


15

Here is a solution using the sf package: library(tidycensus) library(dplyr) library(sf) library(ggplot2) # get data from tindycensus for demonstration (note you need an API key, folow instructions here: https://walkerke.github.io/tidycensus/articles/basic-usage.html) census <- tidycensus::get_acs(geography = "tract", variables = "B19013_001", ...


13

These are Federal Information Processing Standard (FIPS) state codes, and the ones you mention are missing because, from the wiki: "certain numeric codes "are reserved for possible future use in identifying American Samoa (03), Canal Zone (07), Guam (14), Puerto Rico (43), and Virgin Islands (52)", but these codes were omitted from FIPS PUB 5-2 without ...


11

GEOID is the field used to join TIGER/Line geographic data to the demographic data in various American Community Survey products and in the Decennial Census. It is slightly confused by the fact that this field is called GEOID10 only in the TIGER/Line 2010 products (in fact, almost all of the field names in TIGER/Line 2010 end in 10), and by the fact that the ...


9

Basically, you can do the extract using ogr2ogr as long as you give the Census tract ID, so it's really an issue of getting 72,000 ogr2ogr calls. ogr2ogr -where "tract = '<tract_id>'" /dest_folder /source_folder block_shapefile -nln block_shapefile_<tract_id> Notes: You don't have to specify the source format, ogr2ogr will figure it out. You ...


8

Your output field should be a text field. If you set the field calculator's parser to Python, you could enter this for the formula: '{}.{}'.format( !TractField! , !BlockField! ) If they really are double fields storing integer IDs for tracts and blocks, then enter: '{}.{}'.format( int(!TractField!) , int(!BlockField!) ) If they are double fields storing ...


7

They are available on the Census FTP site at http://www2.census.gov/geo/tiger/TIGER2010BLKPOPHU/


7

I recently dealt with the exact same issue with Statistics Canada census profile CSV files containing 2172 columns. You can import your csv into an ESRI File Geodatabase (FGDB) if you have access to ArcGIS. According to ESRI, the FGDB format can handle 65,534 fields in a feature class or table. In my case, I was able to import my 2172 column wide CSV file ...


7

PostgreSQL has a column limit of between 250 and 1600 "depending on column types", and supports spatial data and queries with the PostGIS extension. So I would be inclined to do two things: First, where a column represents a category rather than free text, create a separate table with those categories, and replace the column with an integer ID and foreign ...


7

Very interesting question! I'd note that the geography area calculation is extremely close (in your last example, 3 msq difference over an area of 20000000 msq. (0.000015%) Since we know that census manages its data inside Oracle Spatial, I'd guess that what you're seeing is a very small difference in the implementation of geodetic area calculation between ...


7

What you're trying to do is known as apportionment. This takes a numerical attribute value of a feature and divides it in some way between pieces that feature is split into. There are a number of different solutions. In fact, there's a specific tool for it in the Business Analyst extension, but you probably don't have that. Esri has also published a ...


6

I don't think I can directly link, but what you want is American Community Survey Table S1810 with New Mexico County Geography. Start here: http://factfinder2.census.gov/faces/nav/jsf/pages/searchresults.xhtml?refresh=t Enter S1810 for your search topic. This should give you four tables named Disability Characteristics. From the options on the left, now ...


6

While some districts may follow block group, tract or even county boundaries, many of them (probably most -- might depend on the individual state) were based on block-level data. Since the districts have to be drawn based on Census data, and since the smallest level of geography the Census Bureau uses is the block, there should be no districts that violate ...


6

The Short Answer: No. Census Tracts do change over time (Especially from 2000 to 2010 when in many cases they appear to be essentially redrawn). You can use the faces files to compare the same geographic area with its 2000 Census and 2010 Census assignment using dissolve techniques. You can also check out the relfiles to compare 1990 Census with 2000 Census ...


6

American Fact Finder data (tabular) needs to be joined to TIGER data. The TIGER data delineates the boundaries you are working with, be it State, County, MCD, Census tracts, block groups, blocks, etc. Once you have the data joined through a common field you can start examining the relationships between the data.


6

These data are available from the 2010 Census or the 2014-2018 American Community Survey (ACS) 5-year estimates. (1 year data are only for places with 60,000 population or more, so you need to use the 5-year data or go back to 2010 Census) Probably the easiest way is to check out NHGIS to get table B01003 for all places in the U.S. in a shapefile (point) ...


5

It turns out that I needed to add a WHERE ST_Intersects to my ST_Intersection Query, as follows: SELECT sum(((st_area (st_intersection (p.the_geom,c.the_geom))/st_area(c.the_geom))*ci.pop2000)) AS Parcels_pop FROM parcel_proj p, census_proj c, tgr39035sf1blk ci WHERE ST_Intersects(p.the_geom,c.the_geom) and ci.stfid=c.stfid; This may not matter for others ...


5

Downvote me if I'm wrong, but the GEOID is a concatenation of a bunch of fields as noted in the summary file documentation on page 13. In the geographic area codes, you'll be looking for fields at positions 26-65.


5

The U.S. Census Bureau's TIGER Products have free shapefiles for the American Community Survey 5-Year Estimates and other datasets. These include several measures of income and poverty, including "Poverty Status In The Past 12 Months By Household Type By Age Of Householder", "Poverty Status Of Individuals In The Past 12 Months By Living Arrangement", and "...


5

You might find the TIGER page helpful in finding information. For instance, this shapefile purports to contain basic demographic information by 2010 Census Tract. There are other types of information you might find interesting as well. If you want to find your own Census Tracts, you can use this user interface, known as the American FactFinder to find your ...


5

Try NHGIS, which allows you to get blocks for all states. Set the filters as I have done in the attached screen shot.


5

The Census Bureau has prepared so-called "Relationship Files" that describe the geographic relationships between geographies (including tracts) at different times. This includes a lot of information for going between the two datasets by population or land area (a complete list is below). In your case, since you're trying to go from 1990 to 2010, you will ...


5

From the US Census http://www.census.gov/popclock/ it appears to be updated by the second.


5

Save both files as the same CRS. Load each of the files and then run points in polygon with the newly loaded files. The project on the fly is a visual aid. When running analysis it will use the CRS of the file.


5

I have written a set of scripts to solve this problem and have published the output in case anyone else needs it. See the bottom of this answer for details. GeospatialPython.com pointed me in the right direction. I was able to throw together a quick Python script to do this. I'm sure there are issues here because I'm not very familiar with pyqgis, but it ...


4

I would check that you have downloaded the correct census files. I've spent many wasted days trying to figure out why I could not make an obvious join, only to discover that 2010 census files per state and per county sometimes have different codes. (I know this sounds strange but it is true). Secondly, in the pre-2010 census, you should check how the ...


4

Getting the Data You can acquire US Census Place Data, which includes cities, here. For instance, Rhode Island's 2017 Place Data is available here. For ease of bulk downloading, the full list 2017 place data is: https://www2.census.gov/geo/tiger/TIGER2017/PLACE/tl_2017_01_place.zip https://www2.census.gov/geo/tiger/TIGER2017/PLACE/tl_2017_02_place.zip ...


4

check this link http://www.geohive.com/cntry/ecuador.aspx for the census data (1990,2001 and 2010) this one http://www.fao.org/geonetwork/srv/en/metadata.show?id=1174&currTab=simple for the admin limits By the way it is better to add the info that you are looking for the admin regions called "Parishes". These are third-level administrative units of ...


Only top voted, non community-wiki answers of a minimum length are eligible