Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  878 1708 (N38E30 L=143)   321     42        0 ( 853/1687- 908/1731)  -9.5 E33-E27 N36-N39
      2  926 1567 (N29E20 L=153)   639    188        1 ( 829/1530- 991/1607)  45.0 E24-E16 N24-N33
      3 1678 1265 (N27W40 L=213)  2544    757        2 (1470/ 981-1808/1544) -39.6 W20-W55 N12-N41
      4  634 1460 (N16E35 L=138)   272     83        0 ( 595/1427- 666/1492)  71.6 E36-E33 N12-N18
      5  269 1421 (N06E66 L=107)   542     97        0 ( 257/1377- 277/1476) -40.6 E70-E62 N03-N09
      6 1347 1343 (N24W13 L=186)   375    175        0 (1284/1318-1416/1383)   6.3 W08-W17 N22-N26
      7  602 1332 (N07E33 L=140)   221     48        1 ( 574/1321- 635/1348)   0.0 E34-E30 N06-N08
      8 1384 1290 (N22W16 L=190)   214     54        1 (1363/1265-1400/1320)  90.0 W15-W17 N20-N23
      9  270 1219 (S06E57 L=116)   418     61        1 ( 248/1191- 297/1239)  45.0 E58-E53 S07-S04
     10  520 1029 (S12E31 L=142)   547    141        2 ( 477/ 956- 578/1070)  60.3 E34-E27 S16-S09
     11  268  855 (S28E52 L=121)   854    242        0 ( 224/ 780- 299/ 946)  56.3 E58-E48 S33-S21
     12 1768  863 (N05W54 L=227)   332     74        0 (1731/ 846-1791/ 887)  26.6 W50-W56 N03-N06
     13 1004  864 (S12W02 L=175)   397    123        0 ( 950/ 847-1066/ 885)  18.4 E01-W05 S12-S09
     14  740  800 (S21E13 L=160)   239     93        1 ( 698/ 772- 772/ 837)  59.0 E15-E11 S23-S18
     15 1639  584 (S14W53 L=226)  2633    421       10 (1584/ 426-1706/ 738) -66.8 W48-W62 S24-S03
     16 1512  553 (S19W44 L=217)   487    138        1 (1450/ 532-1570/ 585)  32.0 W40-W48 S21-S16
     17  446  532 (S47E36 L=137)   717     84        1 ( 407/ 499- 485/ 572)  30.3 E42-E30 S50-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