Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1142 1548 (N38W13 L=298)   775    233        1 (1037/1481-1234/1606)  22.4 W04-W21 N33-N41
      2  384 1571 (N43E58 L=227)   545     75        1 ( 344/1550- 425/1590)  11.3 E64-E52 N41-N44
      3  871 1525 (N38E08 L=277)   235     47        1 ( 844/1501- 893/1545) -56.3 E09-E06 N36-N39
      4 1706 1452 (N25W58 L=342)   369     76        1 (1672/1438-1750/1464) -14.0 W54-W62 N23-N26
      5  844 1253 (N20E10 L=275)   299     84        0 ( 806/1224- 887/1282) -38.7 E12-E07 N18-N22
      6 1361 1195 (N13W23 L=307)   417    168        0 (1282/1174-1423/1215)  -6.3 W17-W26 N11-N14
      7  309 1094 (N12E50 L=235)   476    116        1 ( 269/1056- 349/1140) -29.1 E54-E45 N09-N14
      8 1067  999 (N03W02 L=287)   692    265        0 (1019/ 874-1093/1116) -81.9 W00-W03 S04-N10
      9  729 1070 (N09E18 L=267)   271     99        0 ( 709/1029- 767/1095)  36.9 E19-E15 N06-N10
     10 1812  996 (S04W56 L=341)   507    115        0 (1765/ 967-1851/1041)  18.4 W51-W61 S05-S01
     11  213  680 (S14E67 L=217)  1023    106        1 ( 194/ 632- 234/ 741)  66.8 E72-E63 S17-S10
     12  567  647 (S16E33 L=252)   360     53        1 ( 544/ 623- 599/ 670)  26.6 E34-E30 S17-S14
     13  536  599 (S19E37 L=248)   417     66        1 ( 512/ 563- 570/ 636) -53.1 E37-E34 S21-S16
     14  891  531 (S26E12 L=272)   587    224        0 ( 834/ 455- 983/ 592)  36.9 E17-E05 S30-S21
     15 1467  497 (S33W29 L=314)   342     56        1 (1448/ 471-1489/ 521)  71.6 W27-W31 S34-S31
     16 1414  359 (S44W28 L=313)  2740    292       12 (1292/ 288-1528/ 422)  15.9 W17-W39 S49-S37
     17  848  261 (S47E23 L=261)  3522    499        6 ( 648/ 194-1068/ 335)  11.7 E55-E02 S54-S41
  

    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