Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1781 1494 (N11W68 L=  2)  2011    305        7 (1739/1397-1832/1581)  57.5 W62-W72 N05-N16
      2  694 1547 (N42E10 L=284)   375     79        1 ( 664/1510- 729/1578)  18.4 E12-E06 N39-N43
      3  830 1435 (N31E03 L=291)   337    105        0 ( 775/1417- 886/1446) -26.6 E06-W00 N28-N32
      4 1264 1244 (N10W19 L=313)   690    205        5 (1179/1214-1364/1268) -24.4 W13-W24 N06-N11
      5  848 1186 (N16E07 L=287)  1390    359        1 ( 673/1123-1024/1222)  -8.1 E18-W03 N11-N18
      6 1741  977 (S16W45 L=339)   477     82        1 (1723/ 934-1759/1023)  51.3 W42-W46 S18-S13
      7  132  992 (N18E66 L=228)   961    242        1 (  86/ 967- 181/1011) -14.0 E72-E60 N16-N19
      8  376  977 (N13E42 L=252)   379     78        1 ( 342/ 945- 403/1003) -59.0 E44-E39 N10-N15
      9  277  841 (N07E54 L=240)   606    104        2 ( 252/ 811- 307/ 873)  36.9 E56-E50 N05-N09
     10 1231  655 (S23W04 L=298)  1520    473        5 (1112/ 501-1352/ 844)  37.3 E06-W15 S30-S14
     11 1185  746 (S17W03 L=297)   263     51        2 (1159/ 721-1210/ 770) -63.4 W01-W04 S18-S14
     12  530  615 (S11E40 L=254)   871    200        6 ( 457/ 548- 604/ 664) -56.3 E44-E35 S15-S06
     13  423  568 (S11E50 L=244)  1310    230        4 ( 375/ 500- 471/ 623)  45.0 E55-E45 S15-S08
     14  318  530 (S12E62 L=232)  1616    255        5 ( 265/ 467- 364/ 597) -69.4 E66-E58 S15-S07
     15  649  530 (S18E35 L=259)   610    157        1 ( 610/ 476- 692/ 583)  32.0 E39-E31 S20-S14
     16 1068  336 (S41E16 L=278)   540    181        3 (1021/ 293-1097/ 388) -90.0 E18-E13 S44-S35
     17  773  253 (S40E42 L=252)   404     79        2 ( 745/ 221- 806/ 278) -63.4 E43-E39 S42-S37
     18  996  198 (S51E31 L=263)   851    167        2 ( 974/ 159-1026/ 238)  52.1 E35-E26 S55-S46
  

    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