Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1301 1596 (N40W26 L=115)   700    119        3 (1238/1552-1361/1636) -36.9 W22-W30 N36-N43
      2 1515 1487 (N30W40 L=128)   366     83        0 (1475/1465-1556/1509) -31.0 W37-W42 N28-N32
      3 1682 1388 (N22W51 L=139)   542     70        1 (1654/1353-1710/1426) -68.2 W49-W53 N19-N24
      4 1769 1390 (N21W60 L=148)   688    113        1 (1738/1357-1800/1425)  38.7 W56-W64 N19-N23
      5 1797  770 (S17W56 L=144)  1004    165        3 (1745/ 703-1826/ 822)  90.0 W50-W59 S21-S14
      6 1737  696 (S22W51 L=139)   564     85        3 (1708/ 669-1758/ 731)  45.0 W47-W52 S23-S19
      7  454  689 (S14E41 L= 47)   372    124        1 ( 411/ 651- 500/ 716) -45.0 E44-E37 S16-S12
      8  938  664 (S18E07 L= 81)   216     73        1 ( 901/ 645- 977/ 689) -36.9 E09-E05 S18-S15
      9 1759  583 (S30W56 L=144)  1316    158        0 (1698/ 523-1810/ 653)  45.0 W50-W59 S34-S25
     10 1306  577 (S26W17 L=105)   696    155        1 (1255/ 532-1387/ 614)  29.1 W13-W22 S28-S23
     11 1059  374 (S40E02 L= 87)   977    188        2 ( 961/ 334-1123/ 426)  30.3 E09-W03 S42-S35
  

    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