Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1340 1830 (N50W62 L=165)  1243    150        1 (1267/1793-1397/1868)   0.0 W48-W81 N48-N51
      2 1750 1396 (N13W61 L=164)  1101    219        2 (1669/1353-1832/1451) -31.0 W55-W70 N08-N17
      3 1322 1340 (N20W25 L=128)   793    336        0 (1215/1259-1433/1406)  12.5 W16-W34 N16-N21
      4  481 1358 (N36E34 L= 69)   260    111        0 ( 450/1315- 509/1404)  51.3 E36-E31 N33-N38
      5 1875 1036 (S11W63 L=166)  6956    354       16 (1808/ 872-1949/1192)  81.5 W53-W80 S21-S01
      6 1725  842 (S18W44 L=147)   490     91        1 (1697/ 806-1753/ 878)  38.7 W41-W47 S20-S16
      7  223  830 (N05E62 L= 41)   445     58        1 ( 199/ 800- 249/ 857)  23.2 E65-E58 N03-N06
      8 1344  815 (S11W16 L=118)   298     78        1 (1323/ 790-1374/ 855) -68.2 W14-W17 S13-S08
      9 1454  809 (S14W23 L=125)   420     81        1 (1433/ 766-1478/ 851)  68.2 W20-W24 S16-S11
     10 1119  679 (S16E00 L=102)  1332    402        6 (1031/ 511-1186/ 815) -87.4 E03-W03 S27-S05
     11  540  691 (S06E37 L= 66)   316     61        0 ( 512/ 668- 577/ 717) -36.9 E38-E34 S07-S04
     12  229  601 (S09E72 L= 30)  1422    136        1 ( 191/ 546- 272/ 657) -58.0 E76-E68 S13-S05
     13  543  534 (S16E41 L= 61)   566    123        1 ( 526/ 474- 561/ 592)  69.4 E43-E39 S20-S12
     14 1280  483 (S32W07 L=110)   910    132        4 (1219/ 451-1330/ 531)   7.1 W02-W10 S34-S29
     15  841  437 (S27E24 L= 79)   661     99        1 ( 820/ 390- 866/ 481)  45.0 E27-E21 S30-S23
     16  366  438 (S21E66 L= 37)   831     73        1 ( 350/ 390- 383/ 477)  28.3 E73-E60 S25-S18
     17 1458  324 (S51W22 L=125)   736     94        1 (1416/ 296-1501/ 352)   0.0 W17-W26 S54-S49
  

    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