Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1026 1797 (N52W33 L= 18)   328     49        1 (1000/1775-1045/1816)  78.7 W30-W35 N49-N54
      2  769 1701 (N54W05 L=350)  2320    424        5 ( 597/1603- 945/1794)   1.5 E15-W24 N50-N56
      3 1476 1729 (N31W59 L= 44)   680    117        0 (1425/1689-1522/1760) -81.9 W57-W61 N27-N34
      4 1412 1613 (N27W44 L= 29)   293     82        1 (1393/1564-1433/1648)  56.3 W42-W46 N23-N29
      5  710 1468 (N39E08 L=337)   373     62        2 ( 673/1443- 741/1489) -36.9 E10-E05 N37-N40
      6 1503 1442 (N15W41 L= 26)   573    129        0 (1443/1416-1563/1473) -49.4 W38-W44 N11-N18
      7  495 1131 (N23E30 L=315)  1130    304        1 ( 359/1078- 624/1185) -27.8 E40-E21 N18-N29
      8  228 1059 (N25E56 L=289)   771    139        1 ( 209/ 986- 258/1122)  69.4 E58-E51 N20-N28
      9  346 1020 (N20E44 L=301)   550    123        0 ( 300/ 976- 388/1075) -53.1 E47-E40 N16-N24
     10 1533  973 (S11W28 L= 13)   386     91        0 (1501/ 935-1560/1020)  26.6 W25-W31 S13-S09
     11 1234  856 (S10W07 L=352)   866    280        0 (1205/ 743-1289/ 991)  75.1 W05-W11 S17-S02
     12  416  916 (N13E40 L=305)   373    101        1 ( 368/ 890- 462/ 937) -33.7 E43-E37 N10-N14
     13 1655  814 (S24W34 L= 19)  1332    266        2 (1542/ 739-1724/ 908)  18.4 W24-W41 S26-S19
     14  314  807 (N09E52 L=293)   494    124        0 ( 294/ 753- 334/ 871) -90.0 E52-E50 N05-N12
     15 1440  809 (S18W19 L=  4)   282     94        0 (1411/ 776-1471/ 857)  31.0 W16-W21 S19-S16
     16  835  402 (S27E30 L=315)   477    147        0 ( 789/ 359- 900/ 443)  10.3 E35-E24 S28-S25
     17  635  280 (S30E54 L=291)  1632    235        2 ( 583/ 209- 658/ 384)  29.7 E67-E44 S35-S23
  

    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