Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1370 1751 (N48W48 L=202)   428     61        0 (1339/1720-1400/1777)  14.0 W42-W54 N46-N49
      2 1497 1753 (N45W63 L=217)   562     62        0 (1464/1733-1526/1772)   3.8 W56-W71 N44-N46
      3  646 1694 (N56E29 L=125)   449     56        0 ( 626/1659- 668/1723) -80.5 E30-E26 N52-N58
      4  418 1670 (N55E64 L= 90)   916     96        2 ( 370/1649- 465/1691)  -3.6 E83-E51 N53-N56
      5  588 1355 (N30E27 L=126)  1117    219        2 ( 489/1325- 691/1377)  -3.6 E35-E19 N28-N31
      6  496 1107 (N15E34 L=120)   597    190        0 ( 454/1034- 531/1187)  77.5 E37-E31 N11-N20
      7  928  977 (N04E06 L=147)  1186    417        0 ( 844/ 805- 975/1136) -78.7 E10-E03 S06-N14
      8 1621  959 (S06W38 L=191)   605    103        1 (1589/ 917-1659/1003) -56.3 W35-W40 S08-S02
      9  713  768 (S07E22 L=132)   458     85        1 ( 679/ 736- 753/ 793)  21.8 E24-E19 S08-S05
     10  852  763 (S08E13 L=140)   344    130        0 ( 791/ 730- 909/ 781) -26.6 E17-E10 S10-S06
     11 1819  718 (S25W59 L=212)  1010     95        6 (1797/ 668-1842/ 766) -74.1 W55-W61 S28-S21
     12 1279  619 (S22W12 L=166)   973    175        4 (1224/ 556-1332/ 688)  41.2 W07-W15 S26-S18
     13 1749  650 (S28W51 L=205)   443     57        0 (1734/ 618-1765/ 678)  76.0 W49-W53 S30-S26
     14  727  617 (S16E24 L=130)   418     75        3 ( 695/ 591- 766/ 653) -36.9 E26-E21 S18-S14
     15 1432  482 (S35W23 L=177)  3015    486        4 (1275/ 398-1582/ 626)  20.0 W11-W34 S40-S26
     16  814  467 (S27E21 L=132)   473    115        0 ( 774/ 445- 868/ 496) -26.6 E23-E17 S29-S25
     17  603  376 (S33E42 L=112)  1133    204        3 ( 551/ 296- 673/ 444)  31.0 E53-E32 S39-S27
     18 1454  279 (S54W32 L=186)  1018    125        2 (1380/ 259-1507/ 300)  -9.5 W22-W41 S57-S51
  

    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