Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1395 1708 (N53W24 L=217)   369     84        3 (1347/1688-1428/1723)   0.0 W19-W28 N51-N54
      2  473 1590 (N29E51 L=142)   639    161        1 ( 444/1519- 498/1650) -71.6 E52-E48 N24-N33
      3  755 1388 (N20E22 L=170)   579    229        1 ( 725/1285- 784/1492) -77.9 E24-E18 N13-N27
      4 1646 1408 (N33W42 L=234)   311     70        0 (1625/1377-1670/1434) -45.0 W40-W43 N30-N34
      5 1721 1381 (N32W48 L=241)   590     91        3 (1692/1347-1749/1412) -31.0 W46-W51 N29-N33
      6 1566 1053 (N09W35 L=227)   235     51        0 (1544/1027-1592/1080) -26.6 W32-W36 N07-N10
      7 1816 1037 (N10W57 L=249)   598    111        2 (1787/ 991-1837/1074) -59.0 W54-W59 N07-N12
      8  279  847 (S19E51 L=141)  2910    787       15 ( 173/ 710- 394/ 965) -33.0 E63-E41 S26-S12
      9  422  785 (S21E38 L=154)   297     81        1 ( 404/ 743- 449/ 826) -59.0 E39-E36 S23-S18
     10 1765  677 (S12W60 L=252)  4303    590       12 (1676/ 568-1839/ 809)  82.9 W53-W67 S19-S03
     11 1840  685 (S11W69 L=262)   645    128        3 (1813/ 634-1860/ 727)  71.6 W66-W72 S14-S08
     12 1140  674 (S19W12 L=204)   435    154        0 (1058/ 655-1211/ 700)  18.4 W06-W15 S20-S16
     13 1638  673 (S14W47 L=239)   383    100        2 (1600/ 659-1673/ 692)   9.5 W43-W50 S14-S12
     14 1451  572 (S22W35 L=228)   539    111        1 (1398/ 537-1502/ 616)  50.2 W32-W38 S24-S18
     15 1714  522 (S23W60 L=253)   513    122        3 (1688/ 497-1745/ 542)   0.0 W57-W64 S24-S21
     16  811  483 (S37E08 L=184)   428     82        0 ( 782/ 457- 847/ 507)  38.7 E10-E05 S39-S34
  

    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