Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  468 1622 (N19E56 L=254)  1958    379        6 ( 402/1513- 516/1740) -64.8 E63-E52 N10-N27
      2  381 1542 (N12E59 L=250)   617     94        0 ( 344/1500- 423/1586)  66.8 E61-E57 N08-N15
      3  493 1494 (N11E47 L=263)   702    126        2 ( 455/1447- 546/1548)  74.1 E50-E44 N08-N15
      4  350 1344 (S00E53 L=257)   498     98        2 ( 323/1298- 384/1382)  78.7 E54-E50 S03-N02
      5  132 1187 (S13E73 L=237)  1939    168        1 ( 105/1112- 153/1253) -69.4 E77-E69 S17-S09
      6  464 1077 (S13E36 L=274)   588    107        1 ( 432/1022- 503/1119)  74.1 E37-E33 S17-S10
      7  445  820 (S29E33 L=277)   250     62        0 ( 428/ 790- 468/ 840)  53.1 E34-E31 S31-S27
      8  985  749 (S20W05 L=314)   400    128        0 ( 934/ 714-1051/ 776)  45.0 W01-W08 S23-S17
      9  321  678 (S42E44 L=265)   910    249        1 ( 260/ 631- 385/ 717)  28.1 E52-E37 S46-S37
     10 1716  568 (S11W61 L= 11)   918    244        5 (1662/ 539-1776/ 593)  24.4 W55-W69 S13-S08
     11 1423  435 (S27W44 L=354)   605    163        2 (1354/ 412-1488/ 469)  10.3 W38-W49 S27-S25
     12 1369  261 (S39W53 L=  2)   999    226        2 (1343/ 208-1389/ 325) -38.7 W48-W59 S42-S34
  

    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