Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1628 1654 (N34W68 L=315)   928     80        3 (1589/1630-1676/1675)  -3.6 W61-W78 N32-N34
      2  385 1552 (N48E55 L=192)  2158    293        4 ( 339/1433- 425/1657) -23.0 E87-E43 N38-N55
      3  780 1569 (N44E11 L=236)   271     54        0 ( 752/1548- 810/1590) -38.7 E13-E08 N42-N46
      4 1028 1414 (N30W06 L=253)   274     92        1 ( 976/1400-1071/1429) -26.6 W02-W09 N28-N31
      5  610 1404 (N35E25 L=222)  1499    405        3 ( 432/1366- 785/1428)  -9.8 E40-E11 N30-N37
      6 1293 1279 (N17W21 L=267)   525    101        3 (1256/1220-1324/1320)  71.6 W18-W22 N13-N19
      7 1640 1263 (N11W45 L=292)   644    195        0 (1546/1217-1715/1288) -30.3 W38-W50 N06-N13
      8  172  998 (N12E66 L=181)  2746    291        7 (  89/ 907- 235/1113)  77.0 E89-E57 N06-N19
      9 1913 1006 (S11W70 L=317)  1825    122        3 (1889/ 944-1943/1071) -63.4 W65-W77 S14-S06
     10  655  702 (S09E28 L=219)  1834    482        5 ( 558/ 584- 795/ 848)  28.8 E36-E16 S16-S01
     11  217  656 (S10E70 L=177)  2602    272        3 ( 156/ 548- 283/ 759) -69.0 E78-E63 S17-S03
     12 1755  717 (S25W50 L=297)   442     79        1 (1731/ 679-1778/ 755) -50.2 W46-W53 S28-S22
     13  519  673 (S10E38 L=209)   263     61        2 ( 492/ 654- 542/ 699) -33.7 E39-E36 S11-S07
     14  348  650 (S10E53 L=194)   368     58        1 ( 323/ 632- 377/ 670) -36.9 E56-E50 S11-S08
     15 1656  643 (S28W40 L=287)   389     82        1 (1613/ 627-1690/ 670)   0.0 W36-W43 S29-S26
     16  363  480 (S22E61 L=186)  2179    281        4 ( 327/ 368- 392/ 607)  50.2 E75-E54 S31-S12
     17  710  517 (S22E28 L=218)   905    189        5 ( 649/ 469- 760/ 570)  26.6 E34-E24 S24-S18
     18 1714  474 (S44W55 L=302)   995     73        4 (1675/ 446-1743/ 508) -11.3 W48-W63 S46-S40
     19  994  165 (S57E25 L=222)   546     41        1 ( 969/ 143-1020/ 194)  33.7 E30-E21 S60-S53
  

    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