Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  576 1694 (N57E37 L= 75)   674     97        1 ( 527/1671- 620/1712)  -7.6 E46-E30 N54-N58
      2  399 1307 (N29E43 L= 69)   345     99        0 ( 368/1260- 421/1340) -45.0 E46-E40 N26-N31
      3 1133 1248 (N18W10 L=122)   321    119        0 (1101/1192-1149/1301) -90.0 W08-W10 N14-N21
      4 1591 1112 (N03W37 L=149)   742    236        1 (1554/1010-1613/1216)  84.8 W35-W39 S02-N09
      5 1796  691 (S27W56 L=168)  3748    391        5 (1759/ 524-1844/ 863) -55.8 W50-W74 S41-S16
      6  413  599 (S15E48 L= 64)   459     71        0 ( 376/ 575- 453/ 618)  14.0 E52-E44 S15-S13
      7 1389  553 (S29W19 L=131)   869    160        2 (1321/ 509-1463/ 597)  10.3 W13-W24 S31-S26
      8 1072  409 (S35E05 L=107)  1314    383        3 ( 894/ 309-1205/ 496)  19.1 E21-W05 S41-S29
      9  735  266 (S42E38 L= 74)   556     71        0 ( 690/ 252- 774/ 285)   0.0 E43-E34 S43-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