Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1076 1860 (N55W48 L=270)   682     80        0 (1052/1818-1098/1895)  13.2 W40-W57 N52-N57
      2 1704 1593 (N16W70 L=292)  2055    208        3 (1620/1513-1768/1655) -67.4 W65-W77 N09-N21
      3 1536 1594 (N22W54 L=276)   354     89        0 (1489/1568-1570/1612) -51.3 W50-W55 N19-N24
      4  910 1471 (N35W07 L=229)  1374    445        1 ( 732/1428-1139/1504) -30.3 E06-W19 N25-N40
      5 1646 1423 (N09W52 L=274)   570    117        0 (1599/1382-1673/1471)  20.6 W46-W55 N07-N11
      6 1211 1250 (N13W17 L=239)  1265    514        1 ( 981/1150-1381/1308) -41.8 W03-W24 N03-N20
      7  750 1181 (N22E12 L=210)   714    192        0 ( 663/1152- 845/1213) -14.0 E17-E05 N20-N23
      8  279 1140 (N31E50 L=172)   413     62        1 ( 258/1102- 295/1176) -78.7 E52-E48 N28-N33
      9 1815 1118 (S12W54 L=277)   666     95        2 (1777/1085-1862/1154) -24.0 W50-W60 S15-S09
     10  768 1065 (N15E13 L=209)   507    217        0 ( 689/1027- 853/1088) -16.7 E19-E09 N11-N16
     11  257  892 (N16E56 L=166)   907    205        1 ( 166/ 840- 342/ 930) -21.3 E66-E48 N12-N19
     12 1583  820 (S21W28 L=250)   371    123        0 (1550/ 760-1620/ 861)  23.2 W24-W31 S24-S19
     13  159  722 (N07E77 L=145)  1388    104        4 ( 141/ 677- 184/ 778)  81.9 E82-E72 N04-N11
     14 1692  543 (S45W36 L=258)  3584    597        1 (1451/ 387-1844/ 728)  12.8 W10-W55 S48-S37
     15  772  446 (S21E33 L=189)   677    147        2 ( 750/ 362- 791/ 510)  59.0 E36-E28 S26-S16
     16 1018  244 (S43E29 L=193)   364     38        1 ( 990/ 226-1045/ 264)   0.0 E32-E26 S43-S41
     17 1178  218 (S51E21 L=201)  1390    235        2 (1070/ 179-1280/ 247)  -7.8 E32-E10 S54-S48
  

    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