Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  793 1521 (N31E12 L=191)   777    205        1 ( 706/1474- 882/1558) -28.3 E18-E05 N27-N34
      2 1443 1402 (N18W30 L=233)   301     94        0 (1399/1374-1493/1422) -26.6 W27-W33 N16-N19
      3 1679 1177 (N03W44 L=247)   412     64        2 (1653/1148-1706/1206) -90.0 W41-W46 N01-N04
      4 1892 1116 (S02W65 L=268)  1271    161        2 (1854/1074-1932/1158) -26.6 W60-W71 S04-N00
      5  734 1102 (N05E17 L=186)   244     47        1 ( 712/1069- 753/1137)  53.1 E18-E15 N03-N07
      6  398 1055 (N05E40 L=163)   323     76        1 ( 377/1006- 412/1094) -90.0 E41-E38 N02-N07
      7 1832 1009 (S08W57 L=260)   647    118        0 (1797/ 947-1860/1047) -71.6 W53-W59 S11-S05
      8 1771  823 (S19W51 L=254)  2391    453        5 (1664/ 687-1910/1004)  38.0 W42-W66 S26-S08
      9  113  979 (N04E71 L=132)   565    125        1 (  90/ 964- 143/ 999)  -6.3 E76-E66 N02-N04
     10 1766  933 (S12W50 L=253)   450     79        1 (1736/ 899-1793/ 956) -11.3 W47-W52 S14-S10
     11  435  845 (S07E39 L=164)  1248    303        0 ( 347/ 778- 528/ 907)   0.0 E47-E33 S11-S03
     12 1587  789 (S20W35 L=238)   609    161        1 (1521/ 754-1644/ 817)  12.5 W30-W39 S21-S18
     13 1613  734 (S23W38 L=241)   358    111        0 (1573/ 711-1658/ 762) -29.7 W34-W41 S25-S21
     14  532  615 (S22E37 L=166)   259    101        0 ( 509/ 574- 560/ 653)  33.7 E39-E33 S24-S19
     15  853  584 (S27E15 L=188)   454     74        1 ( 814/ 552- 891/ 610) -11.3 E17-E12 S29-S25
     16 1555  406 (S46W40 L=243)   875    133        1 (1504/ 380-1603/ 438) -15.3 W34-W45 S48-S43
  

    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