Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  631 1803 (N31E53 L=122)   751    126        3 ( 601/1782- 660/1833)  -8.1 E56-E49 N29-N33
      2 1315 1657 (N40E00 L=176)   412     39        0 (1284/1627-1343/1683) -26.6 E02-W02 N37-N41
      3  661 1628 (N19E40 L=136)   477     75        0 ( 622/1604- 708/1645)  31.0 E42-E36 N17-N20
      4  754 1516 (N14E29 L=147)   573    188        0 ( 727/1449- 785/1602) -48.4 E33-E25 N10-N19
      5  913  981 (S12E05 L=171)   181     96        0 ( 864/ 968- 954/1007)  36.9 E07-E03 S13-S09
      6 1526  982 (N04W30 L=206)   359    103        1 (1480/ 963-1580/ 998)  33.7 W27-W33 N01-N05
      7 1614  755 (S05W43 L=219)  1543    409        6 (1563/ 618-1652/ 907) -66.4 W37-W47 S13-N03
      8  747  814 (S26E11 L=165)   670    200        0 ( 635/ 783- 836/ 848)  24.8 E18-E05 S29-S21
      9 1423  677 (S15W33 L=209)  1097    226        3 (1328/ 626-1523/ 747)  56.3 W28-W37 S20-S08
     10 1278  598 (S24W28 L=203)   844    200        0 (1198/ 557-1366/ 651)   0.0 W20-W34 S24-S22
     11 1033  439 (S40W20 L=196)  1913    436        5 ( 852/ 406-1222/ 478)  14.5 W06-W33 S44-S35
     12 1540  378 (S27W57 L=233)  2291    418        3 (1374/ 349-1662/ 409)  25.6 W44-W67 S32-S21
  

    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