Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1420 1462 (N24W37 L=199)  1417    269        3 (1307/1371-1494/1551) -71.6 W32-W39 N17-N32
      2 1617 1508 (N22W55 L=217)   461    147        0 (1554/1478-1674/1527) -29.1 W48-W59 N18-N23
      3 1014 1220 (N19W03 L=165)   210    105        0 ( 973/1186-1049/1262) -59.0 W01-W05 N16-N21
      4 1285 1180 (N11W19 L=181)   609    347        0 (1220/1115-1381/1224)   0.0 W14-W25 N06-N13
      5 1392  939 (S05W20 L=182)   327    116        1 (1343/ 920-1451/ 960) -33.7 W17-W24 S07-S03
      6 1663  752 (S24W38 L=200)  1902    523        2 (1516/ 598-1796/ 957)  30.7 W24-W51 S30-S14
      7 1634  923 (S12W37 L=199)   348    142        0 (1601/ 892-1664/ 956)  36.9 W34-W39 S13-S09
      8 1506  845 (S14W26 L=188)   331    100        0 (1481/ 803-1523/ 894)  59.0 W24-W28 S15-S10
      9  449  868 (N07E40 L=122)   282     69        0 ( 420/ 849- 481/ 884) -11.3 E42-E37 N05-N07
     10 1738  684 (S31W46 L=208)  1084    227        1 (1674/ 598-1807/ 763)  26.6 W39-W53 S34-S27
     11  347  504 (S15E63 L= 99)   476     46        1 ( 332/ 476- 363/ 531)  33.7 E65-E59 S17-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