Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1648 1423 (N38W32 L=167)  1655    254        2 (1550/1334-1732/1483)  -9.0 W22-W41 N34-N41
      2 1459 1375 (N28W17 L=151)   856    133        1 (1407/1321-1505/1422) -29.7 W13-W21 N25-N30
      3 1320 1357 (N23W08 L=143)   356     99        0 (1291/1310-1341/1408)  90.0 W07-W08 N19-N26
      4  372 1346 (S02E51 L= 84)  1336    156        3 ( 323/1279- 423/1394) -22.6 E57-E45 S05-N00
      5 1198 1262 (N13W03 L=138)   341    114        0 (1139/1242-1252/1286)  45.0 W00-W06 N10-N16
      6 1804 1190 (N28W50 L=184)   473    116        0 (1761/1135-1833/1230) -20.6 W44-W53 N25-N30
      7 1732 1094 (N20W44 L=178)   590     79        1 (1703/1049-1755/1131) -53.1 W41-W46 N17-N22
      8 1498  905 (N02W31 L=165)   746    219        2 (1472/ 813-1520/1023) -71.6 W28-W33 S03-N09
      9 1852  947 (N16W60 L=194)   415     57        1 (1830/ 925-1872/ 966)   0.0 W56-W62 N13-N16
     10 1238  903 (S06W15 L=150)  1260    300        1 (1120/ 833-1312/ 956)  13.0 W08-W21 S10-S02
     11 1653  910 (N07W41 L=176)   342     80        1 (1612/ 884-1693/ 923)  18.4 W38-W44 N04-N08
     12  539  783 (S31E25 L=110)   241     25        1 ( 517/ 762- 561/ 799)  14.0 E26-E22 S32-S30
     13 1707  530 (S11W64 L=199)  1330    183        1 (1647/ 480-1751/ 610)  84.3 W61-W69 S15-S05
     14 1582  525 (S16W53 L=187)   454     85        0 (1543/ 509-1620/ 543)  26.6 W49-W55 S17-S13
     15 1291  392 (S33W41 L=175)   271     61        1 (1259/ 381-1327/ 404)  36.9 W37-W42 S33-S30
     16 1478  331 (S30W61 L=195)  1832    175        2 (1417/ 286-1555/ 376)   2.4 W51-W75 S31-S27
     17 1157  176 (S52W57 L=191)   847    111        2 (1134/ 131-1186/ 227) -11.8 W46-W71 S54-S48
  

    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