Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1181 1736 (N45W10 L= 98)   280     34        1 (1160/1720-1205/1756) -26.6 W07-W11 N43-N46
      2 1357 1715 (N44W25 L=113)   294     25        0 (1337/1693-1373/1737)  45.0 W22-W26 N42-N46
      3  621 1613 (N33E33 L= 56)   292     82        1 ( 597/1571- 651/1651) -40.6 E36-E29 N29-N35
      4 1397 1494 (N27W24 L=112)   393    148        1 (1346/1450-1447/1557) -54.5 W20-W26 N24-N31
      5 1812 1415 (N25W62 L=150)   948    107        2 (1777/1364-1850/1455) -39.8 W58-W65 N22-N28
      6  258 1285 (N11E56 L= 32)   597     97        0 ( 244/1241- 275/1336) -71.6 E58-E53 N07-N13
      7 1780 1222 (N12W52 L=141)   570     86        2 (1743/1188-1813/1253) -15.9 W48-W55 N10-N14
      8 1732 1143 (N07W47 L=136)   532    119        2 (1708/1083-1754/1193)  81.9 W45-W49 N03-N10
      9 1864 1050 (N03W61 L=149)  2090    220        3 (1825/ 943-1912/1155)  76.0 W57-W67 S03-N09
     10 1804 1092 (N05W54 L=142)   319     57        1 (1787/1072-1828/1119)  26.6 W52-W56 N03-N06
     11  809  876 (S14E13 L= 76)   398    166        1 ( 762/ 814- 872/ 915)  40.6 E15-E08 S17-S11
     12  142  837 (S16E71 L= 18)   866     90        1 ( 102/ 803- 180/ 869)   9.0 E83-E64 S17-S14
     13 1014  776 (S19W00 L= 89)   419     78        1 ( 991/ 740-1042/ 821)  90.0 E01-W02 S21-S16
     14  699  750 (S22E20 L= 68)   336    161        0 ( 619/ 729- 768/ 777)  -5.7 E25-E15 S22-S20
     15  964  733 (S22E03 L= 86)   238     99        0 ( 907/ 712-1001/ 748)  -9.5 E06-W00 S23-S21
     16  851  661 (S27E10 L= 78)   237     68        0 ( 827/ 630- 882/ 695) -53.1 E11-E08 S28-S24
     17 1160  535 (S34W12 L=100)   380     73        0 (1124/ 519-1201/ 557)   0.0 W09-W15 S35-S32
  

    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