In medical imaging, several authors have proposed to characterize roughness of observed textures by their fractal dimensions. Fractal analysis of 1D signals is mainly based on the stochastic modeling using the famous fractional Brownian motion for which the fractal dimension is determined by its so-called Hurst parameter. Lots of 2D generalizations of this toy model may be defined according to the scope. This lecture intends to present some of them. After an introduction to random fields, the first part will focus on the construction of Gaussian random fields with prescribed invariance properties such as stationarity, self-similarity, or operator scaling property. Sample paths properties such as modulus of continuity and Hausdorff dimension of graphs will be settled in the second part to understand links with fractal analysis. The third part will concern some methods of simulation and estimation for these random fields in a discrete setting. Some applications in medical imaging will be presented. Finally, the last part will be devoted to geometric constructions involving Marked Poisson Point Processes and shot noise processes. Definition 2. The distribution of (X t) t∈T is given by all its finite dimensional distribution (fdd) ie the distribution of all real random vectors (X t 1 ,. .. , X t k) for k ≥ 1,t 1 ,. .. ,t k ∈ T. Note that joint distributions for random vectors of arbitrary size k are often difficult to compute. However we can infer some statistics of order one and two by considering only couples of variables. Definition 3. The stochastic process (X t) t∈T is a second order process if E(X 2 t) < + ∞, for all t ∈ T. In this case we define • its mean function m X : t ∈ T → E(X t) ∈ R; • its covariance function K X : (t, s) ∈ T × T → Cov(X t , X s) ∈ R. We have therefore a very simple condition on variogram for Gaussian random fields with stationary increments. Corollary 2. Let X be a centered Gaussian field with stationary increments. If for all ε > 0, there exist c 1 , c 2 > 0, such that for all x ∈ [−1, 1] d , c 1 x 2γ+ε ≤ v X (x) = E((X x − X 0) 2) ≤ c 2 x 2γ−ε , then X has critical Hölder exponent on [0, 1] d equal to γ.