Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1584 1653 (N46W54 L=256)  1690    232        4 (1485/1610-1691/1694)  -8.4 W41-W71 N42-N49
      2 1031 1560 (N38E01 L=201)   423     77        1 ( 988/1538-1073/1589)  21.8 E04-W01 N35-N39
      3  330 1502 (N30E61 L=141)   819     97        2 ( 297/1470- 378/1525)  -4.8 E66-E54 N27-N31
      4 1096 1483 (N32W04 L=206)   334     79        0 (1070/1437-1118/1516) -68.2 W01-W05 N29-N34
      5  500 1490 (N30E42 L=160)   300     62        0 ( 465/1476- 528/1507)  26.6 E44-E39 N29-N31
      6 1090 1343 (N23W03 L=206)   420    119        1 (1064/1291-1137/1377)  45.0 W01-W06 N19-N25
      7  184 1166 (N08E65 L=137)   719    130        0 ( 171/1104- 205/1230) -49.4 E68-E61 N03-N11
      8  757 1194 (N12E17 L=185)   193     59        1 ( 731/1174- 792/1215) -21.8 E19-E14 N11-N13
      9 1484 1121 (N10W29 L=231)  1018    255        1 (1411/1008-1534/1205) -59.7 W24-W32 N02-N14
     10 1645 1065 (N06W41 L=244)   485     97        2 (1616/1030-1679/1113) -51.3 W39-W44 N04-N09
     11 1725  814 (S09W50 L=252)  1079    189        2 (1681/ 742-1763/ 899)  68.2 W46-W52 S13-S03
     12  201  843 (S12E62 L=140)   508     50        0 ( 179/ 822- 223/ 864) -45.0 E65-E59 S13-S10
     13 1113  760 (S13W06 L=208)   211     61        0 (1087/ 733-1141/ 789) -36.9 W04-W08 S14-S11
     14  190  733 (S19E68 L=134)   918     71        0 ( 166/ 693- 209/ 772)  36.9 E75-E63 S22-S16
     15  692  605 (S25E22 L=181)  1081    209        1 ( 620/ 569- 786/ 659)  -9.5 E26-E14 S27-S21
     16 1638  598 (S23W47 L=249)  1478    166        3 (1581/ 539-1696/ 659)  53.1 W41-W51 S27-S18
     17  494  594 (S27E37 L=165)   324     84        0 ( 451/ 578- 527/ 620)  14.0 E41-E33 S27-S24
     18 1170  219 (S55W20 L=222)  2626    348        4 (1035/ 188-1313/ 250)   3.8 W05-W36 S59-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