Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  741 1825 (N41E51 L= 55)  2175    197        2 ( 689/1758- 806/1906) -18.4 E71-E38 N37-N48
      2  442 1538 (N15E53 L= 53)  1237    168        3 ( 407/1466- 471/1628) -57.5 E57-E49 N10-N21
      3  316 1559 (N14E69 L= 38)   569    101        0 ( 275/1519- 354/1588)  50.2 E71-E65 N10-N16
      4  528 1564 (N18E47 L= 59)   296     41        1 ( 507/1546- 556/1582)   0.0 E49-E44 N17-N19
      5  380 1480 (N10E55 L= 51)   343     50        0 ( 364/1457- 405/1499)  56.3 E57-E53 N08-N11
      6  317 1243 (S05E51 L= 55)   917    295        0 ( 241/1125- 368/1354)  74.1 E56-E48 S12-N02
      7  726 1087 (S05E19 L= 88)  1010    303        2 ( 616/1013- 807/1206)  63.4 E24-E15 S11-N03
      8 1195 1021 (N01W10 L=116)  1142    294        3 (1154/ 898-1241/1160) -55.3 W05-W15 S05-N08
      9 1893 1047 (N20W64 L=171)  1916    271        1 (1818/ 934-1943/1153) -21.4 W53-W76 N14-N23
     10 1702  551 (S13W61 L=167)  2942    360        6 (1538/ 478-1830/ 633)  37.3 W49-W72 S21-S05
     11 1623  363 (S26W70 L=176)  1194    103        2 (1575/ 335-1668/ 397)   2.9 W61-W81 S27-S23
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y