Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1079 1683 (N48W20 L=327)   474     63        1 (1051/1651-1117/1705)  14.0 W16-W24 N46-N49
      2 1191 1588 (N39W24 L=331)   765    116        1 (1157/1543-1235/1638) -82.9 W21-W25 N34-N42
      3 1577 1522 (N26W52 L=359)   465     93        0 (1534/1494-1618/1553) -51.3 W49-W54 N23-N28
      4 1100 1511 (N35W14 L=321)   224     70        0 (1066/1480-1121/1528) -45.0 W12-W15 N33-N36
      5 1420 1379 (N21W33 L=340)   223     41        0 (1390/1368-1447/1389) -14.0 W30-W34 N19-N21
      6  170 1151 (N23E69 L=238)   912    143        2 ( 145/1070- 183/1211) -48.8 E75-E65 N18-N26
      7  511 1108 (N18E33 L=274)   787    231        0 ( 410/1048- 612/1165) -31.0 E40-E25 N13-N22
      8 1723  860 (S16W45 L=352)   596    205        0 (1644/ 795-1781/ 932)  21.0 W38-W51 S18-S12
      9 1504  856 (S12W28 L=335)   743    176        2 (1466/ 791-1548/ 928)  53.1 W24-W31 S15-S07
     10 1806  837 (S20W54 L=  1)  1281    165        3 (1786/ 755-1821/ 922) -71.6 W51-W56 S25-S13
     11 1006  798 (S07E04 L=303)   773    196        1 ( 966/ 703-1055/ 875)  59.0 E08-W00 S12-S02
     12 1500  640 (S26W26 L=333)   258     90        0 (1473/ 604-1522/ 681)  45.0 W23-W27 S27-S23
  

    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