I have a point feature layer with ~240,000 parcel centroids. I have a second feature layer containing ~9000 subdivision polygons. I want to assign the name of the associated subdivision as a new field to each point in the parcel centroid layer. The polygon layer is a shapefile containing a single field, the subdivision name. The point layer is a csv file containing 3 fields - the parcel ID, lat, lon and is loaded using layer->Add Delimited Text Layer.
The polygon layer is in WGS84/EPSG:4326. I have explicitly set the point layer to the same CRS. Everything looks sane on the map. I then try to do a spatial join: Vector->Data Management Tools->Join Attributes by Location. I set the target layer to the point layer, the join vector layer to the polygon layer. I set attribute summary to "Take attributes of first feature". I set the output shape file and click OK. I immediately get a warning message that the two layers have different CRS. But I have verified that are both EPSG:4326, as is the project CRS. I click OK and the progress bar immediately shows 15% and then life comes to a halt. The task monitor shows qgis running at 100% of CPU, but after several hours, it's still running at 100%, and the progress bar has not updated. I spatially indexed both layers.
Why did I get the different CRS message?
Why is this taking so long to run? (qgis 2.0.1 Dufour, Windows 7/64bit, Intel i7 2.4GHz)
Is there a better/different approach I should try?
Edit: For what it's worth, ArcMap 10.2 did the spatial join in about 15 minutes (rough estimate - I wasn't actually timing it.)
If anyone wants to try this join, the data can be found here