Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  595 1755 (N70E17 L=108)  1354    105        4 ( 549/1711- 645/1788)   9.5 E26-E07 N66-N75
      2 1207 1718 (N42W40 L=165)   433     60        1 (1174/1697-1235/1738)  -9.5 W37-W43 N40-N43
      3 1588 1622 (N23W61 L=187)   433     84        0 (1555/1587-1620/1651) -90.0 W59-W63 N20-N25
      4 1382 1477 (N22W36 L=161)   458    212        0 (1278/1456-1460/1503) -28.6 W29-W40 N19-N25
      5  897 1359 (N29W02 L=127)   420    164        0 ( 845/1321- 972/1408) -58.0 W00-W05 N25-N33
      6  348 1305 (N40E44 L= 81)   543     74        3 ( 319/1266- 375/1339) -45.0 E47-E40 N36-N42
      7  814 1256 (N25E07 L=119)   644    193        0 ( 735/1225- 911/1301) -35.0 E10-W00 N21-N30
      8  528  940 (N14E32 L= 93)   720    146        2 ( 467/ 918- 603/ 967) -21.8 E36-E26 N11-N16
      9  366  908 (N15E45 L= 80)   320     78        1 ( 340/ 879- 390/ 943)  21.8 E48-E43 N13-N17
     10  771  912 (N07E17 L=108)   260     94        2 ( 748/ 864- 798/ 942) -90.0 E18-E15 N03-N09
     11  213  896 (N17E62 L= 64)   534     74        1 ( 191/ 862- 236/ 929) -53.1 E64-E58 N14-N18
     12 1308  779 (S14W10 L=135)   366     89        0 (1290/ 719-1324/ 822)  63.4 W08-W11 S17-S11
     13  285  755 (N07E57 L= 68)   962    116        1 ( 237/ 729- 346/ 784) -21.8 E62-E51 N05-N09
     14 1592  632 (S33W27 L=152)   461     61        1 (1568/ 600-1611/ 668) -90.0 W24-W28 S34-S29
     15  731  608 (S10E28 L= 97)   287     69        0 ( 714/ 568- 745/ 643)  68.2 E29-E27 S12-S07
     16  532  490 (S13E47 L= 79)   863    159        3 ( 507/ 424- 557/ 563)  56.3 E50-E43 S17-S08
     17  444  459 (S13E56 L= 69)   375     39        0 ( 427/ 434- 463/ 481)  45.0 E58-E53 S14-S11
     18  513  366 (S21E57 L= 68)   684    115        0 ( 456/ 347- 573/ 393)   0.0 E64-E49 S21-S19
  

    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