More stream lines in a ListStreamPlot

Question

To provide a bit of background, I have a set of wind flow data from a two-dimensional computational fluid dynamics simulation in an urban environment. The data that are returned from the simulation are two matrices, u and v, that highlight the wind flow velocity in the horizontal and vertical directions of the plane:

To briefly go over the above images, the first is a false-color plot of the u matrix. In this case, warm colors correspond to regions of high wind flow velocities in the horizontal direction (going from right to left) while cooler colors correspond to either regions with low wind flow velocities or those velocities going in the opposite direction. The light blue polygons in the image are the buildings, which have no horizontal wind flow velocity.

The second image is a false-color plot of the v matrix. In a similar fashion as before, warm colors correspond to regions of high wind flow velocities in the vertical direction (going from top to bottom) while cooler colors correspond to either regions with low wind flow velocities or those velocities going in the opposite direction. The yellow-green polygons in the image are the buildings, which have no vertical wind flow velocity.

As I have done in the past for my computational fluid dynamics plots, I combined both of these matrices into a form that could be understood by Mathemtica as a ListStreamPlot. That is, for each spatial location in my u and v matrix, I pass out an array of values: {{{x_location,y_location},{u_value,v_value}},...}. For this particular simulation, I only considered rows of u and v that were sandwiched between the buildings; a copy of the data contained within a Mathematica notebook can be obtained from here. When I went to go and make a ListStreamPlot,

ListStreamPlot[data, AspectRatio -> 0.5, StreamPoints -> Fine, PerformanceGoal -> "Quality", StreamColorFunction -> "Rainbow"]

however, I ended up with far fewer stream lines than I have been accustomed to seeing:

Consequently, my question is this: how can I ensure that Mathematica includes more streamlines in the plot? (It is apparent, from the above images of u and v, that there should be more streamlines that could be generated.) Moreover, I would like to ensure that the arrows on the ListStreamPlot are included in the final visualization along with the chosen color scheme; the latter of these desires seems to be ignored in my ListStreamPlot, which leads me to believe that there might need to be some post-processing done to the data.

Also, as a final note, I have already looked at the discussion here and the highest rated response from Jens has not worked in the context of my data.

Answer 1

Here is a simple adaptation of my answer mentioned in the question:

data = Table[{-1 - x^2 + y, 1 + x - y^2}, {x, -3, 3, .2}, {y, -3, 
    3, .2}];

t = Map[ListStreamPlot[data, DataRange -> {{-3, 3}, {-3, 3}}, 
     StreamStyle -> "Line", 
     StreamPoints -> RandomReal[{-3, 3}, {#, 2}], ImageSize -> 500] &,
    ConstantArray[50, 20]];

Show[t]

I used a simple data set data from the docs, and added the appropriate DataRange to the ListStreamPlot.

In order to try this with your own data, just do a single ListStreamPlot first and make sure the data range is correct. With the data range as above, RandomInteger would make no sense. But as I mentioned in the other answer, you can replace RandomInteger by some function of your choice. It just has to be constructed to fit the desired amount of lines into the plot range.

To include the color and arrow information, the following works:

t = Map[ListStreamPlot[data, DataRange -> {{-3, 3}, {-3, 3}}, 
     StreamColorFunction -> "Rainbow", 
     StreamPoints -> RandomReal[{-3, 3}, {#, 2}], ImageSize -> 500] &,
    ConstantArray[20, 20]];

Show[t]

Here I added the StreamColorFunction option and removed the line style for the stream lines. Also, to make the arrows recognizable, I reduced the number of stream lines (that's determined by the arguments of ConstantArray[20, 20]).

Edit: the actual data

After looking at the actual data linked in the question, it does look like there is something wrong with the way ListStreamPlot colors the lines. Strangely, this can be fixed by just omitting one element of the list data. It doesn't seem to matter which element I omit. So here I leave out the last element, and everything works. It takes quite long to make the plot because of the large data file:

t = Map[ListStreamPlot[Most[data], StreamColorFunction -> "Rainbow", 
     StreamPoints -> 
      Transpose[{RandomInteger[{1, 583}, #], 
        RandomInteger[{1, 29}, #]}], ImageSize -> 500] &, 
   ConstantArray[20, 20]];
Show[t, AspectRatio -> .5, PlotRange -> {{1, 583}, {1, 29}}]

In addition to adding the different AspectRatio, I also specified the PlotRange according to the data.

Answer 2

Well, right now I have a solution without correct colouring. Trying to get the colour map right has led me into all sorts of trouble, and there's still new problems that I might have to ask another question about. But here's what works so far.

{{x, y}, {u, v}} = Transpose[data, {3, 1, 2}];
{{x0, x1}, {y0, y1}} = Through[{Min, Max}[#]] & /@ {x, y};
aspectRatio = 1/2;
scale = {1/(x1 - x0), aspectRatio/(y1 - y0)};
Show[
 ListStreamPlot[Transpose[{scale {x, y}, scale {u, v}}, {2, 3, 1}], 
   AspectRatio -> aspectRatio, StreamPoints -> Fine, ImageSize -> 500] 
  /. Arrow[pts_] :> Arrow[#/scale & /@ pts], 
 PlotRange -> {{x0, x1}, {y0, y1}}]

Note that you can't just say StreamColorFunction -> "Rainbow" on this plot, because then it'll choose colours based on the norms of the scaled vectors we've fed into ListStreamPlot, and those are not proportional to the norms of the original vectors because the scaling is nonuniform.