Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  786 1549 (N22E23 L=346)   249     74        1 ( 770/1513- 802/1586) -53.1 E24-E21 N20-N24
      2  484 1518 (N17E44 L=325)   837    213        3 ( 408/1495- 595/1539)   3.8 E50-E35 N15-N19
      3  893 1458 (N18E14 L=355)   866    235        3 ( 796/1416- 995/1498)   0.0 E20-E07 N15-N19
      4  834 1370 (N12E16 L=353)   710    207        4 ( 773/1307- 886/1445)  77.5 E19-E13 N07-N16
      5  942 1373 (N13E10 L=359)   243     85        1 ( 915/1345- 980/1392)  26.6 E11-E07 N11-N14
      6 1293 1068 (N00W15 L= 24)   182     64        1 (1264/1041-1314/1091)   0.0 W13-W16 S01-N01
      7 1154 1068 (S02W07 L= 16)   188     91        1 (1104/1059-1196/1076)  11.3 W04-W09 S02-N00
      8  645 1020 (S11E23 L=346)   404    142        1 ( 593/ 968- 678/1074)  66.8 E25-E20 S14-S07
      9 1524 1005 (N00W30 L= 39)   269    111        0 (1512/ 950-1537/1059) -78.7 W29-W31 S02-N03
     10  172  903 (S22E63 L=306)  1628    217        3 (  91/ 834- 223/ 966)  18.4 E80-E56 S25-S17
     11  594  786 (S25E25 L=344)  1369    373        2 ( 435/ 741- 768/ 840)   7.1 E36-E12 S29-S21
     12 1329  686 (S21W24 L= 33)  1109    431        2 (1210/ 566-1495/ 815)  55.6 W19-W32 S30-S10
     13 1066  768 (S20W06 L= 15)   360    199        0 ( 980/ 740-1154/ 786)  15.3 W00-W11 S21-S18
     14 1553  636 (S20W41 L= 50)   741    170        2 (1517/ 570-1576/ 721) -66.0 W37-W43 S24-S15
     15  992  554 (S35W05 L= 14)   894    238        0 ( 870/ 534-1109/ 574)   7.1 E03-W13 S36-S32
     16 1409  287 (S45W52 L= 60)  1880    230        3 (1314/ 244-1494/ 338)  12.5 W39-W63 S46-S40
  

    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