Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  920 1615 (N41W09 L=  8)   359     68        2 ( 887/1596- 959/1637) -36.9 W06-W12 N38-N42
      2 1897 1225 (S07W68 L= 66)  1525    209        7 (1866/1137-1924/1304)  52.1 W64-W74 S12-S02
      3  572 1268 (N27E23 L=336)   222     73        0 ( 546/1247- 606/1298) -45.0 E24-E20 N24-N29
      4 1775 1203 (S05W53 L= 51)   357     58        2 (1745/1175-1805/1225) -18.4 W49-W55 S06-S03
      5 1486 1106 (S03W28 L= 27)  1720    357        2 (1376/ 945-1533/1219) -80.0 W23-W31 S12-N05
      6 1691 1176 (S04W45 L= 43)   592    100        4 (1638/1148-1734/1208) -45.0 W40-W48 S06-S01
      7  540 1162 (N21E26 L=332)   661    188        1 ( 432/1139- 612/1188)  -8.7 E34-E21 N18-N22
      8 1797 1088 (S12W52 L= 50)  1428    224        1 (1758/1004-1838/1171)  41.6 W46-W57 S15-S07
      9  872 1118 (N11E07 L=352)   442    141        1 ( 804/1095- 941/1151) -35.5 E09-E02 N09-N14
     10 1710 1075 (S10W44 L= 42)   798    223        3 (1658/1011-1766/1148)  22.6 W38-W50 S13-S07
     11 1848  850 (S27W55 L= 53)  1216    133        4 (1793/ 823-1907/ 893) -20.6 W48-W64 S30-S23
     12 1712  598 (S40W41 L= 39)  6097    474       14 (1630/ 439-1809/ 856) -86.1 W33-W56 S52-S21
     13 1083  784 (S12E02 L=356)   834    173        5 (1008/ 744-1147/ 840)  -7.1 E07-W01 S14-S09
     14 1136  646 (S21E02 L=356)   974    279        0 (1034/ 574-1192/ 757)  33.7 E10-W02 S24-S15
     15  641  643 (S10E32 L=327)   406    142        1 ( 570/ 611- 692/ 683)   5.7 E37-E27 S11-S08
     16  759  447 (S25E32 L=327)  2511    557        6 ( 609/ 254- 941/ 593)  21.0 E54-E15 S35-S19
     17  408  424 (S19E61 L=298)   534     43        1 ( 389/ 393- 426/ 454)  38.7 E63-E57 S20-S16
  

    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