Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1428 1773 (N55W53 L=163)  1130    112        1 (1398/1731-1471/1813)  13.0 W45-W71 N51-N57
      2  518 1676 (N49E52 L= 58)   504     64        1 ( 490/1651- 545/1705)  68.2 E54-E48 N46-N51
      3 1446 1564 (N38W37 L=147)   307     69        1 (1417/1540-1467/1589)  18.4 W33-W39 N35-N39
      4  716 1435 (N31E21 L= 88)  2419    551        1 ( 523/1283- 909/1567) -29.7 E42-E06 N20-N40
      5 1727 1483 (N30W63 L=173)   447     31        1 (1708/1460-1743/1512)  26.6 W59-W66 N28-N31
      6 1259 1263 (N18W16 L=126)   557    221        0 (1226/1183-1328/1329)  52.1 W13-W21 N13-N22
      7  330 1151 (N12E49 L= 61)   287     70        0 ( 313/1119- 358/1170)  45.0 E50-E46 N10-N13
      8  964 1009 (N03E04 L=106)   584    157        3 ( 911/ 959-1034/1061)  45.0 E07-W00 N00-N06
      9 1655  938 (S04W42 L=152)   334     79        2 (1619/ 913-1686/ 971) -36.9 W39-W45 S05-S01
     10  826  665 (S18E14 L= 96)  1634    378        3 ( 705/ 496- 902/ 798)  53.1 E24-E08 S29-S09
     11 1459  718 (S17W28 L=138)   277    102        0 (1425/ 675-1490/ 760)  51.3 W26-W30 S19-S13
     12 1561  524 (S31W40 L=150)  4788    687        8 (1285/ 357-1763/ 734)  35.5 W19-W55 S42-S17
     13 1147  470 (S33W07 L=117)  1935    523        1 ( 979/ 352-1368/ 621)  35.1 E05-W22 S41-S22
  

    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