Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  707 1838 (N62E25 L=358)   589     49        1 ( 676/1810- 733/1862)  26.6 E29-E21 N58-N64
      2 1340 1519 (N26W27 L= 49)  1008    212        2 (1226/1475-1427/1545) -26.6 W19-W31 N22-N28
      3  668 1378 (N24E20 L=  2)   967    248        1 ( 560/1313- 778/1424) -28.1 E27-E12 N18-N27
      4  339 1407 (N29E48 L=335)   396     75        0 ( 302/1388- 381/1421)   7.1 E51-E43 N28-N30
      5 1230 1405 (N20W17 L= 39)   201    122        0 (1178/1378-1289/1419) -23.2 W13-W20 N17-N20
      6 1701 1383 (N14W51 L= 73)   579    131        3 (1643/1361-1755/1407) -11.3 W46-W56 N12-N15
      7 1244 1173 (N06W15 L= 37)   591    210        0 (1179/1130-1345/1229) -29.1 W11-W20 N02-N09
      8  987 1117 (N05E01 L= 21)   245     89        1 ( 977/1069-1002/1165)  80.5 E02-W00 N01-N07
      9 1829  939 (S13W57 L= 79)   449     62        1 (1795/ 919-1858/ 961)  18.4 W53-W60 S14-S11
     10 1725  818 (S20W47 L= 69)   387    116        2 (1695/ 780-1756/ 844) -45.0 W44-W49 S22-S17
     11 1645  734 (S24W41 L= 63)   874    236        1 (1622/ 632-1679/ 842)  90.0 W38-W43 S31-S17
     12  804  578 (S26E19 L=  3)   983    280        2 ( 707/ 471- 873/ 682)  38.7 E28-E13 S32-S20
  

    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