Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  757 1760 (N53E31 L=159)   446     73        1 ( 715/1739- 792/1776)  23.2 E35-E27 N51-N54
      2 1586 1679 (N48W58 L=248)   692     63        0 (1552/1657-1620/1698)   7.1 W50-W66 N45-N49
      3 1181 1546 (N37W10 L=201)   362     82        0 (1139/1527-1223/1567)   9.5 W07-W13 N35-N38
      4  413 1489 (N30E50 L=140)  1813    370        2 ( 270/1441- 560/1545)  -5.6 E76-E35 N26-N32
      5 1240 1491 (N33W14 L=205)   284     52        1 (1217/1468-1261/1523)  76.0 W12-W16 N31-N35
      6 1259 1334 (N23W15 L=205)   423     88        2 (1240/1285-1279/1380)  71.6 W13-W16 N19-N25
      7  280 1163 (N09E54 L=137)   836    176        1 ( 258/1097- 309/1256) -51.3 E57-E49 N04-N14
      8 1668 1055 (N06W43 L=234)  1438    331        2 (1582/ 927-1715/1179) -61.9 W36-W47 S02-N13
      9  599 1073 (N05E27 L=164)   295    144        0 ( 570/1026- 665/1110)  39.8 E28-E22 N01-N07
     10 1845  834 (S08W63 L=253)   520     35        0 (1825/ 810-1865/ 866) -45.0 W60-W65 S09-S06
     11 1570  759 (S13W37 L=228)   323     60        0 (1544/ 736-1592/ 784)  33.7 W35-W39 S14-S11
     12  261  700 (S21E59 L=132)   775    127        1 ( 212/ 665- 303/ 739)  21.8 E67-E52 S23-S17
     13  584  640 (S23E29 L=161)   357     96        0 ( 569/ 588- 605/ 686)  66.8 E30-E27 S26-S19
     14 1759  591 (S24W61 L=252)  2723    218        5 (1676/ 525-1817/ 683)  50.7 W53-W68 S28-S17
     15  861  579 (S26E10 L=181)   852    218        1 ( 795/ 544- 943/ 635) -21.8 E14-E04 S28-S21
     16  662  591 (S26E24 L=167)   472    144        0 ( 598/ 564- 717/ 628)  21.8 E29-E19 S27-S22
     17 1276  218 (S55W32 L=222)  2674    235        7 (1169/ 185-1399/ 254)   6.1 W20-W49 S59-S51
  

    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