Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1732 (N43W31 L=206)   652    107        0 (1074/1690-1169/1767)   5.2 W25-W36 N41-N44
      2 1186 1389 (N20W19 L=194)   675    225        1 (1146/1291-1235/1487)  36.9 W13-W25 N15-N24
      3  586 1366 (N34E20 L=155)   857    201        2 ( 530/1299- 671/1434)  18.4 E25-E12 N30-N36
      4  823 1383 (N29E03 L=172)   215     62        0 ( 788/1369- 859/1395) -11.3 E05-W00 N27-N29
      5  284 1294 (N36E49 L=127)   357     55        1 ( 253/1275- 309/1313)   0.0 E52-E45 N34-N37
      6 1810 1025 (S16W52 L=227)   802    134        4 (1752/ 984-1876/1054)  -4.1 W46-W60 S17-S14
      7 1660  826 (S24W36 L=211)   499    125        1 (1623/ 776-1679/ 883)  51.3 W32-W37 S26-S20
      8 1009  768 (S10E07 L=168)   535    143        3 ( 960/ 732-1079/ 804) -33.7 E10-E03 S12-S07
      9 1804  729 (S34W50 L=225)   667     87        0 (1767/ 700-1834/ 766)  21.8 W45-W54 S36-S32
     10  828  479 (S23E26 L=149)   581    214        0 ( 772/ 390- 882/ 557)  29.7 E34-E20 S27-S19
     11  605  295 (S30E54 L=121)  1430    215        2 ( 559/ 202- 662/ 370)  51.8 E63-E49 S38-S24
  

    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