Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1528 1701 (N28W63 L=259)   636     84        1 (1505/1668-1555/1739)  56.3 W59-W65 N24-N30
      2 1302 1630 (N32W39 L=236)  1515    380        1 (1138/1583-1478/1690) -46.7 W31-W49 N23-N40
      3  990 1601 (N40W17 L=214)   364     81        0 ( 947/1579-1025/1618) -18.4 W14-W20 N38-N41
      4 1469 1456 (N17W40 L=237)   769    317        0 (1309/1434-1593/1482) -25.3 W29-W48 N11-N21
      5 1077 1342 (N22W12 L=209)   650    200        1 ( 994/1318-1184/1376) -18.4 W07-W19 N20-N24
      6  633 1171 (N24E20 L=177)   265     65        1 ( 597/1145- 658/1189) -45.0 E21-E18 N22-N26
      7  354 1128 (N28E42 L=155)   286     77        0 ( 336/1092- 374/1165) -68.2 E43-E40 N25-N30
      8 1862 1025 (S19W58 L=255)   412     64        1 (1841/ 988-1881/1059)  26.6 W54-W61 S20-S17
      9  502 1001 (N18E32 L=164)   940    344        1 ( 389/ 967- 666/1046) -27.8 E40-E21 N11-N22
     10  354  949 (N18E45 L=152)   283     19        0 ( 328/ 931- 374/ 964) -26.6 E47-E43 N16-N18
     11 1763  918 (S22W45 L=242)   355     91        0 (1719/ 906-1807/ 931) -12.5 W40-W49 S23-S20
     12  329  812 (N10E50 L=146)   292     46        1 ( 315/ 785- 347/ 838)  56.3 E52-E48 N08-N11
     13 1483  799 (S19W21 L=217)   228     65        0 (1462/ 763-1502/ 832)  36.9 W18-W22 S20-S17
     14  350  748 (N06E51 L=146)   486    122        1 ( 307/ 723- 408/ 770) -26.6 E54-E46 N03-N08
     15 1119  675 (S16E04 L=193)   453    168        2 (1080/ 624-1159/ 757) -90.0 E05-E02 S20-S10
     16 1380  711 (S22W12 L=209)   236     86        0 (1359/ 668-1398/ 753)  45.0 W09-W14 S23-S19
     17  489  302 (S24E67 L=130)   860     74        1 ( 463/ 264- 517/ 337)  45.0 E71-E62 S27-S21
     18 1342  297 (S51E03 L=194)   767     82        0 (1298/ 271-1388/ 334)   0.0 E06-W02 S53-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