Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1601 1684 (N48W61 L=249)   580     73        0 (1579/1655-1624/1705)  17.1 W54-W67 N45-N49
      2 1225 1548 (N37W14 L=202)   353     76        0 (1196/1522-1260/1569)  11.3 W11-W16 N35-N38
      3 1286 1483 (N33W18 L=206)   336    110        0 (1254/1449-1338/1508) -18.4 W15-W21 N30-N34
      4  456 1477 (N29E46 L=142)   596    114        0 ( 427/1428- 502/1507)  68.2 E48-E41 N26-N31
      5  552 1459 (N29E36 L=151)   354     76        1 ( 519/1439- 595/1480)  -8.1 E39-E32 N27-N30
      6 1302 1327 (N22W18 L=205)   294     81        1 (1288/1288-1324/1374)  71.6 W16-W19 N19-N25
      7  313 1156 (N08E50 L=137)   553    131        1 ( 291/1094- 340/1221) -53.1 E53-E47 N04-N12
      8 1671 1147 (N11W44 L=232)   463     99        0 (1628/1111-1703/1180) -39.8 W40-W46 N08-N13
      9 1726  989 (N01W48 L=236)   891    184        2 (1689/ 914-1751/1062) -76.0 W45-W50 S03-N05
     10 1863  803 (S11W66 L=254)   517     79        1 (1838/ 769-1882/ 839)  38.7 W63-W68 S12-S08
     11 1612  755 (S13W41 L=228)   348     33        1 (1585/ 735-1635/ 776)  14.0 W38-W42 S14-S11
     12  288  695 (S21E55 L=132)   584    134        0 ( 235/ 663- 333/ 734)  15.9 E63-E49 S22-S17
     13  901  581 (S25E07 L=180)   855    194        3 ( 837/ 544- 981/ 641) -16.7 E11-E01 S27-S21
     14 1768  574 (S25W64 L=251)  1891    197        2 (1694/ 523-1818/ 634)  37.9 W55-W70 S28-S21
     15  700  584 (S26E21 L=166)   376    130        0 ( 637/ 559- 748/ 622)  24.0 E26-E17 S27-S23
     16 1295  218 (S55W34 L=222)  1836    312        5 (1165/ 189-1408/ 261)  16.3 W19-W48 S58-S50
  

    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