Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1193 1787 (N52W41 L=247)   358     40        0 (1168/1767-1219/1806)   0.0 W37-W45 N51-N53
      2  897 1769 (N60W10 L=216)   693    114        1 ( 864/1736- 941/1790)  11.3 W05-W15 N57-N62
      3  267 1323 (N36E58 L=148)   566     55        4 ( 234/1303- 298/1343)   0.0 E63-E53 N34-N37
      4  816 1193 (N20E10 L=196)   616    251        1 ( 711/1163- 917/1229) -17.1 E16-E03 N17-N23
      5 1730 1009 (S08W47 L=253)   394     58        1 (1709/ 985-1764/1035) -36.9 W44-W49 S09-S06
      6  191  956 (N13E64 L=142)   430     58        1 ( 167/ 935- 218/ 981) -15.9 E67-E60 N11-N14
      7  956  776 (S07E08 L=198)   303     31        0 ( 926/ 753- 981/ 798) -56.3 E09-E06 S08-S05
      8 1282  735 (S15W11 L=217)   306     94        1 (1250/ 694-1321/ 772)  38.7 W08-W13 S17-S13
      9  948  711 (S11E10 L=196)   254     83        1 ( 923/ 687- 977/ 735)  14.0 E11-E07 S12-S09
     10  245  658 (S06E66 L=140)  2169    283        5 ( 203/ 566- 286/ 731)  70.0 E70-E61 S11-N00
     11 1769  611 (S36W53 L=259)  1291    154        2 (1702/ 578-1829/ 652)   0.0 W43-W62 S38-S34
     12 1183  491 (S30W00 L=206)  1157    300        1 (1110/ 399-1300/ 600)  36.3 E06-W09 S35-S24
     13 1565  559 (S34W30 L=236)   363     30        2 (1544/ 538-1588/ 578) -26.6 W28-W32 S35-S32
     14  651  530 (S18E34 L=172)   716    198        1 ( 581/ 502- 729/ 572) -11.3 E37-E27 S19-S14
     15  858  517 (S22E20 L=186)   330     98        1 ( 833/ 487- 902/ 550)  26.6 E22-E16 S23-S19
     16  669  326 (S33E43 L=163)  3111    324        5 ( 623/ 182- 759/ 439)  24.5 E76-E30 S47-S24
  

    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