Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1441 1780 (N55W56 L=164)  1136    112        0 (1413/1735-1488/1815)  11.0 W47-W78 N51-N57
      2  531 1681 (N49E50 L= 57)   518     68        1 ( 501/1651- 557/1712)  90.0 E52-E47 N46-N51
      3 1474 1579 (N38W40 L=148)   298     45        0 (1451/1559-1498/1596)   9.5 W37-W43 N37-N39
      4  753 1431 (N31E19 L= 89)  2405    542        4 ( 544/1289- 943/1558) -28.5 E39-E04 N20-N39
      5 1290 1266 (N18W18 L=126)   528    212        0 (1257/1190-1356/1331)  52.1 W15-W23 N13-N22
      6 1392 1205 (N14W24 L=132)   309     84        3 (1367/1177-1418/1236)  33.7 W22-W26 N12-N15
      7  995 1005 (N03E02 L=106)   490    120        3 ( 949/ 970-1051/1050)  38.7 E04-W01 N00-N05
      8  844  664 (S18E13 L= 95)  2204    643        2 ( 642/ 429- 940/ 879)  50.2 E31-E05 S34-S04
      9 1496  726 (S16W31 L=139)   308    103        0 (1470/ 682-1530/ 773)  56.3 W29-W33 S19-S13
     10 1591  529 (S31W43 L=151)  4809    673        7 (1346/ 367-1782/ 735)  38.0 W25-W57 S42-S17
     11 1141  445 (S35W07 L=115)  1716    348        2 (1006/ 352-1267/ 535)  35.8 E03-W15 S41-S28
  

    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