Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1521 1602 (N32W49 L= 33)   539    150        0 (1470/1577-1553/1640)   5.7 W44-W54 N30-N33
      2 1216 1473 (N30W21 L=  5)   494    141        1 (1144/1435-1282/1502)  -5.7 W16-W26 N28-N32
      3 1739 1415 (N15W60 L= 44)  1484    231        3 (1673/1384-1804/1442) -26.6 W53-W66 N11-N18
      4 1229 1192 (N13W15 L=359)   202     85        0 (1183/1167-1261/1210) -36.9 W12-W16 N11-N14
      5 1312  949 (S03W16 L=  0)  2206    511        7 (1253/ 740-1358/1156)  71.6 W09-W19 S14-N10
      6  695 1088 (N15E20 L=324)   770    157        2 ( 624/1046- 771/1131) -33.7 E23-E14 N11-N18
      7 1585  950 (S08W34 L= 18)   719    203        3 (1528/ 889-1663/1000) -48.0 W30-W39 S13-S03
      8 1482  711 (S21W24 L=  8)  1321    357        2 (1446/ 581-1516/ 899)  79.2 W21-W27 S29-S08
      9 1314  539 (S29W11 L=355)   267    118        1 (1262/ 523-1373/ 579)   7.1 W07-W15 S30-S27
     10 1085  374 (S37E08 L=336)  1657    364        4 ( 966/ 284-1158/ 501)  34.4 E19-W00 S42-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