Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1201 1748 (N52W01 L=160)   462     45        1 (1171/1722-1231/1767)  26.6 E02-W03 N50-N54
      2  599 1727 (N39E51 L=108)   518     63        1 ( 563/1712- 644/1739)  15.9 E54-E46 N37-N40
      3  906 1522 (N29E16 L=143)  2038    425        3 ( 731/1424-1093/1599)  -8.1 E29-E01 N25-N32
      4 1709 1548 (N45W55 L=215)   882     84        1 (1667/1531-1755/1576)   9.5 W47-W65 N42-N47
      5 1216 1324 (N21W08 L=167)   237     73        0 (1182/1301-1249/1351) -31.0 W05-W10 N19-N22
      6 1123 1249 (N15W03 L=162)   472    217        0 (1051/1157-1176/1297) -32.0 E01-W07 N10-N17
      7  193 1225 (N01E65 L= 94)  1485    131        4 ( 160/1163- 239/1283)  45.0 E70-E60 S03-N04
      8 1733 1144 (N17W47 L=207)   358     46        0 (1710/1124-1764/1158)  21.8 W45-W50 N15-N18
      9 1149  793 (S12W11 L=170)   661    235        1 (1061/ 749-1252/ 852)  32.5 W05-W16 S15-S07
     10  489  782 (S22E32 L=127)   311     79        2 ( 456/ 759- 516/ 812)  45.0 E34-E30 S23-S19
     11 1222  626 (S21W19 L=178)  2055    605        2 ( 953/ 572-1404/ 769)  34.7 W01-W27 S27-S09
     12  449  589 (S36E37 L=122)   297     43        0 ( 435/ 566- 467/ 616)  63.4 E38-E35 S37-S33
     13 1130  253 (S50W27 L=187)   950    190        1 (1044/ 210-1212/ 291)   0.0 W17-W38 S53-S48
  

    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