Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  606 1484 (N31E25 L=192)   442    118        1 ( 562/1452- 650/1509) -35.5 E28-E21 N28-N33
      2 1590 1484 (N22W44 L=261)   313     91        0 (1554/1469-1626/1506)   8.1 W40-W47 N21-N23
      3 1246 1373 (N18W17 L=234)   344    158        2 (1190/1347-1307/1402) -23.2 W13-W20 N15-N19
      4  466 1322 (N21E35 L=183)   317     73        2 ( 446/1295- 489/1360) -68.2 E36-E32 N19-N24
      5 1472 1162 (N03W29 L=246)   429    106        0 (1454/1102-1491/1214)  80.5 W27-W30 N00-N06
      6 1731 1111 (S02W48 L=265)   888    132        2 (1692/1080-1786/1164) -29.7 W44-W52 S03-N01
      7  510 1086 (N07E32 L=186)   358     44        0 ( 481/1060- 535/1111)  26.6 E33-E29 N05-N08
      8 1663 1025 (S06W41 L=259)   231     51        1 (1643/1003-1679/1046) -33.7 W39-W42 S07-S04
      9  344  962 (N01E45 L=172)   288     73        1 ( 307/ 940- 364/ 987)  21.8 E48-E43 N00-N02
     10 1737  931 (S12W47 L=265)   714    154        2 (1705/ 892-1767/ 986) -80.5 W44-W49 S14-S08
     11 1557  852 (S16W33 L=250)   892    258        1 (1453/ 726-1619/ 938)  47.7 W25-W37 S22-S10
     12 1555  744 (S22W33 L=250)  1305    343        2 (1442/ 593-1681/ 903)  46.6 W25-W42 S31-S13
     13  312  757 (S11E51 L=166)  3066    572        9 ( 198/ 570- 419/ 895) -74.1 E62-E44 S22-S01
     14  667  843 (S09E23 L=194)   220     52        0 ( 652/ 806- 681/ 875) -90.0 E24-E22 S11-S07
     15  649  561 (S26E30 L=187)   455     53        0 ( 614/ 534- 680/ 583) -11.3 E32-E27 S27-S24
     16 1573  417 (S46W41 L=258)   844    119        1 (1551/ 377-1600/ 470)  90.0 W38-W44 S48-S41
  

    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