Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  945 1558 (N32W00 L=270)   397    113        0 ( 897/1525- 999/1579) -33.7 E03-W03 N29-N34
      2 1714 1519 (N23W58 L=327)   497     69        2 (1688/1495-1753/1544) -18.4 W54-W61 N20-N24
      3 1744 1475 (N20W59 L=328)   437     50        2 (1722/1448-1768/1495) -63.4 W56-W62 N17-N21
      4 1687 1445 (N18W51 L=321)   552    104        5 (1645/1404-1732/1491) -71.6 W48-W54 N15-N21
      5 1201 1168 (N06W12 L=281)   432    108        1 (1158/1137-1258/1199) -38.7 W09-W14 N03-N07
      6 1742 1084 (S04W48 L=317)   283     30        0 (1721/1062-1763/1108)  45.0 W46-W49 S05-S02
      7 1815  921 (S14W55 L=324)   721    102        0 (1780/ 869-1836/ 962) -50.2 W50-W56 S17-S11
      8 1781  839 (S19W52 L=321)   357     37        1 (1762/ 812-1797/ 868)  90.0 W49-W53 S20-S16
      9  140  805 (S05E70 L=199)  1503     97        2 ( 103/ 775- 186/ 842) -20.0 E77-E64 S07-S02
     10 1672  792 (S21W42 L=311)   266     35        0 (1654/ 768-1692/ 818)  33.7 W40-W43 S22-S19
     11 1484  698 (S25W27 L=297)   280     66        0 (1459/ 672-1515/ 726)  36.9 W25-W29 S26-S23
     12  304  601 (S19E58 L=212)  2666    322        9 ( 237/ 526- 392/ 667)  14.7 E67-E48 S23-S15
     13  927  495 (S33E13 L=257)   406    113        0 ( 870/ 458- 959/ 521) -33.7 E16-E10 S35-S30
     14  728  467 (S33E28 L=242)   865    212        1 ( 641/ 411- 823/ 511)  13.2 E36-E19 S35-S30
     15 1172  416 (S41W04 L=274)   620    211        2 (1084/ 375-1275/ 443)   7.6 E03-W12 S43-S39
     16 1417  329 (S52W28 L=297)  1451    212        3 (1336/ 301-1513/ 358)   3.0 W19-W38 S54-S48
     17 1198  280 (S53W06 L=275)   560     72        1 (1155/ 254-1239/ 311)  14.0 W01-W09 S55-S50
  

    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