Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1789 1267 (N12W54 L= 87)   326    113        0 (1743/1246-1822/1288) -14.0 W49-W57 N10-N13
      2  687 1256 (N11E21 L= 12)   219    103        0 ( 638/1240- 739/1273) -18.4 E23-E17 N10-N12
      3 1405 1150 (N04W23 L= 56)   268     75        0 (1370/1113-1437/1193) -51.3 W21-W25 N01-N06
      4  283  871 (S10E51 L=342)   956    195        0 ( 268/ 776- 302/ 975)  79.7 E53-E49 S15-S04
      5 1904  902 (S10W67 L= 99)   517    155        0 (1865/ 885-1934/ 922)   5.7 W61-W71 S10-S08
      6  590  836 (S14E28 L=  5)   363    172        1 ( 542/ 780- 662/ 863)  18.4 E31-E22 S16-S11
      7  139  809 (S13E69 L=324)  1296    235        1 ( 111/ 772- 165/ 847) -29.7 E74-E64 S15-S10
      8  594  656 (S25E30 L=  3)   389     81        1 ( 570/ 614- 615/ 690) -90.0 E31-E28 S27-S22
      9  371  564 (S30E50 L=343)  1440    308        5 ( 309/ 490- 419/ 642)  33.7 E58-E43 S34-S24
     10 1681  519 (S34W51 L= 84)   414     73        0 (1647/ 497-1702/ 547)  36.9 W48-W53 S36-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