Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1022 1507 (N32E06 L=287)   418    130        3 (1006/1449-1043/1565) -76.0 E07-E04 N27-N35
      2 1268 1210 (N15W13 L=307)   272    102        2 (1244/1163-1291/1262) -63.4 W11-W15 N12-N18
      3 1628 1125 (N13W40 L=333)   817    272        2 (1566/1002-1688/1228) -55.3 W34-W44 N06-N19
      4 1847 1074 (N12W62 L=355)  1783    262        3 (1827/ 948-1864/1195)  78.7 W59-W64 N04-N19
      5  132 1015 (S09E73 L=220)   763     87        0 (  97/ 990- 167/1041)  -6.3 E84-E66 S10-S08
      6  757  868 (S11E15 L=278)   648    264        2 ( 701/ 746- 811/ 965)  71.6 E18-E12 S19-S04
      7  239  836 (S19E57 L=236)  3057    366       14 ( 152/ 755- 342/ 920)  16.5 E73-E46 S25-S13
      8  274  743 (S25E55 L=238)  1230    152        4 ( 246/ 657- 300/ 810) -84.3 E58-E51 S30-S20
      9 1525  584 (S21W41 L=334)   318    118        0 (1480/ 553-1553/ 629)  45.0 W37-W42 S23-S18
     10 1231  440 (S34W23 L=316)   495    163        0 (1214/ 370-1255/ 525) -70.0 W21-W26 S39-S28
  

    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