Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  412 1487 (N31E51 L=235)   377     75        0 ( 386/1460- 450/1518) -18.4 E54-E45 N29-N33
      2 1646 1455 (N30W50 L=335)  1015    156        2 (1591/1383-1707/1505) -58.0 W44-W53 N25-N33
      3  396 1020 (N02E42 L=243)   256     46        0 ( 377/ 992- 416/1047) -45.0 E43-E40 N00-N03
      4 1773  969 (S01W53 L=338)   478    107        1 (1736/ 944-1809/1014)  14.0 W49-W57 S02-N01
      5 1706  963 (S01W47 L=332)   943    123        1 (1672/ 896-1741/1006) -71.6 W43-W50 S05-N01
      6  163  873 (S09E68 L=217)  1740    174        3 ( 133/ 813- 192/ 975) -78.7 E76-E64 S12-S02
      7 1838  856 (S08W62 L=347)   735     84        3 (1802/ 825-1864/ 886)  26.6 W57-W65 S10-S06
      8  266  793 (S13E56 L=229)  1101    217        5 ( 222/ 735- 348/ 852) -28.6 E60-E48 S16-S09
      9  215  680 (S21E67 L=218)  1233    104        0 ( 173/ 648- 254/ 730)  15.9 E81-E59 S23-S17
     10 1214  688 (S18W12 L=298)   226     56        0 (1196/ 660-1232/ 720) -63.4 W11-W13 S19-S15
     11 1581  656 (S20W40 L=325)   342     52        1 (1554/ 640-1606/ 677)  14.0 W37-W41 S21-S18
     12 1711  609 (S24W54 L=339)   374     59        1 (1692/ 576-1726/ 635)  90.0 W52-W55 S26-S22
     13  872  565 (S26E10 L=275)   358    109        0 ( 826/ 534- 912/ 606)  39.8 E13-E07 S28-S23
     14 1779  550 (S29W68 L=353)   850    116        0 (1739/ 517-1817/ 593)  32.0 W62-W72 S30-S25
     15 1146  516 (S29W09 L=294)  1657    341        4 ( 983/ 462-1289/ 578)  19.3 E02-W18 S33-S24
     16 1508  411 (S38W42 L=327)   401     79        0 (1466/ 399-1545/ 432)  15.9 W37-W45 S39-S36
  

    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