Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  690 1627 (N38E34 L=343)  1051    167        2 ( 619/1592- 771/1657)  -3.8 E40-E25 N36-N40
      2  300 1521 (N27E69 L=307)   604     68        0 ( 272/1497- 326/1539)  -4.4 E76-E63 N25-N27
      3 1255 1395 (N27W13 L= 29)  1406    350        2 (1177/1304-1313/1510) -69.0 W06-W16 N20-N34
      4  884 1474 (N29E14 L=  3)   252     41        1 ( 861/1454- 905/1497)  45.0 E15-E12 N27-N31
      5  440 1364 (N18E45 L=332)   570    173        0 ( 385/1296- 479/1405) -24.4 E50-E39 N14-N20
      6  877 1121 (N07E10 L=  7)   787    253        2 ( 804/1052- 988/1185)  45.0 E13-E03 N01-N11
      7  664  867 (S11E22 L=355)   207     52        0 ( 647/ 849- 694/ 890) -33.7 E22-E19 S11-S09
      8 1248  568 (S25W19 L= 36)   377    101        1 (1196/ 548-1297/ 591)  -8.1 W15-W22 S26-S23
      9 1115  386 (S40W14 L= 30)   308     64        1 (1079/ 372-1148/ 409)  18.4 W10-W16 S41-S38
  

    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