Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1303 1777 (N51E08 L=213)   374     47        0 (1264/1761-1330/1790)  31.0 E11-E05 N48-N52
      2  750 1696 (N26E39 L=182)   362    105        0 ( 726/1659- 779/1745) -56.3 E41-E36 N23-N29
      3  521 1640 (N16E52 L=169)  1070    173        0 ( 501/1566- 537/1719) -45.0 E60-E47 N11-N21
      4 1515 1676 (N51W15 L=236)   864    117        3 (1461/1638-1568/1714) -10.3 W08-W20 N49-N53
      5  158 1191 (S17E68 L=153)  1713    167        2 ( 130/1114- 193/1261)  77.5 E73-E64 S21-S12
      6  318 1049 (S22E47 L=175)  1852    362        4 ( 150/ 934- 441/1124)  25.8 E66-E36 S31-S16
      7 1686  999 (N12W40 L=262)   750    111        0 (1639/ 951-1725/1035)  -7.1 W36-W44 N10-N13
      8 1698  914 (N08W44 L=265)   271     51        0 (1675/ 892-1724/ 941) -11.3 W41-W46 N07-N08
      9 1400  891 (S03W24 L=246)   504    116        0 (1364/ 845-1445/ 936) -15.9 W21-W28 S03-N00
     10  766  751 (S28E08 L=214)  1116    277        0 ( 723/ 670- 829/ 838)  84.8 E10-E01 S32-S21
     11 1541  602 (S13W44 L=266)   492    159        1 (1468/ 577-1621/ 627)  31.0 W39-W49 S15-S09
     12 1614  523 (S15W55 L=276)  1170    230        3 (1540/ 454-1664/ 590)  84.3 W49-W59 S20-S10
     13 1506  491 (S20W48 L=270)   334     92        0 (1480/ 455-1531/ 526) -78.7 W46-W49 S22-S17
     14  800  420 (S48W06 L=228)  1660    389        2 ( 633/ 352- 963/ 505)   3.7 E10-W21 S50-S44
     15 1130  220 (S49W47 L=268)   371     53        1 (1103/ 208-1167/ 240)  33.7 W42-W50 S51-S47
     16 1230  177 (S48W61 L=282)  1138     91        4 (1192/ 137-1265/ 220)  -5.4 W51-W73 S49-S45
  

    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