Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  925 1545 (N28E09 L=224)   252     83        0 ( 905/1497- 948/1578) -51.3 E11-E07 N24-N29
      2  796 1488 (N23E17 L=216)   599    144        0 ( 710/1460- 862/1506)  -5.2 E23-E12 N21-N24
      3  708 1412 (N18E22 L=212)   238     52        0 ( 675/1389- 736/1433) -26.6 E24-E20 N16-N19
      4  597 1362 (N14E29 L=205)   372    109        1 ( 581/1317- 615/1421) -81.9 E30-E27 N11-N18
      5 1897 1220 (N14W64 L=298)  1999    329        1 (1831/1132-1968/1317) -22.8 W57-W77 N09-N19
      6 1207 1264 (N11W10 L=244)   415     82        2 (1164/1225-1251/1308) -38.7 W07-W12 N08-N13
      7 1585 1213 (N10W35 L=269)   511    156        1 (1520/1193-1654/1240)   0.0 W30-W40 N08-N11
      8  401 1181 (N03E41 L=193)   250     49        1 ( 376/1161- 430/1195)  11.3 E43-E38 N01-N03
      9  452 1043 (S05E37 L=197)   751    127        2 ( 430/ 978- 478/1113)  90.0 E38-E34 S09-S01
     10  857 1059 (S03E10 L=224)   176     34        1 ( 832/1046- 881/1080)   0.0 E11-E08 S03-S01
     11 1865  981 (S01W61 L=295)   448     49        1 (1847/ 949-1887/1005)  45.0 W58-W63 S02-N00
     12  156  865 (S16E65 L=169)  2369    333        8 ( 122/ 765- 188/ 970)  71.6 E69-E61 S21-S09
     13 1594  914 (S07W37 L=271)   298     41        1 (1566/ 889-1617/ 936)  18.4 W35-W38 S08-S05
     14  835  888 (S13E11 L=223)   207     38        0 ( 811/ 872- 863/ 903)  33.7 E12-E09 S14-S12
     15 1138  756 (S20W09 L=242)  1178    273        0 (1033/ 659-1289/ 840)  40.9 W02-W17 S26-S13
     16 1783  551 (S27W63 L=296)  1693    283        2 (1699/ 510-1849/ 625)  39.3 W54-W68 S30-S21
     17  915  522 (S36E05 L=229)   490    161        1 ( 836/ 476- 976/ 570)  24.4 E11-W00 S40-S32
  

    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