Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1395 1690 (N33W45 L=296)   297     40        0 (1370/1668-1420/1707)   9.5 W41-W47 N32-N34
      2 1460 1327 (N10W32 L=284)   226     54        1 (1433/1312-1490/1342)   0.0 W30-W34 N09-N11
      3  247 1282 (N32E53 L=199)   767    124        2 ( 214/1226- 283/1330) -35.0 E57-E47 N28-N35
      4  377 1252 (N28E39 L=213)   407     99        1 ( 345/1209- 411/1300) -54.5 E41-E35 N24-N31
      5  182 1184 (N27E60 L=191)   476     95        2 ( 143/1144- 217/1223)  16.7 E65-E55 N25-N29
      6 1347 1191 (N05W22 L=273)   611    210        2 (1246/1166-1423/1220) -10.3 W15-W26 N02-N07
      7  716 1093 (N12E16 L=235)   347    115        0 ( 688/1067- 763/1143) -80.5 E17-E13 N09-N15
      8  916 1028 (N04E06 L=245)   307     89        0 ( 868/ 993- 954/1058) -45.0 E08-E04 N01-N06
      9  233 1002 (N15E54 L=197)   646    196        2 ( 161/ 954- 293/1045) -12.1 E62-E48 N11-N18
     10 1855  763 (S31W58 L=309)  2463    390        0 (1732/ 614-1971/ 884)  18.4 W44-W80 S38-S26
     11  468  744 (S03E40 L=212)   421    166        1 ( 398/ 710- 514/ 769) -20.6 E45-E36 S06-S01
     12  502  625 (S11E41 L=210)  1289    369        4 ( 383/ 575- 642/ 705) -29.4 E48-E31 S14-S04
     13 1188  471 (S35E02 L=250)  1122    466        2 (1011/ 363-1379/ 591)  13.1 E17-W13 S37-S30
     14  472  519 (S16E48 L=204)  1046    166        2 ( 417/ 482- 539/ 576) -54.5 E51-E43 S19-S12
     15  787  498 (S24E27 L=225)   332    144        0 ( 776/ 428- 808/ 567)  52.1 E29-E22 S28-S19
     16  376  482 (S17E59 L=193)  1279    211        3 ( 334/ 401- 414/ 552)  68.2 E62-E53 S22-S12
     17  622  234 (S38E57 L=194)  1150    132        0 ( 588/ 190- 645/ 303)  19.2 E70-E47 S41-S33
  

    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