Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  922 1593 (N41E08 L=353)   279     40        0 ( 893/1575- 947/1608)   0.0 E10-E06 N39-N42
      2 1484 1376 (N25W33 L= 34)   350    125        0 (1451/1323-1516/1423) -66.8 W31-W34 N21-N28
      3  576 1282 (N19E31 L=330)   944    309        0 ( 512/1172- 655/1405) -47.0 E37-E23 N12-N27
      4 1191 1171 (N13W11 L= 11)   710    202        0 (1117/1106-1260/1242) -48.4 W06-W14 N08-N17
      5  971 1103 (N08E03 L=357)   295    163        0 ( 951/1031- 987/1175) -82.9 E04-E02 N04-N12
      6 1685 1041 (N03W45 L= 46)   366     97        1 (1664/ 995-1701/1088)  80.5 W43-W46 N00-N06
      7 1212  870 (S06W12 L= 12)   205     65        0 (1197/ 826-1221/ 904) -90.0 W10-W12 S08-S03
      8 1067  793 (S11W03 L=  3)   641    196        5 (1003/ 746-1145/ 855)  35.5 E01-W07 S13-S06
      9  954  748 (S14E05 L=356)   262    128        0 ( 931/ 688- 973/ 794)  74.1 E06-E03 S17-S10
  

    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