Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  926 1593 (N41E08 L=353)   297     30        1 ( 896/1575- 951/1608)   0.0 E10-E05 N39-N42
      2 1486 1384 (N26W33 L= 34)   326     89        0 (1460/1341-1512/1425) -80.5 W31-W34 N22-N28
      3  584 1276 (N19E30 L=331)   980    292        0 ( 514/1171- 660/1398) -45.0 E37-E23 N12-N26
      4 1197 1172 (N13W11 L= 11)   828    208        0 (1122/1107-1266/1245) -45.0 W06-W15 N08-N17
      5  978 1105 (N09E03 L=357)   353    165        0 ( 959/1034- 996/1174)  90.0 E04-E01 N04-N12
      6 1690 1044 (N04W46 L= 46)   403     71        1 (1668/1006-1707/1085)  68.2 W43-W47 N01-N06
      7 1531  899 (S05W33 L= 33)   258     39        2 (1509/ 874-1548/ 922)  90.0 W31-W34 S06-S03
      8 1214  870 (S06W12 L= 12)   197     67        1 (1196/ 827-1223/ 904) -78.7 W10-W12 S08-S03
      9 1077  793 (S11W03 L=  4)   811    217        5 (1008/ 721-1165/ 861)  32.0 E01-W08 S15-S06
     10  959  748 (S14E04 L=356)   305    124        1 ( 936/ 690- 980/ 796)  81.9 E05-E02 S17-S10
  

    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