Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1012 1717 (N45E16 L=144)   353     60        1 ( 982/1702-1044/1733)   0.0 E18-E13 N43-N47
      2 1176 1634 (N42E00 L=160)   261     24        0 (1153/1619-1203/1649)  14.0 E02-W02 N40-N43
      3  374 1477 (N17E55 L=105)  2094    372        6 ( 306/1351- 421/1614) -80.5 E59-E51 N08-N26
      4  574 1301 (N10E33 L=127)   585    143        2 ( 537/1255- 619/1347)  -9.5 E35-E29 N06-N12
      5 1114 1246 (N14W02 L=162)   274     66        0 (1074/1221-1142/1268) -26.6 W00-W04 N13-N15
      6  344 1061 (S08E46 L=114)   272     73        0 ( 313/1035- 375/1086)  38.7 E48-E43 S10-S06
      7  609  935 (S12E24 L=135)   925    167        2 ( 571/ 846- 659/1020) -52.1 E27-E20 S16-S06
      8 1301  947 (S01W18 L=178)   329     95        1 (1254/ 917-1349/ 974)  31.0 W15-W20 S03-N01
      9  243  629 (S38E57 L=103)   588    144        0 ( 221/ 602- 260/ 664)  39.8 E60-E54 S40-S35
     10 1313  405 (S35W33 L=193)   461    115        1 (1278/ 384-1372/ 435)  -6.3 W29-W38 S35-S32
  

    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