An extended comment...
The word "correlated" is many times used to suggest some sort of statistical dependence of one variable on another (especially outside the field of Statistics). But there are many forms of "statistical dependence" so no summary statistic like a correlation coefficient always gives an adequate picture of what's going on. Hence, one looks at a variety of forms of dependence. Fortunately you seem to have lots of data which gives you a lot of leeway (but also the potential for more complication).
One aspect of "statistical dependence" is that if there is some sort of dependence, then the distribution of values for one variable will depend on the values of another variable. So with lots of data we could look at a few histograms where the range of one of the variables is restricted. Here are two such slices:
(* Generate some data *)
data = Transpose[{RandomVariate[LogNormalDistribution[3.5, 1], 50000],
RandomVariate[LogNormalDistribution[3.5, 1], 50000]}];
(* Slice 1 *)
data1 = Select[data, 15 < #[[1]] < 20 &][[All, 2]];
Histogram[data1, {10}, "PDF", PlotRange -> {{0, 250}, {0, 0.02}}]

(* Slice 2 *)
data2 = Select[data, 100 < #[[1]] < 105 &][[All, 2]];
Histogram[data2, {10}, "PDF", PlotRange -> {{0, 250}, {0, 0.02}}]

The first slice is centered around the horizontal variable being around 17.5 and the second slice is centered around 102.5. I don't see very much difference between the two histograms (displayed as PDF's). (Of course, the two variables generated were by definition independent so the histograms ought to look very similar.)
An alternative is to use some sort of Generalized Additive Model (gam) (which is a nonparametric regression approach). While some folks have created their own gam packages in Mathematica, there is no built-in Mathematica function to do so (yet?). A very, very crude approximation is to take slices across one of the variables and calculate the means of the other variable for each slice. Then plot the means with the data in the background.
bw = 10 (* bandwidth *)
means = Table[{x,
Mean[Select[data, x - bw/2 < #[[1]] < x + bw/2 &][[All, 2]]]},
{x, bw/2, 250, bw}];
Show[ListPlot[data],
ListPlot[means, Joined -> True, PlotStyle -> Red]]

Here, too, there does not seem to be evidence of the vertical variable statistically dependent on the horizontal variable.
SmoothDensityHistogram[]
? $\endgroup$x1 = RandomVariate[LogNormalDistribution[3.5, 1], 50000]; x2 = RandomVariate[LogNormalDistribution[3.5, 1], 50000]; ListPlot[Transpose[{x1, x2}], PlotRange -> {{0, 250}, {0, 250}}, AspectRatio -> 1, AxesStyle -> White, TicksStyle -> Black]
. $\endgroup$x1
and between 100 and 105 forx1
and compare the resulting histograms forx2
. If they differ much in any observed way, then further investigation of the correlation structure is warranted. $\endgroup$SmoothDensityHistogram[]
function is what I was looking for. $\endgroup$