Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1777 1406 (N18W66 L=159)  1037    109        3 (1739/1364-1807/1454)  17.1 W60-W74 N15-N20
      2  564 1305 (N27E29 L= 64)   390    120        1 ( 520/1255- 603/1326) -39.8 E33-E26 N23-N28
      3 1439 1296 (N17W31 L=124)   261     91        1 (1396/1271-1485/1320) -38.7 W28-W33 N15-N19
      4  394 1253 (N25E44 L= 49)   356    131        1 ( 366/1210- 448/1278)  23.2 E46-E38 N22-N26
      5 1249 1162 (N11W15 L=109)   224     97        1 (1203/1147-1290/1181)   0.0 W12-W18 N10-N12
      6  167 1055 (N13E68 L= 25)  1203    106        2 ( 132/ 997- 189/1109)  74.1 E74-E64 N09-N16
      7  990  830 (S06E04 L= 89)  1090    217        3 ( 953/ 711-1024/ 931)  85.6 E05-E02 S13-N00
      8 1752  854 (S15W50 L=143)   538     46        0 (1727/ 830-1777/ 876) -31.0 W47-W52 S16-S12
      9  548  637 (S14E36 L= 57)   256     62        0 ( 521/ 619- 582/ 651)  11.3 E38-E33 S15-S13
  

    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